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ABSTRACT 

This  thesis  contains  programs  written  in  APL  and 
documentation  for  performing  certain  nonpar ametr ic 
tests  and  computing  nonpar ametr ic  confidence  intervals. 
These  methods  of  inference  are  particularly  useful  in 
dealing  with  Department  of  Defense  related  problems  as 
illustrated  in  the  several  military  examples  worked  in 
Appendix  C.  The  following  nonpar ametr ic  tests  are 
considered:  Sign  Test,  Wilcoxon  Signed-rank  Test,  Mann- 
Whitney  Test,  Kruskal-Wallis  Test,  Kendall's  B, 
Spearman's  R,  and  Nonpar ametr ic  Linear  Regression.  The 
tests  are  based  on  the  exact  distributions  of  the 
respective  test  statistics  unless  a  large  sample 
approximation  is  determined  to  provide  at  least  a  three 
decimal  place  accuracy.  The  software  consists  of  two 
APL  workspaces;  one,  which  is  designed  for 
microcomputers  (IBM  PC's  or  compatibles)  and  is  menu- 
driven,  and  the  other,  without  menus,  is  designed  for 
the  mainframe  computer  (IBM  3033)  at  the  Naval 
Postgraduate  School. 
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I . INTRODUCTION 

Although  nonpar ametric  procedures  are  powerful 
tools  to  the  analyst,  they  are  currently  underused  and 
often  avoided  by  potential  users.  Perhaps  one  reason 
for  this  is  the  difficulty  in  generating  the  exact 
distributions  of  the  test  statistics,  even  for  moderate 
sample  sizes.  Consequently,  tables  of  these 
dis^r icuxions  are  only  available  for  very  small  sample 
sizes  and  normal  theory  based  approximations  must  then 
be  used. 

The  purpose  of  this  thesis  is  to  make  a  variety  of 
nonpar ametric  procedures  quick,  easy  and  accurate  to 
apply  using  menu  driven  computer  programs  in  APL.^ 
These  programs  use  enumeration,  recursion,  or 
combinatorial  formulas  to  generate  the  exact  null 
distribution  of  the  various  nonpar ametric  test 
statistics.  This  allows  hypothesis  testing  and 
confidence  interval  estimation  to  be  based  on  exact 
distributions  without  the  use  of  tables.  For  larger 
sample  sizes,   the  normal,   F,   and  T  distributions  are 


"•  APL  was  chosen  because  it  is  an  interactive 
language  that  is  especially  powerful  at  performing 
calculations  dealing  with  rank  order  statistics  and 
vector  arithmetic.  Menus  are  not  included  in  the 
workspace  designed  for  the  mainframe. 


used   to   approximate   the  distributions   of   the   test 
statistics  with  three  decimal  place  accuracy. 

Section  II  addresses  workspace  design  issues,  to 
include,  workspace  requirements  and  assumptions 
regarding  its  use.  Section  III  discusses  the  methods 
used  to  assess  the  accuracy  of  different  asymptotic 
approximations,  and  the  sample  sizes  required  for  an 
approximation  to  yield  three  decimal  place  accuracy. 
Section  IV  gives  background  information  and  discusses 
programming  methodology  for  nonpar ametr ic  tests  based 
on  single  and  paired  sample  data.  In  Section  V, 
nonpar ametr ic  xests  for  two  or  more  independent  samples 
are  considered.  Section  VI  discusses  nonpar ametr ic 
tests  for  association;  and,  Section  VII  deals  with 
nonparametr ic  simple  linear  regression.  Section  VIII 
recommends  other  nonparametr ic  tests  that  may  be  added 
to  the  workspace  and  areas  for  further  work. 

To  show  application  of  nonparametr ic  statistical 
methods  to  Department  of  Defense  problems,  several 
military  examples  are  worked  in  Appendix  C. 
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II .WORKSPACE  DESIGN  ISSUSS 

This  section  presents  a  brief  overview  of  the 
design  considerations  used  in  developing  the  APL 
workspace  for  both   the  mainframe  and  microcomputer. 

A.  EQUIPMENT  AND  SOFTWARE  REQUIREMENTS 

The  microcomputer  must  be  an  IBM  PC  or  AT 
compatible,  equipped  with  512  kilobytes  of  RAM  and  the 
APL* PLUS/ PC  system  software,  release  3.0  or  later,  and 
IBM's  DOS,  version  2.00  or  later.  ">  The  8087  math 
coprocessor  chip  is  not  required  to  run  this  software, 
but  will  increase  the  computational  speed. 

B.  KNOWLEDGE  LEVEL  OF  THE  USER 

The  user  is  expected  to  have  had  some  exposure  to 
APL  and  a  working  knowledge  of  nonparametr ic 
statistics.  Fsuniliarity  with  microcomputers  or  the 
Naval  Postgraduate  School  mainframe  computer  is 
assumed. 


"I  The  APL  system  software  requires  14-4  kilobytes  of 
RAM  while  the  NONPAR  workspace  requires  an  additional 
190  kilobytes. 
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C.  SELECTION  OF  TESTS 

The  nonpar ame trie  tests  chosen  for  this  workspace 
are  some  of  the  more  widely  known,  and  are  considered 
basic  material  for  any  nonpar ametric  statistics  course. 
More  information  about  the  tests  can  be  found  in  any  of 
the  textbooks  that  are  referred'  to  in  this  document. 

D.  MENU  DISPLAYS 

The  microcomputer's  workspace  is  designed  around 
the  use  of  menus.  This  was  accomplished  using  the 
software  package  ?C  TOOLS  from  STSC .  These  menus  are 
designed  to  guide  the  user  through  the  selection  of  the 
tests  without  an  excessive  amount  of  prompting.  The 
main  menu  displays  the  choices  available  in  the 
workspace,  while  the  test  menus  give  the  background 
information  and  options  available  for  each  test-  Help 
menus  to  provide  additional  information  about  the 
tests  are  also  available. 

E.  ORGANIZATION  OF  WORKSPACE  DOCUMENTATION 
Separate    documentation   is   included    for    the 

microcomputer's  and  mainframe  computer's  workspaces 
(see  Appendices  A  and  B,  respectively).  These 
appendices  explain  the  organization  and  operation  of 
the  workspaces.  Appendix  C,  which  provides  example 
problems  for  each  nonpar ametric  test,  is  applicable  to 
both  workspaces. 
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Ill .GENERAL  SAMPLE  SIZE  CONSIDERATIONS  AND  ASYMPTOTIC 
APPROXIMATIONS 

In  this  thesis,  the  term  alpha  value  is  used  in 
a  general  sense,  and  refers  to  the  probability  of 
rejecting  a  true  null  hypothesis.  The  term  ?-value 
refers  xo  the  probability  that  a  test  statistic  vill 
exceed  (or  not  exceed  in  the  lower-tailed  test)  the 
computed  value,  when  the  hypothesis  being  tested  is 
~rue  . 

For  selected  values,  the  exact  cumulative 
distribution  functions  (C.D.F.)  of  the  test  statistics 
are  compared  with  those  obtained  from  normal  based 
asymptotic  approximations.  The  results  of  the  compari- 
sons are  used  as  a  basis  for  assessing  the  accuracy  of 
the  approximations.  In  those  cases  where  more  than  one 
asymptotic  approximation  has  been  suggested  in  the 
literature,  the  accuracy  of  each  approximation  is 
compared  over  a  range  of  desired  C.D.F.  values  and 
sample  sizes.  From  the  results,  the  most  consistently, 
accurate  approximation,  and  the  sample  size  for  which 
"hat  approximation  provides  at:  least  three  decimal 
place  accuracy  is  determined. 

Once  the  accuracy  comparisons  were  completed  for  a 
specific  nonpar ametric  test,  microcomputer  capabilities 


13 


were  considered-  In  some  cases,  generation  of  the 
exact  distribution  up  to  the  desired  sample  size  took 
too  long  or  was  not  possible  on  the  PC.  When  this 
occurred,  the  mainframe  computer  was  used  to  generate 
the  required  distributions  with  the  results  stored  in 
numerical  matrices  for  quick  recall  by  the 
nonpar ametr  ic  tesz    prograjns. 
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IV. TESTS  ?0R  LOCATION  BASED  ON  SINGLE  AND  PAIRED- 
SAMPLE  DATA 

The  tests  assume  that  the  data  consists  of  a 
single  set  of  independent  observations  Xj.  or  paired 
observations  (Xi_,Yj_),  i  =  1,2,...,N,  from  a  continuous 
distribution.  For  the  single  and  paired-sa.Taple  cases, 
the  null  hypotheses  are  concerned  with  the  median  of 
the  Xj_  and  the  median  of  the  differences  Xj_-  Yj_, 
respectively.  The  tests  considered  are  the  Ordinary 
Sign  Test  and  zhe    Wilcoxon  Signed-Rank  Test. 

A.    ORDINARY  SIGN  TEST 

The  Sign  test  can  be  used  to  test  various 
hypothesis  about  the  population  median  (or  the  median 
of  the  population  of  differences).  Confidence 
intervals  for  these  parameters  can  also  be  constructed. 
As  a  final  option,  nonparajnetr ic  confidence  intervals 
for  the  quantiles  of  a  continuous  distribution  are 
offered. 

1 .    Computation  of  the  Test  Statistic 

For  single-sample  data,  the  test  statistic  K 
is  computed  as  the  number  of  observations  Xi  greater 
than  the  hypothesized  median  Mq .  For  the  paired-sample 
case,  K  is  the  number  of  differences  Xj^-  Y^  that  exceed 
M0 .    All  observations  X^  (or  Xj^-  Y±)    that  are  equal  to 
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M0  are  ignored  and  the  sample  size  decreased 
accordingly.  As  long  as  the  number  of  such  ties  is 
small  relative  to  the  size  of  the  sample,  the  test 
results  are  not  greatly  affected.  Gibbons  [Ref.  2:pp. 
1  08]  . 

2 .    The  Null  and  Asymptotic  Distribution  of  K 

The  null  distribution  of  K  is  binomial 
with  p  =  .5.  In  Table  1 ,  the  exact  values  of  the 
C.D.F  are  compared  with  the  corresponding  approximate 
values  using  a  normal  approximation  with  and  without 
continuity  correction. 

TABLE  1  .   C.D.F.  COMPAx^ISONS  FOR  THE  SIGN  TEST 


PROBCX  1  kJ;  FOR  SAMPLZ  3IZZ  EQUAL  TO  24. 
-I 1 1 1 1 


TEST  STAT.  VALUE 

5 

6 

7 

8 

9    1 

10   1 

1 

11 

IXACT  CD.?. 

.00331 

.01133 

.03196 

.07579  1 

.15373  1 

1 

1 

.27063  i 

.41941 

ER.50R;  NORMAL 

.00117 

.00417 

.01134 

.02456  1 

.04339  1 

I 

.06352  1 
1 

.07736 

ERROR;  NORM.  U/CC 

".00063 

".00104 

".00114 

".00073  I 

1 
.00001  1 

1 
.00043  1 

.00023 

1 
PROBCK  < 

.1  . 
k] 

;  FOR 

1        1 
SAMPLE  SIZE 

1 
EQUAL  TO 

1 
25. 

1 

TEST  STAT.  VALUE 

1    6 

1 

1 
1 

7 

1        1 
1    8   1 

1        1 

1 
9   1 

1 

1 
10   1 

1 

1 
11    1 

1 

12 

EXACT  C.D.F. 

1  .00732 

1 

-I — 
1 

.02164 

1        1 

1  .05333  1 
1        1 

.11476  1 

1 

1 
.21213  1 

1 
.34502  1 

. 50000 

iSROS;  MORMAL 

1  .00266 

1 

.00774 

'     ■    1 
!  .31795  ! 

. J3400  i 

,-,  q  ■■;  c  ^  1 

.07077  i 

.07926 

ZRRQ2;  :<CRM,  'J/ CC 

1  '.00033 

1 

1 

.00111 

1  ".00092  1 
1        1 

'.00031  1 
1 

.00032  1 

1 

.00044  i 
1 

. 00000 
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As  can  be  seen,  for  sample  sizes  greater  than 
25,  a  normal  approximation  with  continuity  correction 
is  accurate  to  at  least  three  decimal  places. 

3 .  Hypothesis  Testing 

P-values  are  computed  for  three  basic 
hypotheses  comparing  the  median  of  the  population  or 
the  median  of  the  population  of  differences  M  with  a 
hypothesized  median  Mg .  P-values  are  taken  from  the 
cumulative  distribution  of  the  binomial  for  the 
following  tesTis  of  hypoTihesis. 

a.  One-sided  Tests 

(1)  H0 :  M  =  Mg  Versus  HI:  M  <  Mg .  The 
P-value  equals  Pr  [K  _<_  k]  ,  where  k  is  the  computed 
value  of  the  test  statistic. 

(2)  H0:  M  =  M0  Versus  HI:  M  >  Mq .  The 
P-value  equals  Pr  [K  2.  k]  • 

b.  Two-sided  Test. 

(1)  H0:  M  =  M0  Versus  HI:  M  ^  Mq .  The 
P-value  equals  twice  the  smaller  value  of  a( 1 )  or 
a(2),  but  does  not  exceed  the  value  one. 

For  sample  sizes  greater  than  25,  a  normal 
approximation  with  continui*;y  correction  is  used. 

4 .  Confidence  Interval  Sstimation 

Confidence  intervals  for  the  population 
median  are  based  on  the  ordered  observations  in  the 
sample.    For  paired-sample  data,  confidence  bounds  are 
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obtained  from  the  ordered  differences  of  the  pairs  of 
data.  A  100(1-  ct  )%  confidence  interval  is  determined 
in  the  following  manner.  Let  k  be  the  number  such  that 
Pr[K  1  i^]  1  (  ^/2).  Then,  the  (k+1)th  and  (N-k)th 
order  statistics  constitute  the  end  points  of  the 
confidence  interval.  Gibbons  [Ref.  1:  pp.104-]. 

For  computing  confidence  intervals  when 
sample  size  N  is  greater  than  25,  a  normal 
approximation  with  continuity  correction  is  used. 

Also  included  under  this  test  is  an  option  to 
generate  nonparametr ic  confidence  intervals  for  any 
specified  quantile  given  a  sample  size  N  from  a 
continuous  distribution.  The  end  points  of  the 
intervals   are  sample  order  statistics. 

B.    WILCOXON  SIGNED-RANK  TEST 

The  signed-rank  test  requires  the  added  assumption 
that  the  underlying  distribution  is  symmetric.  This 
test  uses  the  ranks  of  the  differences  Xi-  M0  (or  Xi- 
Yi-  M0 )  together  with  the  signs  of  these  differences  to 
determine  the  test  statistic.  Confidence  intervals  for 
the  median  can  also  be  constructed. 

1 .    Computation  of  the  Test  Statistic 

For  single-sample  data,  the  test  statistic  W 
is  computed  as  follows. 
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(  1   if  Xi-  M0  >  0 
Let  Zj_  = 

(  0   if  Xi-  Mg  <_  0 

N 
and  let  ri  =  rank(|Xj_-  Mq  |  )  .   Then,  W  =  ^Z±r±, 

i=1 

For  paired  sample  data,  W  is  calculated  in 
the  same  manner ,  except  the  differences  to  be  ranked 
are  the  paired-differences  minus  the  hypothesized 
median.  Zero  differences  are  ignored  and  the  sample 
size  is  decreased  accordingly.  When  ties  occur  between 
ranks,  the  average  value  of  the  ranks  involved  are 
assigned  to  the  tied  positions.  It  has  been  shown  that 
a  moderate  number  of  ties  and  zero  differences  has 
little  effect  on  the  test  results.^ 

2 .    The  Null  and  Asymptotic  Distribution  of  W 

The  exact  null  distribution  of  W  is  given  by: 
Pr[W  =  w]  =  un(w)/2N,  w  =  0 , 1 , 2 , . . . ,N( N+1 ) /2 ,  where 
un(w)  is  the  number  of  ways  to  assign  plus  and  minus 
signs  to  the  first  N  integers  such  that  the  sum  of  the 
positive  integers  equals  w.  It  can  be  shown  (see 
Gibbons  [Ref.  1:pp-  112])  that  un(w),  for  successive 
values  of  N,  can  be  computed  using  the  recursive 
relationship : 

un(w)  =  UN-i('*^-N)  +  UN_i(w) 


iFor  more  information  on  the  effects  that  zeros  and 
tied  ranks  have  on  the  Wilcoxon  Signed-Rank  Test,  see 
Pratt  and  Gibbons  [Ref.  3]. 
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Exact  C.D.F.  values  were  compared  with  those 
obtained  using  the  following  asymtotic  approximations: 
student's  T  with  (N-1)  degrees  of  freedom  (T), 
student's  T  with  continuity  correction  (TC),  normal 
(Z),  normal  with  a  continuity  correction  (ZC),  the 
average  of  T  and  Z  as  suggested  by  Iman  [Ref.  4],  and 
the   average  of  TC  and  ZC. 

As  can  be  seen  in  Table  2  below,  the  average 
of  TC  and  ZC  gives  the  most  consistently  accurate 
results  with  three  decimal  place  accuracy  when  the 
sample  size  exceeds  9. 


TABLE 


C.D.F.        COMPARISONS    FOR       THE       WILCOXON 
SIGNED-RANK    TEST 


PROBCy  i  u»];    FOR  SAMPLE  SIZE  EQUAL  TO     9. 


TEST  STAT.  VALUE 

r — 1 

3 

5 

r -■   I 

6    1 

3 

1 ■ 

9 

12 

1    14 

EXACT 

C.D.F.  . 

.00977 

.01953 

1 

.02734  1 

1 

.04383 

.06445 

.12500 

.17969 

ERROR; 

NORMAL 

-.00067 

.00046 

.00204  1 

1 

,00591 

.00958 

.01324 

.02272 

ERROR; 

NORM,  y/cc 

■.00243 

-.00247 

1 
-.00167  1 

1 

.00023 

.00269 

.00693 

.00306 

ERROR; 

T  DI3T 

.00513 

.00608 

.00673  1 

1 

.00701 

.00317 

.00326 

.00318 

ERROR; 

T  U.'CC 

.00355  t 

.00276 

1 
.00233  1 

1 

.00012 

-.00010 

-.00437 

".00725 

E.^RuR; 

AVE  T/Z 

.00225  1 

.00327 

.00433  1 

.00646 

.00833 

.01325 

.01545 

ZREuR; 

AVE  :c/zc 

,00056  1 

.00014 

.00033  ! 

.00017 

.00123 

.  ^v'OiuS 

.00041 
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TABLE  2.   (Continued) 


PPOBCy  i  w];  FGR  SAMPL2  SIZI  EQUAL  TO  10. 


TEST  3 

TAT.  i^ALUI 

1    5 

7 

T r 

1    8   1 

10 

12 

1 1 

15 

1 

17 

EXACT 

Z.j.T. 

1  .0097? 

.01355 

i  .02441  I 

.04199 

.06543 

r— — —  ■ 
.11621 

.16113 

E?.?G?. ; 

NORMAL 

i  ".00115 

.00023 

1  .00099  1 

.0047b 

.00337 

.01490 

.01333 

E2R0K; 

NOPM.  U/CC 

1  ".00270 

■.00219- 

1  ■.00193  1 

.00043 

.00229 

.00558 

.00710 

iPFOS: 

T  DI3T 

!  .00399 

.00512 

1  ,00525  1 

.00665 

.00661 

.00661 

.00631 

£2?.0R; 

T  U/CC 

1  .00249 

[ 

.00243 

1  .00132  1 

.00151  1 

■.00052 

".00376  1 

".00571 

E2HQ2; 

AUE  T/Z 

1 

!  .00142 

.00267 

1  .00312  ! 

.00571  1 

.00749 

.01075 

.01284 

IP.P.GS: 

:)tjE  TC/ZC 

:  ".00010 

i_  — 

.00012 

1  -.00008  1 

.00097 
.1 

.00039 

.00091 
-. i 

.00070 

3 .    Hypothesis  Testing 

P-values  are  computed  for  three  basic 
hypotheses  comparing  the  median  of  the  population  or 
the  median  of  the  population  of  differences  M  with  a 
hypothesized  mediaji  Mg  as  shown  below. 

a.  One-sided  Tests 

(1)  H0:  M  =  M0  Versus  HI:  M  <  Mg .  The 
P-value  equals  Pr  [W  <_  w]  ,  where  w  is  the  computed 
value  of  the  test  statistic  W. 

(2)  H0:  M  =  M0  Versus  HI:  M  >  Mg .  The 
P-value  equals  Pr  [W  2.  ^1  • 

b.  Tvo-sided  Test 

(1)  H0:  M  =  M0  Versus  HI:  M  ?t  Mg .  The 
P-value  equals  twice  the  smaller  value  of  a( 1 )  or 
a(2),  but  not  exceeding  the  value  one. 
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For  sample  sizes  greater  than  9,  an  average 
of  the  normal  and  student's  T  approximations,  each  with 
continuity  correction,  is  used.  Computations  of  the  P- 
value  for  each  of  the  alternative  hypotheses  are: 

a.  HI  :  M  <  Mg 

Let  Pzc  =  P^  [  Z  _<_  (w  +.5  -  Um)  /    ^w]  S-^d 


let   Ptc  =  P^ 


T(N-1 )  1  _ 


w  -  /; 


w 


-  .5 


Ncrv2  [|w  -  ^^\     -     .5]2 


N  -  1 


.5 


where  Z  is  standard  normal,  T(n-i)  has  a  student's  T 
distribution  with  (N-1)  degrees  of  freedom,  a/^  = 
N(N+1  )/4  and  <J^'2-  =  (  N(  N+1  )  (  2N+1  ) /24- )  )  .  Then,  the 
P-value  for  the  test  is  ( Pzc  +  ( 1  -  Ptc))/2  if  w  is 
less  than  ^^^  and  (Pzc  +  Ptc)/2.  otherwise-  The  above 
formulas  are  obtained  from  those  given  by  Iman  [Ref .  4] 
after  inclusion  of  a  continuity  correction. 

b.  HI :  M  >  M0 

The  P-value  equals  ((1  -  Pzc)  +  Ptc)/2  if 
w  is  less  than  ^^  and  ((1  -  Pzc)  +  (  1  -  Ptc))/2, 
otherwise.  The  computation  of  Pzc  a-^d  Ptc  is  similar 
to  the  above  except  the  sign  of  the  continuity 
correction  is  changed. 

c .  H 1  :  M  ?t  M0 

The  P-value  equals  twice  the  smaller 
value  of  a  or  b  above,  but  not  exceeding  the  value  one. 
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4 .    Confidence  Interval  Estimation 

For  single-sample  data,  the  confidence 
interval  for  the  population  median  is  based  on  the 
ordered  averages  of  all  pairs  of  observations  (X2.-r 
Xi  )/2  such  that  i  <  i.  A  100(1-  a  )?&  confidence 
interval  is  determined  in  the  following  manner.  Let  w 
be  the  number  such  -chat  Pr  [W  _<_  w]  <_  (  a/2).  Then,  the 
(w+1 )th  and  (m-w)th  order  statistics,  where  m 
N(N+1 )/2  or  the  total  number  of  paired-averages, 
constitute  the  end  points  of  the  confidence  interval. 
A  confidence  interval  for  paired-sample  data  is 
computed  in  the  same  manner ,  except  the  end  points  are 
taken  from  the  paired-averages  of  the  differences  Xi- 
Yi .  Gibbons  [Ref.  1:pp.  114-118]. 

For  computing  confidence  intervals  when 
sample  sizes  are  greater  than  9,  a  normal  approximation 
with  continuity  correction  is  used. 
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V. TESTS  BASED  ON  TWO  OR  MORE  SAMPLES 

The  tests  assume  that  the  data  consists  of 
independent  random  samples  from  two  or  more  continuous 
distributions.  The  general  null  hypothesis  is  that  the 
samples  are  drawn  from  identical  populations.  The 
Mann- Whitney  and  Kruskal-Wallis  tests  are  considered. 

A.    MANN-WHITNEY  TEST 

The  Mann-Whitney  test  is  based  on  the  distribution 
of  the  test  statistic  U,  which  can  be  used  to  compare 
the  equality  of  the  population  medians  or  variances 
for  two  samples.^  The  Mann-whitney  test  with  a 
modified  ranking  scheme  can  be  used  to  test  for 
equality  of  variances  if  the  population  means  or 
medians  are  assumed  to  be  equal  (Conover  [Ref.  5:pp. 
229-230]  )  .  If  the  medians  differ  by  a  known  amount ,  the 
data  can  be  adjusted  before  applying  the  test.  A 
confidence  interval  for  the  difference  in  the  medians 
of  the  two  populations  can  also  be  estimated. 

.1  .    Computation  of  the  Test  Statistic 

For  the  comparison  of  population  medians,  the 
test  statistic  U  is  computed  from  the  combined   ordered 


"I  The   test   statistic   U  and   the    method   used   to 
compute  it  are  taken  from  Gibbons  [Ref.  1:pp.  140-14-1]. 
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arrangement  of  observations  Xi  and  Y-j  ,    i  =  1,2,...,N; 

j   =   1,2,...,M.    Let   r-j_  =  rank(Xj_)  in   the   combined 

N 
ordered   sample   and   Rx  =  S  ^i-  Then, 

i  =  1 

U  =  Rx  -  M(M+1 )/2. 

For  testing  the  equality  of  variances,  the 
computation  of  U  is  similar  except  for  the  method  of 
assigning  ranks  to  the  ordered  sample.  This  method 
ranks  the  smallest  value  1,  largest  value  2,  second 
largest  value  3,  second  smallest  value  4.  and  so  on  , 
by  two's,  until  the  middle  of  the  combined  ordered 
sample  is  reached.  For  either  test,  tied  ranks  for  the 
combined  sample  are  assigned  the  average  value  of  the 
ranks  involved.  A  moderate  number  of  ties  has  little 
effect  on  the  test  results  . 

2 .    The  Null  and  Asymptotic  Distribution  of  U 

The  exact  null  distribution  of  U  is 
determined  using  a  recursion  algorithm  due  to  Harding 
[Ref.  6]. 

Exact  C.D.F.  values  were  compared  with 
approximate  values  obtained  from  the  following 
asymtotic  distributions:  student's  T  with  (n-2)  degrses 
of  freedom  where  n  =  N  +  iVI,  the  total  number  of 
observations  in  both  samples  (T),  student's  T  with 
continuity  correction  (TC),  normal  (Z),  normal  with 
continuity   correction   (ZC),   the  average  of  T   and   Z 
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( Iman  [Ref.   7]),   and  the  average  of  TC  and   ZC.    The 
results  for  various  sample  sizes  are  given  in  Table  3. 


TABLE  3.    C.D.F.  COMPARISONS  FOR  THE  MANN-WHITNEY 
TEST 


PROBCa  i 

u];  702 

SAMPLE  3IZI3  N  EQUAL  TO  9 

AND  M  EQUAL 

TO  9. 

TI3T  ; 

jTat.  value 

1 

1   14 

1   17 

1   13 

1   21 

1   23   1 

1 
27   1 

29 

IXACT 

Z.D.J. 

1  .00933 

1  .01999 

1  .02515 

1  .04696 

1  .06796-  i 

.12904  1 

.17005 

12202; 

NORMAL 

1  ".00026 

1  .00100 

1  .00163 

1  .00441 

1  .00632  1 

•               1 

.01243  1 

.01511 

E220R; 

H02M.  lil-'CC 

1  ".00146 

1  -.00114 

1 -.00033 

1 

1  .00026 

1  .00130  t 

.00354  1 

.00436 

I22G2; 

:  3137 

i  .00237 

1  .00337 

1  .00372 

1  .00464 

1  .00525  1 

.00676  1 

.00779 

E22Q2; 

t  y/cc 

i  .00119 

1  .00108 

1  .00095 

1  ,00007 

1  -.00073  1  - 

.00259  i 

-.00323 

12202; 

At;s  T/Z 

1  .00105 

1  .00219 

1  .00270 

1  .00453 

i  ,00603  ! 

.00959  1 

.01145 

Z22Q2; 

AVE  TC/ZC 

1  ".00014 

1  -.00003 

1  .00004 

1  .00017 

f 

1  .00026  1 

1        1 

.00043  I 

t 

.00054 

P20BCU  i 

u];  EOR 

SAMPLE  SIZES  N  EQUAL  TO  7 

AMD  M  EQUAL 

TO  12. 

TZ3T  S 

TAT.  VALUE 

1 

1    14 

1 

1   17 

19 

1   21 

i   24   1 

1        t 

27   i 

30 

EXACT 

CD. 7. 

.00353 

p 
.01792 

p 
.02732 

1 — 
.04156 

r 1 — 

1  .07111  1 
1        1 

.11342  1 

.17012 

E2R02; 

NORMAL 

-.00045 

.00062 

.00137 

.00359 

1  .00701  1 

1        1 

.01097  1 

.01437 

E2E0R; 

NORM.  U/CC 

-.00152 

-.00123 

-.00079 

-.00003 

1  .00154  1 

1        I 

.00322  1 

.00453 

ERROR; 

T  DIST 

.00199 

.00292 

.00356 

.00422 

1  .00527  1 
1        1 

.00643  1 

.00796 

ERROR; 

T  y/cc 

.00094 

.00092 

.00063 

.00026 

1        1 
1  -.00066  1  - 

1        1 

.00173  1 

-.00262 

I22QR: 

AVI  T/Z    ! 

.00077 

.00177  1 

.00272 

.00391 

1        1 
1  .00614  1 

00370  1 

.01142 

Z  ?.  2  U  2 : 

^vi  70/ z: 

-.  .0029  i 

•.J0013  , 

".00005  1 

.00011 

i  . J0044  ; 
i \ — 

00072 

.  J00'?3 
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TABLE    3 


( Continued) 


P2Q3CU  i  ui;    F02  3AMPLI  3IZI3*  M  iQUAL  TO     5  AND   M  SQ.UAL  TO   17. 


i;al'js 


1 1 1 

13   !   20   I   23   I   27 


IXACT 

CO.?. 

1  .00963 

i  .01933 

i  .02916  1 

1 

.04245  1 

-.07013 

i  .12440 

.17935 

£2202; 

M02fiAL 

i  ".00077 

i  .00039 

1  .00170  1 

.00349  1 

.00639 

t  -  — 

.01210 

.01364 

12202; 

M02H.  WCC 

■;  -.00190 

i  ".00130 

-.00037  ! 

r 

.00007  1 

,00133 

. 00444 

\— — 
.00573 

12202; 

T  0I3T 

1  .00130 

1  .00222 

.00299  i 

.00390  ! 

1 

.00344 

.00774 

.00953 

E2202; 

T  U/CC 

1  .00017 

1  .00023 

.00023  1 

1 

.00021  1 

I 

.00007 

".00025  1 

! 

".00051 

12202; 

,4LZ  '/Z 

•  .00027 

.00130 

.00233  ; 

1 
•  .00363  ! 

.00616 

.00992 

.01260 

12202; 

Aui  :■:/:: 

,  ".00036 

".00063 

",00032  i 

.00014  1 

.00097 

.00209  1 

.00264 

P20BCU  i 

uJ:  r02 

3AMPLZ  31 

L          .1 

ZI3  M  IQUA 

L  TO  3 

1        1 
AND  n  EQUAI 

TO  27. 

TI3T  ■; 

TAT.  VALUE 

i    7 

!   10 

1    12   I 

15 

1    19    1 

1 

23   1 

25 

EXACT 

c.]>.y. 

1  .00764 

1  .01650 

1  .02512  1 

1       1 

.04236 

1  .07734  1 
1        1 

.12660  1 

1 

.17433 

EB2Q2; 

N02MAL 

1  ".00265 

1 -.00099 

1  .00072  1 

1 

.00339 

1        1 

1  .00374  1 

1        I 

1 
.01342  1 

.01630 

E2202; 

NORM,  y/cc 

".00363 

1 -.00254 

1 -.00133  1 

1       1 

.00089 

1        I 
1  .00405  1 

I        1 

i 
.00665  I 

1 

.00331 

E2202; 

T  DI3T 

".00121 

1  .00031 

1  .00173  1 

I       1 

.00414 

1        1 
1  .00741  1 

1        1 

1 
.01026  1 

1 

.01253 

E2202; 

T  WZZ 

".00219 

1 -.00129 

1       1 
1 -.00042  1 

1 

.00097 

1        1 
1  .00250  1 

1       1 

.00323  1 

1 

.00390 

Z22Q2; 

Aiv'I  T/Z 

".00193 

1  -.00034 

1 
.00122  1 

I 

.00402 

1  .00303  1 

1        1 

1 
.01134  1 

1 

.01466 

E2R02; 

- —  ■  t 

AVE  TC/ZC   1 

1 

".00291 

i  -.00192 

1 
-.00037  1 

.00093 

1        1 
1  .00327  ! 

'        1 

.00496  1 

1 

.00610 

As  can  "be  seen  from  the  tables.  the  average 
oi  ZC  and  TC  gives  xhe  mosi:  consistently  accurate 
results.  For  sample  sizes  NxM  >  80,  nearly  three 
decimal  place  accuracy  is  obtained  in  all  cases. 
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3 •    Hypothesis  Testing 

P-values  are  computed  for  three  basic 
hypotheses  comparing  the  medians  or  variances  of  the 
two   populations  as  shown  below. 

a.  One-sided  Tests 

( 1 )  H0 :  Mx  =  My  Versus  HI;  Mx  <  My  or 
H0 :  Vx  =  Vy  Versus  HI :  Vx  >  Vy .  The  P-value  equals 
Pr[U  _<_  u]  ,  where  u  is  the  observed  value  of  the  test 
statistic. 

(2)  H0  :  Mx  =  ^^Y  Versus  HI  :  Mx  >  My  or 
H0 ;  Vx  =  Vy  Versus  HI :  Vx  <  Vy .  The  P-value  equals 
?r  [U  2.  u]  . 

b-  Two-sided  Test 

(  1  )  H0 :  Mx  =  My  Versus  HI  :  Mx  9^  My  or 
H0 :  Vx  =  Vy  Versus  HI  :  Vx  i"  Vy .  The  P-value  equals 
twice  the  smaller  value  of  a( 1 )  or  a(2),  but  not 
exceeding  the  value  one. 

For  ssimple  sizes  NxM  greater  than  80,  the 
average  of  the  normal  and  student's  T  approximations, 
each  with  continuity  correction,  is  used.  Computations 
of  the  P-value  for  each  alternative  hypothesis  are: 
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a.  HI  :   Mx  <  My  or  Vx  >  Vy 

Let  Pzc  =  ?!•[  Z  £  (u  +.5  -  ;iu)/  ^u]  a^d. 


let  PTif  =  ?r 


T(n-2)  1 


|u    - 

/^u      -    •  5 

(N+M-1  )'^u- 

.5 

N+iM-2 

N+M-2 

. 

* 

where  Z  is  standard  normal,  ^(^-2)  ^^.s  a  student's  T 
distr ibuxion  with  in-2)  degrees  of  freedom,  /i^  =  NxM/2 
and  (7^2  =  (N(M)(N+M+1  )  )/12.  Then  the  P-value  for  the 
test  is  (Pzc  +  (1  -  Ptc))/2  for  u  less  than  /i^  and 
(Pzc  +  Ptc)/2,  otherwise.  The  above  formulas  are 
obtained  from  those  given  by  Iman  [Ref .  7]  after 
inclusion  of  the  continuity  correction. 

b.  HI :  Mx  >  My  or  Vx  <  Vy 

The  P-value  equals  ((1  -  Pzc)  +  Ptc)/2  if 
u  is  less  than  /y^  and  ((1  -  Pzc)  +  (1  -  Ptc))/2, 
otherwise.  The  computation  of  Pzc  ^^^  ^TC  is  similar 
to  the  above  except  the  sign  of  the  continuity 
correction  is  chainged. 

c  .  HI  :  Mx  ?t  My  or  Vx  ?t  Vy 

The   P-value   equals   twice   the    smaller 
value  of  a  or  b  above,  but  not  exceeding  the  value  one. 
4- .    Confidence  In~erval  Ss'imation 

Confidence  intervals  for  the  difference  in 
medians,  (My  -  Mx ) ,  are  based  on  the  ordered 
arrangement  of  the  differences  (Yj-  Xi ) ,  j  =  1,2,...,M; 
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i  =  1,2,...,N  for  all  i  and  j.  A  100(1-a)% 
confidence  interval  is  determined  in  the  following 
manner .  Let  u  be  the  number  such  that  Pr  [U  _<_  u]  <_ 
(  a     12).  Then,    the   (u+1)th    and   (m-u)th    order 

statistics,  where  m  =  NxM  or  the  total  number  of 
possible  differences,  constitute  the  end  points  of  the 
confidence  interval . 

For  computing  confidence  intervals  when 
sample  sizes  NxM  are  greater  than  80,  a  normal 
approximation  with  continuity  correction  is  used. 

3.    KRUSKAL-WALLIS  TEST 

The  Kruskal-Wallis  test  is  a  nonparametr ic  analog 
of  the  one-way  classification  analysis  of  variance  test 
for  equality  of  several  population  medians.  Gibbons 
[Ref .  1  ;pp.  1 99]  . 

1 .    Computation  of  the  Test  Statistic 

Calculations   of  the  test  statistic  H   center 

around  the  ordered  arrangement  of  the  combined   samples 

from  which  the  sum  of  ranks  for  each  sample  is  derived. 

Let  Xij ,   j=1 ,2, . . . ,ni  and  i=1 ,2, . . . ,k,   be  independent 

random  samples  from  k  populations.  Let  r^-j  =  rank(Xi-j), 

Hi  k 

Ri=  23rii,andN=  2^i*   Then , 
j=1  i=1 

(  ^         ) 

H  =  (12/(N(N+1  ))   X;  (Ri2/ni)}-  3(N+1) 

(i=1  ) 
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If  ties  occur  in  the  combined  sample,  they 
are  resolved  by  assigning  the  average  value  of  the 
ranks  involved.  A  correction  based  on  the  number  of 
observations  tied  at  a  given  rank  and  the  number  of 
rsuiks  involved,  is  included  in  the  calculations.  A 
complete  description  of  the  correction  factor  is  given 
in  Gibbons  [Ref.  2:pp.  178-179]. 

2 .    The  Null  and  Asymptoxic  distribution  of  H 

The  null  distribution  of  H  is  generated  by 
enumeration.  Each  possible  permutation  of  ranks  is 
listed  for  the  combined  sample,  and  the  corresponding  H 
value  computed.  The  frequency  distribution  of  H.is  the 
total  number  of  occurrences  of  each  distinct  H  value. 
The  H  values  are  arranged  in  increasing  order  while 
maintaining  the  frequency  pairings.  The  null 
distribution  is  obtained  by  dividing  the  cumulative 
frequencies  by  n-|  !  n2  !  .  .  .  ni^  !  /N  !  . 

Due  to  computer  limitations,  generation  of 
the  exact  distribution  of  H  was  only  possible  for  k  =  3 
populations  with  n  =  4  observations  in  each,  and  4 
populations  with  3  observations  in  each.  Most  of  the 
distributions  were  generated  on  the  mainframe  computer 
and  saved  in  matrices  for  quick  recall  by  the  Xruskal- 
Wallis  test  program. 

Exact  C.D.F.  values  were  compared  with  the 
corresponding   approximate   values  using  the   following 
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distributions:  chi-square  with  (k-1)  degrees  of  freedom 
(C),  F  distribution  with  (k-1)  and  (N-k)  degrees  of 
freedom  (F),  and  F  with  (k-1)  and  (N-k-1 )  degrees  of 
freedom  (F1).  The  chi-square  distribution  uses  the 
Kruskal-Wallis  H  statistic,  while  the  F  and  F1 
distributions  use  a  modified  H  statistic,  HI  = 
( (N-k)H)/(k-1 )( (N-1 )-H) ;  see  Iman  and  Davenport  [Ref. 
8].  As  can  be  seen  in  Table  4,  F1  gives  the  mosx 
consistently  accurate  estimates. 


TA3LS   4..    C.D.F.   COMPARISONS   FOR  THE   XRUSKAL- 
WALLIS  TEST 

PRO  sen   >h::  FOR   .1    GROUP  OF   3  SAMPLES   CONSISTING  OF   ^ ,    u,  AND    3  Q3S. 

TEST   STAT.    VALUE    i  T.IU39  I  5.7121  T  5.1313  I  5.5985  I  5.J53d  1  a. 2121  I  3.5935 
EXACT   C.D.F.  I  .00970  i  .01905  !  .02961  I  .0U355  I  .07310  |  .12913  |  .1773U 

■ ^ + «  + + + + u 

ERROR:    CSISQUARS    I". 013^0  |~. 01532  I". 01535  I". 01220  T.OOlSu.  |  .007U6  |  .0124-1 

+- * + + -r -T + 

ERROR;    F  DIST           I  .0030^4  |  .00736  |  .00836  |  .01113  |  .01820  I  .01596  |  .00990 
+ + + + + + + 

ERROR;    F  W/'X   DF    I  .00084  j  .001+26  |  .00403  |  .00544  |  .01138  |  .00895  |  .00172 

PROBCH   2  h];  FOR  A    GROUP  OF   3  SAMPLES   CONSISTING  OF   ^ ,    4,  AND    4  OBS. 

TEST  STAT.    VALUE   I  7.6538  J  6.9615  I  6.5000  I  5.6923  I  4.9615  I  4.2692  I  3.5759 
+ + + + + + + 

EXACT   C.D.F.               I  .00762  I  .01939  |  .02996  |  .04366  |  .08000  |  .12190  |  .17299 
. + + + + + + + 

ERROR;    CBISQUARE    I". 01415  I". 01139  I". 00882  I". 00941  |~. 00368  |  .00361  |  .00577 
+ + + + + + + 

ERROR;    F  DIST           |  .00290  |  .00839  I  .01204  |  .01100  |  .01272  |  .01225  |  .00253 
+ + + + + + + 

ERROR;    F  »/~l   DF    I  .00149  |  .00600  I  .00890  |  .00647  |  .00708  |  .00592  I". 00384 

PROBCH   2  hi:  FOR  A   GROUP  OF   4  SAMPLES  CONSISTING  OF   3 .    3,  2,  AND   2   OBS. 

TEST   STAT.    VALUE   ]    7.6364  ]    7.1818  I  7.0000  I  6.5273  I  6.0182  T  5.3818  I  4.8727 
+ + + + + + + 

EXACT   C.D.F.               I  .01000  I  .01921  I  .02921  |  .04921  |  .07984  |  .12984  |  .17952 
+ + + + + + + 

ERROR;    CBISQUARE    I". 04415  I". 04712  j". 04259  |". 03939  |". 03089  I". 01604  |". 00183 

+ + H +  --- + -  + + 

ERROR;    F   DIST  !  .00294  |  .00250  I  .00730  1  .00330  I  .01094  |  .01123  I  .00902 

-  —  ----■.  —  --_---..--  — __------^--  —  ...  — _..^— _.-._  —  _-.^— _--.-  —  --.k  — ----  —  ---;-—-----  —  -  —  —  —  -  —  —  — — 

IRP.OR:    7  7/":  JF    ,".00172  !".0a4gu  !". 00147  |".J0307  i". 00375  l~. 00570  I". 00843 

PROBLB    2  h];  FOR  A    GROUP  OF   4  SAMPLES   CONSISTING  OF   3 ,    3,  3,  AND    2    OBS. 

TEST   STAT.    VALUE    I  8.0152  I  7.6354  ]    7.1515  I  5.7273  I  6.1970  I  5.4597  I  4.9697 
+ + + + + + + 

EXACT   C.D.F.               I  .00961  I  .01831  |  .02974  |  .04948  I  .07805  I  .12740  |  .17571 
+ + ^ + + + + 

ERROR;    CBISQUARE    j". 03609  j". 03534  |~. 03748  |". 03154  j". 02435  |". 01306  |  .00158 
+ + + + + + + 

ERROR;    F  DIST           I  .00215  1  .00481  |  .00441  |  .00920  I  .01185  I  .01035  |  .01214 
_ + + + + + + + 

ERROR;    FW/    IDF    I". 00133  |~. 00019  |". 00259  |  .00030  |  .00098  |". 00230  |". 00091 
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A  final  accuracy  comparison  "between  the  C  and  F1 
approximations  was  conducted  by  computer  simulation  for 
5  populations  with  8  observations  each.  Initially, 
30,000  permutations  of  the  4-0  ranks  were  randomly 
generated  (no  tie  ranks  allowed)  and  the  H  statistic 
calculated  for  each  permutation.  Then  the  empirically 
determined  percentiles  Hp  for  selected  values  of  p 
between  .01  and  .18  were  compared  with  the 
approximations  given  by  the  C  and  F1  distributions. 
The  results  are  shown  in  Table  5.  It  cam  be  seen  that 
the  F1  approximation  compares  well  with  the  simulated 
results,  giving  three  decimal  place  accuracy,  while 
the  C  approximation  is  less  accurate. 


TABLE   5.    C.D.F.   COMPARISONS  FOR   THE   KRUSKAL- 
WALLIS  TEST  USING  COMPUTER  SIMULATION 

PROBZE   2  h];  BASED   ON   100  00  GENERATED   B'S  FOR    5  SAMPLES   OF   8  OBS.    EACH. 

TEST  STAT.    VALUE   ]    12.229  I  11.055  I  10.2U8  I   9.212  I   8.129  T   7.030  J   5.232 
+ + + + + + + 

C.D.F.    VALUE              I  .01000  i  .02000  |  .03000  |  .05000  |  .08000  |  .13000  |  .18000 
+ + + + + + + 

ERROR:    C3ISQUARE    I". 00573  [".OOSSU  r.005U6  I". 00601  I". 00595  [".00432  r.002i+6 
+ + + + + + + 

ERROR;    F  W/'l   DF    \     .00081  |  .00231  I  .00259  j  .00350  |  .00151  |  .00109  I". 00103 

PR0BL3   2  h];  BASED   ON   20000  GENERATED   H'S   FOR    5  SAMPLES   OF   8  OBS.    EACH. 

TEST  STAT.    VALUE   ]    12.315  I  ll.OSu  ]    10.18U  J   9.153  I   3.129  I   7.034  J   5.220 
+ + + + + + + 

C.D.F.    VALUE              I  .01000  |  .02000  j  .03000  |  .05000  I  .08000  |  .13000  |  .18000 
+ + + + + + + 

ERROR;    CEISQUARE    |". 00516  I". 00596  r.007U5  I". 00716  |". 00696  [".00413  |~. 00334 
+ + + + + + + 

ERROR:    T   V/"!  DF    \     .00126  I  .00221  I  .00155  I  .00235  |  .00161  I  .00130  [",00199 

?R03Z:i    >  "- :  ;  3A.ZED    JN    3  0  0C0  ■JEr^ERATED    .7 '  S  70R    5  SAMPLES    OF    3  OBS.    EACH. 

TEST   STAT.    VALUE    ~\     12.305  [  10.975  I  10.147  [   9.179  I   8.158  j   7.072  I   6.265 
+ + + + + + + 

C.D.F.    VALUE              I  .01000  (  .02000  [  .03000  [  .05000  [  .08000  [  .13000  [  .18000 
+ + + + + + + 

ERROR;    CEISQUARE    [".00522  [".00684  [".00802  (".00577  [".00553  [".00213  [".00020 
+ + + + + + + 

ERROR;    F  V/~l   DF    \     .00121  [  .00143  I  .00111  |  .00273  [  .00301  |  .00344  [  .00142 
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3  .    Hypothesis  Testing; 

P-values  for  the  test  H0 :  the  population 
medians  are  all  equal  versus  HI :  at  least  two 
population  medians  are  not  equal,  are  computed  as: 
Pr  [H  2.  h]  ,  where  h  is  the  value  of  the  observed  test 
statistic . 

For  three  or  more  populations  with  at  least  4 
observations  in  each,   the  F1  approximation  is  used. 
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VI. TESTS  FOR  ASSOCIATION  IN  PAIRED-SAMPLES 

The  tests  described  herein  assume  that  the  data 
consists  of  independent  pairs  of  observations  (X±,Y±) 
from  a  bivariate  distribution.  The  general  null 
hypothesis  is  that  of  no  association  between  X  and  Y. 
Kendall's  3  and  Spearman's  R  are  considered. 

A .    KENDALL 'S3 

1  .    Computation  of  the  Tes":  Statistic 

The  test  statistic  is  computed  by  comparing 
each  observation  (Xj_,Yt_)  with  all  other  observations 
( Xj ,Yj )  in  the  sample.  If  the  changes  in  X  and  Y  are 
of  the  same  sign,  sgn(X-j  -  X^)  =  sgn(Yj  -  Y^  )  ,  the  pair 
(Xi,Yi)  and  (Xj,Y-j)  is  "concordant"  and  a  +1  is  scored. 
If  the  signs  are  different,  the  pair  is  "discordant" 
and  a  -1  is  scored.  Any  ties  between  either  the  X's  or 
the  Y's  scores  a  zero  for  that  pair.  The  sum  of  all 
scores  divided  by  the  total  number  of  distinguishable 
pairs,  (N(N-1))/2,  gives  B.  If  zeros  are  scored,  the 
denominator  is  reduced  by  a  correction  factor  which  is 
based  on  the  number  of  observations  tied  ax  a  given 
rank  and  the  number  of  ranks  involved  in  each  of  the  X 
and  Y  samples.  A  complete  description  of  the 
correction   for   ties  is  given  in  Gibbons   [Ref.   2:pp. 
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289] .  The  value  of  B  ranges  between  1 ,  indicating 
perfect  concordance,  and  -1,  for  perfect  discordance. 
Gibbons  [Ref.  1:pp.  209-225]. 

2 .    The  Null  and  Asymptotic  Distribution  of  B 

.  The   null  distribution  of  B  is  derived   from 
the   following  recursive  formula  given  in  Gibbons  [Ref. 
1  :pp  .216]  . 

u(N+1,P)  =  u(N,P)  +  u(N,P-1)  +  u(N,P-2)  +  ...+  u(N,P-N) 

where  u(N,P)  denotes  the  number  of  P  concordant 
pairings  of  N  ranks.  This  formula  is  used  to  generate 
the  frequency  with  which  the  possible  values  of  P 
occur.  Division  by  N!  results  in  the  probability 
distribution  of  P.  Since,  B  =  ( 4-P/ (N(N-1  )  )  )-1  ,  the 
null  distribution  of  B  is  easily  determined. 

Exact  C.D.F.  values  were  compared  with  those 
obtained  using  a  normal  approximation,  with  and  without 
a  continuity  correction  factor  (CO  =  6/N(n2-i),  pro- 
posed by  Pittman  [Ref.  11]  for  the  Spearman's  R  test). 
The  results  for  various  sample  sizes  are  provided  in 
Table  6.  As  can  be  seen,  for  sample  sizes  greater  than 
13,  a  normal  approximation  with  continuity  correction 
provides  three  decimal  place  accuracy. 
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TABLE    6.        C.D.F.     COMPARISONS    FOR    KENDALL'S    3 
PP03C3  2   b];    JQ2  3AMPLI  3122  ZQUAL  TO  13. 


TI3T  STAT.  '-^ALUS 

-T 

!  0.5123 

1  0.4bi5 

1                { 

1  0.4103  1 

1 
0.3590  1 

1 

0.3233  i 

:  0.2564  ! 

i  0.2303 

IXACT  :.2.i. 

1  .00743 

I 

1  .01524 

1  .02363  1 

.04999  ! 

1 

.06443  1 

•   1 

.12593  1 

.15309 

Z220S;  NORMAL 

1  .00014 

1  .00121 

1  .00313  1 

.00620  1 

1 

.00303  ! 

1 

.01473  1 

.01703 

Z220?.;  MOP.M.  WCC 

;  ".00013 

1  .00073 

;  .00239  ! 

.00497  ; 

.00653  1 

,01223  i 
i 

,01415 

PP0BC3  > 

b];  FOR  3 

AMPLE  31 ZI 

IQUAL  TO 

14. 

TEST  STAT.  VAL'JS 

r  

!  0.4725 

i 

!  0.4236 

!  0.4066  ! 

0.3626  ! 

0.2967  ! 

0.2527  1 

0.2033 

IXACT  :.J.7, 

.  w  \^  -  0  -^ 

:    .Ol'T-'S 

!  _  02359  ' 

.03973  1 

.07353  i 

.11656  ' 

.165+1 

I/IHOK;  NORMAL 

1  .00035 

1  .00140 

!  .00213  1 

.00431  1 

.00339  1 

.01256  1 

.01627 

iRSO?.:  NORM.  'J/'CC 

1  .00003 

i  .00095 

1  .00152  ! 

.00345  1 

,00742  i 

.01057  1 

.01371 

3.    Hypothesis  Testing 

P-values  for  tests  of  no  association  between 
X  and  Y  are  computed  for  three  types  of  alternative 
hypotheses.  Because  the  distribution  of  B  is  symmetric, 
all  probabilities  can  be  taken  from  the  upper  tail 
using  the  absolute  value  of  b,  the  observed  value  of 
the  test  statistic.  Linear  interpolation  is  used  when 
b  lies  between  tabulated  values.  The  P-values  are 
computed  as  follows. 

a.  One-Sided  Alternatives 

The   one-sided  alternative  tested   depends 
on  the  sign  of  b.   A  positive  b  will  automatically  test 
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for  direct  association  or  concordance,  while  a  negative 
b  will  test  for  indirect  association  or  discordance. 
The  P-value  equals  Pr[B  >_    |b|]. 

b.  Two-Sided  Alternative 

The   P-value  equals  twice  the   probability 
computed  for  the  one-sided  hypothesis. 

For  sample  sizes  greater  than  12,  a  normal 
approximation  with  continuity  correction  is  used.  The 
approximate  P-values  are  then: 

1  -  Pr[  Z  ±    ((|b|  -  CO  -  M^)/    <7>3]  ,   where   Z 
is  standard  normal,   CC   is  the   continuity  correction, 
/ih)  =  0,    and  <^^~    =    (4N  +  10)/9N(N-1),   for   the   one- 
sided  test   and  twice  this  P-value  for   the   two-sided 
test  . 

B.    SPEARMAN'S  R 

The  Spearman's  R  Test  requires  the  added 
assumption  that  the  underlying  bivariate  distribution 
is  continuous.  The  test  measures  the  degree  of 
correspondence  between  rankings,  instead  of  the  actual 
variate  values,  and  can  be  used  as  a  measure  of 
association  between  X  and  Y.  Gibbons  [Ref.  1:pp.  226]. 

1 .    Computation  of  the  Test  Statistic 

The  test  statistic  R  is  computed  in  the 
following  manner.  Let  ri  =  rank(Xi)  and  Si  =  rank(Yi) 
and   Di  =  r^  -  Si .   Then, 
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N 
6EDi2 

R   =  1  -  i=j 

N(n2  -  1 ) 


where  N  is  the  size  of  the  sajnpLe.  If  ties  occur  in  X 
or  Y,  they  are  resolved  by  assigning  the  average  value 
of  the  ranks  involved.  A  correction  factor,  based  on 
the  number  of  observations  tied  at  a  given  rank  and  the 
number  of  ranks  involved,  is  included  in  the 
calculations.  A  complexe  description  of  the 
correction  factor  is  given  in  Gibbons  [Ref.  2:pp.  279]. 
The  value  of  R  ranges  between  1 ,  indicating  perfect 
direct  association,  and  -1 ,  for  perfect  indirect 
association.  Gibbons  [Ref.  1:pp.  226-235]. 

2 .    The  Null  and  Asymptotic  Distribution  of  R 

The  null  distribution  of  R  for  a  given  sample 
size  N  is  generated  by  enumeration.  The  method,  as 
presented  in  Kendall  [Ref.  9],  involves  generation  of 
an  N  by  N  array  of  all  possible  squared  differences 
between  any  two  paired  ranks  of  X  and  Y.  All  N! 
permutations  of  N  ranks  are  used  to  index  values  from 
the  array.  The  sum  of  these  indexed  values  for  each 
permuxation  gives  rise  to  M!  sum  of  squared  differences 
which  are  then  converted  to  the  R  statistic.  The 
frequency  distribution  of  R  is  the  total  number  of 
occurrences  of  each  distinct  value  of  R  divided  by  N! . 
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Due  to  mainfrajne  computer  memory  limitations 
in  the  APL  environment,  generation  of  the  distribution 
of  R  was  limited  to  sample  sizes  of  7  or  less.^  Using 
tables,  provided  by  Gibbons  [Ref.  2:pp.  417-418]  to 
supplement  computer  computations,  a  numerical  matrix, 
called  PMATSP,  was  created  ,to  store  the  cumulative 
distributions  of  R  for  sample  sizes  less  than  1 1 .  This 
matrix  allows  for  quick  recall  of  cumulative 
probabilities  by  the  Spearman's  R  Test  program. 

Exact  C.D.F.  values  were  compared  with  those 
obtained  using  a  student's  T  approximation  with  (N-2) 
degrees  of  freedom  (see  Glasser  and  Winter  [Ref-  10]), 
and  a  normal  approximation.  Both  normal  and  T 
approximations  were  computed  with  and  without  a 
continutity  correction  factor,  CC  =  6/N(n2-i)  (Pittman 
[Ref.  11]).  From  the  results  presented  in  Table  7,  the 
most  consistently  accurate  approximation  is  given  by 
the  T  distribution  with  a  correction. 
3 .    Hypothesis  Testing 

P-values  for  tests  of  no  association  between 
X  and  Y  can  be  computed  for  three  types  of  alternative 
hypotheses.  Because  the  distribution  of  R  is  symmetric, 
ail   probabilities  are  taken  from  the  upper  tail   using 


1  The   memory  capacity  of  the  mainframe   computer   in 
the  APL  environment  is  limited  to  2.5  megabytes. 
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TABLE  7.   C.D.F.  COMPARISONS  FOR   SPEARMAN'S  R 
P20BC2  >  r];  705?  3AMPLI  SIZ2  IQUAL  TO  9. 


TZ3T  STAT.  i-'ALUS 

1  0.7333 

1 

I  0.7167 

-: r 

1  0.6667  1 

0.6000  1 

0.5333 

1  0.4333 

!  0.3500 

iXACT  C.3.?. 

1  .ooasi 

1  .01343 

i  .02944  1 

.04340  1 

.07375 

!  .12496 

1  .17929 

12202;  MORflAL 

I  '.00475 

1  ".00230 

1  -.00023  1 

.00356  ! 

.00305 

!  .01479 

1  .01319 

12202;  N02M.  WOZ 

1  ■.00553 

1  ".00413 

1  ".00135  ! 

.00123  1 

.00493 

1  .01029 

!  .01236 

12202;  T  0137 

!  .00235 

i  .00352 

1  .00451  i 

.00459  1 

.00415 

1  .00293 

!  .00133 

Z22QH;  T  U/CC 

i  .00153 

1  .00203 

1 

1  .00251  1 

.00175  1 

.00042 

I  -.00212 

!  ".00479 

?2C3C?.  1 

rJ:  702  3 

AMPIZ  3IZ2 

IQUAL  73 

10. 

TT3T  STAT.  '^ALUZ 

1  0.7455 

1  0.6727 

1  0.5364  i 

0.5636  1 

0.4909  i 

1  0.4061 

i  0.3333 

;  .00370 

1  .01343 

!  .02722  ! 

.04314  1 

.07741  ! 

!  .12374 

1  .17437 

12202;  .M02f-.fiL 

.  ".  :0393 

i  ".00220 

i  -.00091  1 

1               » 

.00271  1 

.00700  1 

I 

.01216  1 

i  .01572 

12202;  NORM.  U/CC 

i -.00457 

I 

1  ".00327 

1 

1  -.00210  1 
1       1 

.00095  1 

.00451  1 
1 

.00867  1 
1 

1  .01123 

12202;  T  DIST 

1 

1  .00204 

1 

1 

1  .00236 

1 

1       1 
1  .00326  1 

1       1 

.00323  1 

1 
.00253  1 

1 
.00159  1 

I 

.00107 

12202;  T  U/CC 

1  .00144 

1 

1  .00135 

1       i 
1  .00134  1 
1       1 

.00114  1" 

i 
'.00035  1 

1 

".00230  1 

1 

-.00361 

the  absolute  value  of  r,  the  observed  value  of  the  test 
statistic.  The  P-values  are  computed  as  follows, 
a.  One-Sided  Alternatives 

The   one-sided  alternative  tested   depends 

on  the  sign  of  r .  A  positive  r  vill  test  for  direct 
association,  wnile  negative  r  tests  for  indirect 
association.   The  P-value  equals  Pr  [R  2.  1 1"  I  ]  • 
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b.  Two-Sided  Alternative 

The   P-value  equals  twice  the   probability 
computed  for  the  one-sided  hypothesis. 

For  sample  sizes  greater  than  10,  an 
approximation  based  on  the  student's  T  distribution 
with  (N-2)  degrees  of  freedom  and  continuity 
correction,  is  used-   The  P-value'S  are: 

1  -  Pr[  T(N-2)  1  ((kl  -  CO  -  ^^)/  (7p], 
where  T(n-2'>  denotes  the  T  distnibution  with  (N-2) 
degrees  of  freedom,  CC  is  the  continuity  correction, 
A/p  =  0,  and  Or--  =  (1  -(  I  r  ! -CC  )2  )  /  ( N-2  )  ,  for  the  one- 
sided test,  and  twice  this  P-value  for  the  two-sided 
test.   Gibbons  [Ref.  1:pp.  218]. 
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VII.NONPARAMETRIC  SIMPLE  LINEAR  REGRESSION^ 

Nonparametric  Linear  Regression  assumes  that  the 
data  consists  of  independent  pairs  of  observations  from 
a  bivariate  distribution  and  that  the  regression  of  Y 
on  X  is  linear.  The  program  estimates  linear 
regression  parameters  based  on  the  data  samples.  It 
then  allows  the  user  to  input  X  values  to  predict  the  Y 
values.  Hypo~hesi3  testing  and  confidence  interval 
estimation  for  the  slope  of  the  regression  equation  is 
offered.  If  the  estimated  slope  lies  outside  the 
confidence  interval ,  an  alternate  regression  equation 
is  offered  with  an  opportunity  to  input  X  values  to 
predict  the  corresponding  Y  values. 

A.    COMPUTATION  OF  THE  ESTIMATED  REGRESSION  EQUATION 

The  least  squares  method  is  used  to  estimate  A  and 
B  in  the  regression  equation  Yi  =  A  +  BXi  +  e± 
( i= 1 , 2 , . . .N ) ,  where  ei  ( unobservable  errors)  are 
assumed  to  be  independent  and  identically  distributed. 
A  and  B  are  computed  from  the  following  equations: 


^Except  for  program  design  considerations,  the 
information  and  concepts  provided  in  the  section  are 
paraphrased  from  Conover  [Ref.  12:pp.  263-271]. 
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N  N  N 

NE^iYi       -        E  Xi  E  Yi 

B    =    i^J i=  1        i=  1 

N  /   N        \2 

nE:<i2     -        Exi 

i  =  1  ^i  =  1      ' 


N  N 

EYi       -       BEXi 

A    =  i  =  l i  =  1 

N 


B.  HYPOTHESIS  TESTING 

P-values  for  testing  hypotheses  about  the  slope  of* 
the  regression  equation  are  based  on  the  Spearman's 
rank  correlation  coefficient  R  between  the  X^  and  Uj_  = 
Yt_  -  B0Xj_  ,  where  3^  is  "che  hypothesized  slope.  The 
appropriate  one-sided  test  of  hypothesis,  H0 :  B  =  Bg 
versus  HI :  B  <  Bq  or  HI :  B  >  Bq ,  is  automatically 
chosen  based  on  the  sign  of  the  computed  test  statistic 
r  (positive  r  tests,  HI :  B  >  B0 ;  negative  r  tests,  HI : 
B  <  B0 )  .  The  P-value  is  computed  as:  Pr  [R  2.  1 1"  I  ]  • 
P-values  for  two-sided  tests,  H0 :  B  =  B0  versus  HI  :  B  ^ 
B0 ,  are  also  presented. 

For  sample  sizes  N  greater  than  10,  P-values  are 
approximated  using  a  T  distribution  with  (N-2)  degrees 
of  freedom  and  continuity  correction. 

C.  CONFIDENCE  INTERVAL  ESTIMATION 

100(1-  a  )%  confidence  bounds  for  the  slope 
parameter  B  are  determined  as  follows.   The  n   possible 
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slopes,  Si_i  =  ( Yj_-Y-j  )/ (  Xj_-X-i  )  ,  are  computed  for  all 
pairs  of  data  (Xi,Yj_)  and  (X-j,Yj)  such  that  i  <  3  and 
Xi  #  Xj  and  rearranged  in  increasing  order  to  give  s(  "•  ) 
<_  s(2)  _<_...  _<_  s(^).  Let  w  be  the  (1-  a/2)  percentile 
of  the  distribution  of  Kendall's  statistic  with  sample 
size  n . ''  Let  d  be  the  largest  integer  less  than  or 
equal  to  (n-w)/2  and  u  the  smallest  integer  greater 
than  or  equal  to  (n+w)/2  +  1.  Then  s(d-)  and  s(^)  are 
the  desired  lower  and  upper  confidence  bounds, 
respectively . 

For  sample  sizes  larger  than  13,  a  normal 
approximation  with  continui.ty  correction  is  used  to 
estimate  the  confidence  intervals. 

If  the  slope  of  the  estimated  regression  equation 
does  not  lie  within  the  computed  confidence  interval, 
the  prograjn  automatically  calculates  a  new  regression 
equation  where  the  slope  is  the  median  of  the  two-point 
slopes  Si-j  and  the  intercept  is  the  difference  of  the 
medians  of  the  X  and  Y  samples,  My-  Mx . ^ 


"•Kendall's  statistic  is  defined  here  as  Nq  -  N(j[, 
where  Nq  is  the  number  of  concordant  pairs  of 
observations  and  N4  is  the  number  of  discordant  pairs. 
Conover  [REF.  12:pp.  256]. 

^This  procedure  is  recommended  by  Conover  [REF. 
12:pp.  256]. 
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VTII .AREAS  FOR  FURTHER  WORK 

To  create  a  more  versatile  and  powerful  software 
package,  the  NONPAR  workspace  could  be  expanded  to 
include  some  or  all  of  the  following  nonpar ametric 
xests:  tests  for  randomness  based  on  runs,  Chisquare 
and  Kolmogorov-Smirnov( K-S )  Goodness-of -f i t  tests, 
Chisquare  and  K-S  general  two  sample  distribution 
nesTs,  Chisquare  test  for  independence,  and  the 
Friedman  test  for  association. 
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APPENDIX  A 
DOCUMENTATION   FOR  THE  MICROCOMPUTER  WORKSPACE 

1 .  General  Information 

This  appendix  describes  the  organization  and 
operation  of  the  I3M-PC  (or  compatible)  version  of  the 
workspace.  Appendix  C  continues  from  where  this 
appendix  leaves  off,  to  walk  the  user  through  each  test 
by  working  practical  examples. 

Before  proceeding  any  further,  the  user  should 
refer  to  section  II  (Workspace  Design  Issues)  for 
general  information  about  workspace  requirements  and 
assumptions  regarding  its  use. 

To  get  started,  enter  the  APL  environment  in  the 
usual  manner  and  load  the  NONPAR  workspace. 

2 .  Workspace  Menus 

This  workspace  is  designed  around  the  use  of 
menus.  They  guide  the  user  through  the  selection 
process  of  choosing  a  nonparametr ic  test  and  a  test 
option.  Three  types  of  menus  are  used;  the  main  menu, 
test   menus,   and  help  menus. 

a.    The  Main  Menu 

Within  moments  of  loading  the  NONPAR 
workspace,  the  main  menu  will  appear.  It  is  titled 
Nonparametr ic   Statistical  Tests.    This  menu   presents 
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general  information  abcut  the  workspace.  Its  primary 
purpose  is  to  list  the  choices  of  nonpar ametric  tests 
available  and  provide  an  option  which  allows  the  user 
to  exit  the  main  menu  into  APL  to  copy  data  into  the 
workspace  or  return  to  DOS.  Each  test  choice  is  listed 
with  some  information  about  the  test's  area  of 
application.  To  make  a  selection  from  the  menu,  move 
the  cursor  (using  the  cursor  keys)  to  highlight  the 
desired  choice,  and  press  enter.  As  a  reminder  to  the 
user.  a  footnote  at  the  bottom  of  the  screen  describes 
the  procedure  for  entering  a  choice.  Once  a  test  has 
been  selected  from  the  main  menu,  a  sub-menu 
appropriate  to  the  test  appears.  To  exit  from  any  menu 
back  to  the  main  menu,  press  the  Escape  key. 
b .    Test  Menus 

The  title  of  the  test  menu  is  the  name  of  the 
nonparajnetr ic  test  chosen.  The  text  portion  of  the 
menu  gives  a  general  overview  of  the  test,  to  include, 
the  method  used  to  compute  the  test  statistic,  and  a 
description  of  the  various  options  that  may  be 
exercised.  The  third  section  consists  of  the  list  of 
test  options  available.  These  options  include 
returning  to  the  main  menu  or  choosing  the  help  menu. 
Test  menus  may  have  options  listed  in  single  or 
multiple-paged  formats.  The  comment  in  the  final  block 
of   the   menu  lets  the  user  know  if  a  certain   menu   is 
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multiple-paged  or  not.  To  make  a  selection  from  a 
multiple-paged  qjenu,  use  the  page-up  or  page-down  key 
to  locate  the  desired  option.  Proceed  with  the  scroll 
keys  to  highlight  the  choice.  and  press  enter.  Once  a 
test  option  is  entered,  the  user  is  prompted  to  input 
the  data  required  to  run  the  test.  When  the  option  for 
more  information  is  selected,  the  help  menu  is 
displayed. 

c.    Help  Menus 

The  -i:itle  of  the  help  menu  usually  begins 
with  the  words  "More  Information  About..."  followed  by 
the  title  of  the  nonpar  ainetr  ic  test.  The  text  portion 
of  the  menu  explains  the  test  and  its  options  in 
greater  detail.  No  choices  are  offered  in  the  menu. 
To  return  to  the  test  menu,  press  any  key. 
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APPENDIX  3 
DOCUMENTATION   FOR   THE  MAINFRAME  COMPUTER  WORKSPACE 

1 .  General  Information 

This  appendix  describes  the  organization  and 
operation  of  the  mainframe  computer  workspace.  To  load 
a  copy  of  the  NONPAR  workspace  from  the  APL  library, 
enter  the  APL  environment  and  type:  )LOAD  9  NONPAR. 
Within  a  few  moments  the  variables  LIST  and  DESCRIBE 
are  displayed  on  the  screen.  These  variables  provide  a 
description  of  the  workspace. 

2.  The  NONPAR  Works-pace 

The  NONPAR  workspace  consists  of  seven 
programs  which  call  several  subprograms  during  their 
execution.  The  exact  syntax  for  each  test  and  its 
corresponding  nonparametr ic  test  name  is  given  in  the 
following  format: 

SYNTAX:  Nonpar ametric  Test  and  Application. 

a.  SIGN:   Ordinary   Sign  Test   for  Location  in  Single 
and  Paired-sample  Data. 

b.  WILCOX:   Wilcoxon   Signed-rank  Test   for   Location 
in  Single  and  Paired-sample  Data. 

c.  MANNWHIT:   Mann-Whitney  Test   for  Equal  Medians  or 
Variances  in  Two  Independent  Samples. 

d.  KRUSKAL:   Kruskal-Wall is  Test  for  Equal  Medians  in 
K  Independent  Samples. 
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e.  KENDALL:   Kendall's  B;   Measure  of  Association  for 
Paired-sample  Data. 

f.  SPEARMAN:   Spearman's   R;   Measure   of  Association 
Between  Rankings  of  Paired  Data. 

g.  NPSLR:    Nonpar ametric   Simple  Linear   Regression; 
Least  Squares. 

The  list  presented  above  can  be  displayed  at  any 
time  by  typing:  LIST. 

Eor  each  test  program,  there  exists  a  HOW  variable 
that  gives  a  full  description  of  the  test  and  the 
various  options  "chat  may  be  exercised-  To  display  any 
of  ohe  HOW  variables,  Just  enter  trie  test  program's 
name  with  the  suffix  HOW  appended   (i.e.   3IGNH0W). 

A  test  is  run  by  entering  the  program's  name.  The 
user  is  immediately  prompted  to  inpux  data.  Enter 
numerical  data  separated  by  spaces  or  as  a  variable  to 
which  the  numbers  have  been  previously  assigned. 
Several  of  the  tests  require  a  considerable  amount  of 
prompting  before  all  the  necessary  data  has  been 
entered. 


53 


APPENDIX  C 
WORKSPACE  FAMILIARIZATION  THROUGH  PRACTICAL  EXAMPLES 

1 .  General  Information 

This  appendix  applies  to  both  the  mainfrajtne  and 
microcomputer  workspaces.  Its  purpose  is  to  acquaint 
the  user  with  the  organization  of  the  programs  and  the 
type  of  prompts  to  be  expected- 

Extensive  error  checking  has  been  included  in  the 
programs  to  ensure  that  the  data  is  of  the  proper  form. 
Should  a  prograjn  become  suspended,  clear  the  state 
indicator  by  entering:  )RESET,  check  over  the  data  for 
errors,  and  restart  the  program.  350  kilobytes  of 
computer  memory  are  needed  to  load  APL  and  the  NONPAR 
workspace;  to  avoid  filling  up  the  remaining  workspace 
area,  the  user  should  minimize  data  storage  in  the 
NONPAR  workspace.  To  exit  a  program  at  any  time,  press 
the  Control  and  Escape  keys,  simultaneously. 

2 .  Practical  Examples 
a.  Sign  Test 

( 1 )  Description  of  Problem  1 .  A  Sinclair  mine 
is  manufactured  to  have  a  median  explosive  weight  of 
not  less  than  16  ounces.  The  explosive  weights  of  15 
mines,  randomly  selected  from  the  production  line,  were 
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recorded  as  follows:   16.2  15.7  15.9  15-3  15.9  16   16.1 
15.3  15.9  16  16.1  15.7  15.8  15.9  15.8. 

(a)  Is  the  manufacturing  process  packing 
enough  explosives  in  the  mines? 

(b)  What  range  of  values  can  be  expected 
for  the  median  of  the  explosive  weights  90^  of  the 
time  . 

(2)  Solution .  To  see  if  the  manufacturing 
process  is  meeting  the  specifications,  we  test  the 
hypothesis  Ho:  M  =  16  versus  HI  ;  M  <  16. 

(3)  Workspace  Decision  Process. 

(a)  Microcomputer:  Choose  the  Sign  Test 
from  the  main  menu,  and  the  option.  Single  Sample; 
Test  H0 :  M  =  Mo  versus  HI :  M  <  Mo,  from  the  test  menu. 
Skip  to  the  Program  Interaction  section  below. 

(b)  Mainframe:  Enter  SIGN  at  the  keyboard 
and  receive  the  prompt: 

DID  YOU  ENTER  THIS  PROGRAM  FOR  THE  SOLE 
PURPOSE  OF  GENERATING  CONFIDENCE  INTERVALS  FOR  A 
SPECIFIED  SAMPLE  SIZE  AND  QUANTILE?  (Y/N). 

Enter  N  ( If  Y  is  entered,  the  user  will 
2,0  directly  to  this  last  option  of  the  test).  The  next 
prompt  is: 

THE  NULL  HYPOTHESIS  STATES  -  THE 
POPULATION  MEDIAN  (M)  IS  EQUAL  TO  THE  HYPOTHESIZED 
MEDIAN  (Mo);  H0 :  M  =  Mo.   WHICH  ALTERNATIVE  DO  YOU  WISH 
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TO  TEST?    ENTER:   1  FOR  HI :  M  <  Mo;   2  FOR  HI :  M  >  Mo; 
3  FOR  HI  :  M  i^    Mo. 

Enter  1 .   The  next  prompt  is: 
ENTER:     1  FOR  SINGLE-SAMPLE  PROBLEM;  2 
FOR  PAIRED-SAMPLE  PROBLEM. 

Enter  1 . 
(4)  Program  Interaction.   The  prompt  is: 

ENTER  THE  DATA  (MORE  THAN  TWO  OBSERVATIONS 
ARE  REQUIRED) . 

Enter  the  data  separated  by  spaces  or  as  a 
variable  to  which  the  data  has  been  previously 
assigned.   The  next  prompt  is: 

ENTER  THE  HYPOTHESIZED  MEDIAN. 

Enter  16.   The  following  is  dislayed. 

COMPUTATIONS  ARE  BASED  ON  A  SAMPLE  SIZE  OF: 
13. 

THE  TOTAL  NUMBER  OF  POSTIVE  SIGNS  IS:  3. 

THE  P-VALUE  FOR  H0 :  M  =  16  Versus  HI: 
M  <  1 6  IS:   .0461  . 

Consider  a  significance  level  of  .05. 
Since  the  P-value  of  .0461  is  less  than  .05,  we  reject 
H0 :  M  =  16  in  favor  of  HI :  M  <  16  and  conclude  that  the 
manufacturing  process  is  not  packing  enough  explosives 
in  the  Sinclair  mine.   The  next  prompt  is: 

WOULD  YOU  LIKE  A  CONFIDENCE  INTERVAL  FOR 
THE  MEDIAN?  (Y/N) . 
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Enter  Y  ( If  N  is  entered,  the  progam  asks 
if  confidence  intervals  for  a  quantile  are  desired). 
The  next  prompt  is: 

ENTER   THE  DESIRED  CONFIDENCE  COEFFICIENT; 
FOR  EXAMPLE:  ENTER  95,  FOR  A  95^  CONFIDENCE  INTERVAL. 

Enter  90.   The  following  is  displayed. 

A  90%  CONFIDENCE  INTERVAL  FOR  THE  MEDIAN  OF 
THE  POPULATION  IS:   (  1  5  . 3  j<_  MEDIAN  _<_  16  ). 

The  next  prompt  is: 

WOULD  YOU  LIKE  CONFIDENCE  INTERVALS  FOR  A 
SPECIFIED  QUANTILE?  (Y/N). 

To  see  the  form  of  the  results,  we  generate 
confidence  intervals  for  the  30th  quantile.  Sample 
size  is  automaticly  set  at  the  number  of  data  points 
entered  earlier.  Enter  Y  (If  N  is  entered,  the 
mainfrajne  program  ends;  or,  the  Sign  test  menu 
reappears).   The  next  prompt  is: 

ENTER  DESIRED  QUANTILE;  FOR  EXAMPLE:  ENTER 
20,  FOR  THE  20TH  QUANTILE. 

Enter  30.   The  following  is  displayed. 


ORDER  STATISTICS 
3       8 
2       9 
1       10 


COEFFICIENTS 
.8231 60 

.949490 
.  991 600 


*****  THIS  TABLE  GIVES  CONFIDENCE  COEF- 
FICIENTS FOR  VARIOUS  INTERVALS  WITH  ORDER  STATISTICS  AS 
END  POINTS  FOR  THE  30TH  QUANTILE. 
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The  mainfrajne  program  ends.  The  menu- 
driven  microcomputer  program  pauses  for  input  from  the 
keyboard  by  prompting: 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  the  Sign  test  menu 
reappears . 

b.  Wilcoxon  Signed-rank  Test 

(1)  Description  of  Problem  2.  A  special  train- 
ing program  is  being  considered  to  replace  the  regular 
training  that  Radio  Telephone  Operators  receive.  In 
order  to  evaluate  the  effectiveness  of  the  new  training 
program,  proficiency  tests  were  given  during  the  third 
week  of  regular  training.  Twenty-four  trainees  were 
chosen  at  random  and  grouped  into  twelve  pairs  based  on 
proficiency  test  scores.  One  member  of  each  pair 
received  specialized  training  while  the  other  member 
received  regular  training.  Upon  graduation,  the 
proficiency  tests  were  given  again  with  the  following 
results . 

Specially  Trained  Group  (X):  60  50  55  71 
43  59  64  49  61  54  47  70 

Regularly  Trained  Group  (Y):  40  46  60  53 
49  57  51  53  ^5  59  40  35 

(a)  Does  the  special  training  program 
ensure  higher  scores? 
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(b)  By  what  range  of  values  can  the 
scores  of  the  two  groups  be  expected  to  differ  95%  of 
the  time? 

(2)  Solution .  To  test  the  hypothesis  that  the 
special  training  program  raises  profficiency  scores,  we 
test  H0:  M(X-Y)  =  0  versus  HI:  M(X-Y)  >  0. 

( 3 )  Workspace  Decision  Process 

(a)  Microcomputer:  Choose  the  Wilcoxon 
Signed-rank  Test  from  the  main  menu,  and  the  option. 
Paired-sample;  Test  H0 :  M  =  Mo  versus  HI :  M  >  Mo,  from 
the  test  menu.  Skip  to  the  Program  Interaction  section 
below. 

(b)  Mainframe:  Enter  WILCOX  at  the  key- 
board and  receive  the  prompts: 

THE  NULL  HYPOTHESIS  STATES  -  THE 
POPULATION  MEDIAN  (M)  IS  EQUAL  TO  THE  HYPOTHESIZED 
MEDIAN  (Mo);  H0 :  M  =  Mo.  WHICH  ALTERNATIVE  DO  YOU  WISH 
TO  TEST?  ENTER:  1  FOR  HI :  M  <  Mo;  2  FOR  HI :  M  >  Mo; 
3  FOR  HI  :  M  jt  Mo. 

Enter  2.   The  next  prompt  is: 

ENTER:  1  FOR  SINGLE-SAMPLE  PROBLEM;  2 
FOR  PAIRED-SAMPLE  PROBLEM. 

Enter  2. 

(4)  Program  Interaction.   The  prompt  is: 
ENTER   X  DATA  (MORE  THAN   TWO   OBSERVATIONS 

ARE  REQUIRED) . 
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Enter  the  X  data  separated  by  spaces.    The 
next  prompt  is: 

ENTER  Y  DATA  (NUMBER  OF  Y  ENTRIES  MUST 
EQUAL  NUMBER  OF  X  ENTRIES). 

Enter  the  Y  data.   The  next  prompt  is: 

ENTER  THE  HYPOTHESIZED  MEDIAN  FOR  THE 
DIFFERENCES  OF  THE  PAIRED  DATA. 

Enter  0.   The  following  is  dislayed. 

COMPUTATIONS  ARE  BASED  ON  A  SAMPLE  SIZE  OF: 
12. 

THE  TOTAL  SUM  OF  POSITIVE  RANKS  IS:   60.5. 

THE  P-VALUE  FOR  COMPARING  THE  MEDIAN  OF  THE 
POPULATION  OF  DIFFERENCES  TO  THE  HYPOTHESIZED  MEDIAN, 
H0:  M(X-Y)  =  0   Versus  HI:  M(X-Y)  >  0,  IS:   .0505. 

Consider  a  significance  level  of  .05. 
Since  the  P-value  of  .0505  is  greater  than  .05,  we  do 
not  reject  the  null  hypothesis  that  the  two  training 
cources  are  equally  effective.  However,  due  to  the 
closeness  in  values,  the  choice  of  rejecting  or  not 
rejecting  the  null  hypothesis  is  strictly  a  judgement 
call.   The  next  prompt  is: 

WOULD  YOU  LIKE  A  CONFIDENCE  INTERVAL  FOR 
THE  MEDIAN?  (Y/N) . 

Enter  Y  ( If  N  is  entered,  the  mainframe 
progaun  ends;  or,  the  Wilcoxon  test  menu  reappears). 
The  next  prompt  is: 
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ENTER  THE  DESIRED  CONFIDENCE  COEFFICIENT; 
FOR  EXAMPLE:  ENTER  95,  FOR  A  95%  CONFIDENCE  INTERVAL. 

Enter  95.   The  following  is  displayed. 

A  95%  CONFIDENCE  INTERVAL  FOR  THE  MEDIAN  OF 
THE  POPULATION  OF  DIFFERENCES  IS: 

(  -1  <_   MEDIAN(X-Y)  <.  16.5  ). 

The  mainframe  progrsun  ends.  The  menu- 
driven  microcomputer  program  pauses  for  input  from  the 
keyboard  by  prompting: 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  the  Wilcoxon  test  menu 
reappears . 

c.  Mann-Whitney  Test  for  Equality  of  Medians 

( 1  )  Description  of  Problem  3.  A  group  of  Army 
and  Navy  officers  were  given  the  Defense  Language 
Aptitude  test.  From  the  results,  14  Army  and  17  Navy 
officers'  scores  were  randomly  selected.  These  scores 
are  listed  below. 

Army  (X):  35  30  55  51  28  25  1 6  63  60  44  20 
42  47  38. 

Navy  (Y):  54  26  41  43  37  34  39  50  46  49  45 
33  29  36  38  42  34. 

(a)  Is  there  sufficient  evidence  to  claim 
that  Navy  officers  score  higher  on  this  test  than  Army 
officers? 
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(b)  By  what  range  of  values  can  the 
scores  between  the  two  groups  be  expected  to  differ  90^ 
of  the  time. 

(2)  Solution.  To  see  if  Navy  officers  score 
higher  on  the  exajn,  we  test  H0 :  Mx  =  My  versus 
HI :  Mx  <  My. 

( 3 )  Workspace  Decision  Process. 

(a)  Microcomputer:  Choose  the  Mann- 
Whitney  Test  from  the  main  menu,  and  the  option,  Test 
H0 :  Mx  =  My  versus  HI :  Mx  <  My,  from  the  test  menu. 
Skip  to  the  Program  Interaction  section  below. 

(b)  Mainframe:  Enter  MANNWHIT  at  the 
keyboard  and  receive  the  prompts: 

DO  YOU  WISH  TO  COMPARE  THE  MEDIANS  OR 
VARIANCES  OF  THE  POPULATIONS?  ENTER:  1  TO  COMPARE 
MEDIANS;  2  TO  COMPARE  VARIANCES. 

Enter  1.   The  next  prompt  is: 

THE    NULL    HYPOTHESIS   STATES    -  THE 
MEDIANS   OF   X   AND   Y   ARE   EQUAL;   Mx   =   My.    WHICH 
ALTERNATIVE   DO   YOU   WISH   TO   TEST?     ENTER: 
1  FOR  HI :  Mx  <  My;   2  FOR  Mx  >  My ;   5  FOR  Mx  +    My. 

Enter  1  .. 
(4-)  Program  Interaction.   The  prompt  is: 

ENTER  X  DATA  (MORE  THAN  ONE  OBSERVATION  IS 
REQUIRED) . 
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Enter  the  X  data  separated  by  spaces.  The 
next  prompt  is: 

ENTER  Y  DATA. 

Enter  the  Y  data.  The  following  is 
displayed. 

THE  SUM  OF  THE  X  RANKS  IS:  224.  THE  U 
STATISTIC  EQUALS:   119. 

THE  P-VALUE  FOR  H0 :  Mx  =My  versus  HI: 
Mx  <  My  IS:   .5078. 

We  do  not  reject  the  hypothesis  of  equal 
population  medians  and  conclude  the  median  of  all  Army 
scores  is  equal  to  the  Navy's.   The  next  prompt  is: 

WOULD  YOU  LIKE  A  CONFIDENCE  INTERVAL  FOR 
THE  SHIFT  IN  LOCATION  (My  -  Mx)?  (Y/N). 

Enter  Y  (If  N  is  entered,  the  mainframe 
program  ends;  or,  the  Mann-Whitney  test  menu 
reappears).   The  next  prompt  is: 

ENTER  THE  DESIRED  CONFIDENCE  COEFFICIENT; 
FOR  EXAMPLE:  ENTER  95,  FOR  A  95^  CONFIDENCE  INTERVAL. 

Enter  95.   The  following  is  displayed. 

A  95%  CONFIDENCE  INTERVAL  FOR  THE  SHIFT  IN 
LOCATION  BETWEEN  POPULATIONS  X  AND  Y  IS: 

(-10   _<_  My-Mx  _<_  1  0  )  . 

The  mainframe  program  ends.  The  microcom- 
puter program  pauses  for  input  from  the  keyboard  by 
prompting: 
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PRESS  ENTER  WHEN  READY. 

Press  Enter  and  the  Mann-Whitney  test  menu 
reappears. 

d.  Mann-Whitney  Test  for  Equality  of  Variances 

(1 )  Description  of  Problem.-  Referring  to 
problem  3  in  section  c(1).  Is  there  sufficient 
evidence  to  claim  that  Army  scores  vary  more  than  Navy 
scores? 

(2)  Solution .  To  see  if  Army  scores  vary 
more,  we  test  H0 :  Vx  =  Vy  versus  HI :  Vx  >  Vy . 

( 3 )  Workspace  Decision  Process. 

(a)  Microcomputer:  Choose  the  Mann- 
Whitney  Test  from  the  main  menu,  and  the  option,  Test 
H0:  Vx  =  Vy  versus  HI :  Vx  >  Vy ,  from  the  test  menu  and 
receive  the  prompt: 

ENTER  THE  DIFFERENCE  OF  THE  MEANS  OR 
MEDIANS  (Mx  -  My ) . 

Because  we  believe  the  population 
medians  to  be  approximately  equal,  We  enter  0.  Skip  to 
the  Program  Interaction  section  below. 

(b)  Mainframe:  Enter  MANNWHIT  at  the 
keyboard  and  receive  the  prompts: 

DO  YOU  WISH  TO  COMPARE  THE  MEDIANS  OR 
VARIANCES  OF  THE  POPULATIONS?  ENTER:  1  TO  COMPARE 
MEDIANS;  2  TO  COMPARE  VARIANCES. 

Enter  2.   The  next  prompt  is: 
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THE  TEST  TO  COMPARE  VARIANCES,  REQUIRES 
THE  TWO  POPULATION  MEANS  OR  MEDIANS  TO  BE  EQUAL.  IF 
THEY  DIFFER  BY  A  KNOWN  AMOUNT,  THE  DATA  CAN  BE  ADJUSTED 
BEFORE  APPLYING  THE  TEST.  ENTER  THE  DIFFERENCE  OF 
MEDIANS  (Mx  -My)  OR  900  TO  QUIT. 

We  enter  0.   The  next  prompt  is: 
THE    NULL   HYPOTHESIS    STATES    -  THE 
VARIANCES   OF   X   AND  Y  ARE   EQUAL;   Vx   =   Vy .    WHICH 
ALTERNATIVE  DO  YOU  WISH  TO  TEST?   ENTER: 
1  FOR  HI  :  Vx  <  Vy;   2  FOR  Vx  >  Vy ;   3  FOR  Vx  jt  Vy . 

Enter  2. 
(4)  Program  Interaction.   The  prompt  is: 

ENTER  X  DATA  (MORE  THAN  ONE  OBSERVATION  IS 
REQUIRED) . 

Enter  the  X  data  separated  by  spaces.  The 
next  prompt  is: 

ENTER  Y  DATA. 

Enter  the  Y  data.  The  following  is 
displayed. 

THE  SUM  OF  THE  X  RANKS  IS:  166.  THE  U 
STATISTIC  EQUALS:   61  . 

THE  P-VALUE  FOR  H0 :  Vx  =  Vy  versus  HI: 
Vx  >  Vy  IS:    .0112. 

Consider  a  significance  level  of  .05. 
Since  a  P-value  of  .0112  is  less  than  .05,  we  reject 
the    null  hypothesis   of  equal  variances   in  favor   of 


65 


Vx  >  Vy  and  conclude  that  Army  scores  do  vary  more  than 
Navy  scores. 

The  mainframe  program  ends.  The 
microcomputer  program  pauses  for  input  from  the  key- 
board by  prompting: 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  the  Mann-Whitney  test  menu 
reappears . 

e.  Kruskal-Wallis  Test 

(1  )  Description  of  Problem  4.  During  a  recent 
Monster  Mash  involving  four  Navy  SEAL  Teams,  one  of  the 
events  consisted  of  the  number  of  pushups  a  man  could 
do  in  2  minutes.  Eight  men  were  chosen  randomly  from 
each  Team.   The  following  scores  were  recorded. 

SEAL  1:  90  96  102  85  65  77  88  70. 

SEAL  2:  64-79  99  95  87  74  69  97. 

SEAL  3:  101  66  93  89  71  60  76  98. 

SEAL  4:  72  78  73  81  83  92  94  86. 

Are  the  different  Seal  Teams  considered  to 
be  equally  fit? 

(2)  Solution .  To  see  if  the  Seal  Teams  are 
equally  fit,  we  test  the  hypothesis  that  all  the 
population  medians  are  equal. 
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( 3 )  Workspace  Decision  Process. 

(a)  Microcomputer:  Choose  the  Kruskal- 
Wallis  Test  from  the  main  menu;  and,  once  the  test  menu 
is  displayed,  press  Enter. 

(b)  Mainframe:  Enter  KRUSKAL  at  the 
keyboard. 

(4-)  Program  Interaction.   The  prompt  is: 

ENTER  THE  NUMBER  OF  POPULATIONS  TO  BE 
COMPARED  (MUST  BE  GREATER  THAN  TWO). 

Enter  4.   The  next  prompt  is: 

ENTER  YOUR  FIRST  SAMPLE. 

Enter  the  SEAL  1  data  separated  by  spaces. 
The  next  prompt  is: 

ENTER  YOUR  NEXT  SAMPLE. 

Enter  the  SEAL  2  data.  The  next  prompt  is: 

ENTER  YOUR  NEXT  SAMPLE. 

Enter  the  SEAL  3  data.  The  next  prompt  is: 

ENTER  YOUR  LAST  SAMPLE. 

Enter  the  SEAL  4  data.  The  following  is 
displayed. 

THE  H  STATISTIC  EQUALS:   .1335. 

THE  P-VALUE  FOR  H0 :  THE  POPULATION  MEDIANS 
ARE  EQUAL  versus  HI :  AT  LEAST  TWO  POPULATION  MEDIANS 
ARE  NOT  EQUAL  IS:   .98893. 
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We  do  not  reject  the  null  hypothesis  that 
the  population  medians  are  equal  and  conclude  that  the 
SEAL  Teams  are  equally  fit. 

The  mainframe  program  ends.  The  microcom- 
puter program  pauses  for  input  from  the  keyboard  by 
prompting: 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  the  Kruskal-Wallis  test 
menu  reappears. 

f.  Kendall's   B 

(1 )  Description  of  Problem  5.  In  order  to 
determine  if  cold  weather  affects  target  marksmanship. 
Naval  Special  Warfare  recorded  small  arms  marksmanship 
scores  and  corresponding  air  temperatures  for  a  period 
of  one  year.  20  men  were  chosen  at  random,  and  their 
scores  averaged  for  different  air  temperatures.  The 
average  score  for  each  air  temperature  is  shown  below. 

Air  temperature  (X):  50  55  20  50  65  55  30 
52  4-0  60. 

Average  scores  (Y):  210  200  165  165  260 
215  175  191  180  235. 

Can  it  be  said  that  colder  temperatures 
have  an  effect  on  marksmanship  scores?  Is  that  effect 
positive  or  negative? 
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(2)  Solution .  We  test  the  null  hypothesis 
that  no  association  exists  between  cold  temperatures 
and  marksmanship. 

( 3 )  Workspace  Decision  Process. 

(a)  Microcomputer:  Choose  Kendall's  B 
Test  from  the  main  menu;  and,  once  the  test  menu  is 
displayed,  press  Enter. 

(b)  Mainframe:  Enter  KENDALL  at  the 
keyboard. 

(4)  Prograjn  Interaction.   The  prompt  is: 

ENTER  X  DATA  (MORE  THAN  TWO  OBSERVATIONS 
ARE  REQUIRED) . 

Enter  the  X  data  separated  by  spaces.  The 
next  prompt  is: 

ENTER  Y  DATA  (NUMBER  OF  Y  ENTRIES  MUST 
EQUAL  NUMBER  OF  X  ENTRIES). 

Enter  the  Y  data.  The  following  is 
displayed. 

KENDALL'S  B  EQUALS:    .7817. 

THE  P-VALUE  FOR  H0 :  NO  ASSOCIATION  EXISTS 
versus:  HI:  DIRECT  ASSOCIATION  EXISTS  IS:   .00045. 

THE  P-VALUE  FOR  THE  TWO-SIDED  TEST  OF 
HYPOTHESIS  IS:  .0009. 

Since  the  P-value  for  the  one-sided  test 
equals  .00045,  we  reject  the  null  hypothesis  that  no 
association     exists     between     temperatures     and 
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marksmanship  in  favor  of  direct  association.  We 
conclude  that  colder  temperatures  tend  to  cause  lower 
marksmanship  scores. 

The  mainframe  program  ends.  The  microcom- 
puter program  pauses  for  input  from  the  keyboard  by 
prompting: 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  Kendall's  B  test  menu 
reappears . 

g.  Spearman's  R 

(1 )  Description   of   Problem  6.    When   fitness 

reports   are  written,   officers  of  the  same   grade   are 

ranked  against  each  other  based  upon  their  demonstrated 

level   of   performance.     Last   marking   period,    the 

Commanding   and  Executive  Officers  separately  ranked   9 

Ensigns  as  shown  below. 

Ensigns 
ABCDEFGHI 
CO(X):  64  1  528379 
XO  (Y)  :  56341  9728 

Does   any  association  exist  between  the  two 

sets  of  rankings? 

(2)  Solution .    We   test   the   null   hypothesis 
that  no  association  exists. 
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( 3 )  Workspace  Decision  Process. 

(a)  Microcomputer:  Choose  Spearman's  R 
Test  from  the  main  menu;  and,  once  the  test  menu  is 
displayed,  press  Enter. 

(b)  Mainframe:  Enter  SPEARMAN  at  the 
keyboard. 

(4)  Program  Interaction.   The  prompt  is: 

ENTER  X  DATA  (MORE  THAN  TWO  OBSERVATIONS 
ARE  REQUIRED) . 

Enter  the  X  data  separated  by  spaces.  The 
next  prompt  is: 

ENTER  Y  DATA  (NUMBER  OF  Y  ENTRIES  MUST 
EQUAL  NUMBER  OF  X  ENTRIES). 

Enter  the  Y  data.  The  following  is 
displayed. 

SPEARMAN'S  R  EQUALS:    .5500. 

THE  P-VALUE  FOR  H0 :  NO  ASSOCIATION  EXISTS 
versus:  HI:  DIRECT  ASSOCIATION  EXISTS  IS:     .0664. 

THE  P-VALUE  FOR  THE  TWO-SIDED  TEST  OF 
HYPOTHESIS  IS:     .1528. 

Consider  a  significance  level  of  .05. 
Since  a  P-value  of  .0664.  exceeds  .05,  we  lo  not  reject 
the  null  hypothesis  that  no  correspondence  exists 
between  the  two  sets  of  rankings. 
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The  mainframe  program  ends.  The  microcom- 
puter program  pauses  for  input  from  the  keyboard  by 
prompting: 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  the  Spearman's  R  test   menu 

reappears . 

h.  Nonparametric   Simple   Linear  Regression;  Least 
Squares 

(1)  Description  of  Problem  7.    Battery-powered 

Swimmer   Proplusion   Units  are  sometimes  used   to   aide 

swimmers   during  long  underwater  swims.    Recent   tests 

have   shown   that  a  nearly  linear   relationship   exists 

between   water   temperature  and  battery  life  for   these 

units.    The   following   17  data  points   were   randomly 

selected  from  the  test  results. 


/ater  temper 

ature 

(X) 

Batt 

ery  life  (Y) 

70 

3 

65 

2.75 

50 

1  .8 

40 

1  .2 

60 

2.4 

55 

- 

1  .9 

52 

1  .75 

50 

1  .7 

43 

1  .6 

40 

1  .  1 

72 

2.75 

55 

2 

48 

1  -5 

35 

.9 

70 

3.3 

68 

3 

57 

2.3 

(a)  Find  the  fitted  regression  equation 
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(b)  For  the   following  water   temperatures, 
predict  the  battery  life  of  the  units:  61  52  4-6  36.  . 

(c)  Can  we  determine  with  any  certainty 
if  the  slope  of  the  regression  line  equals  .05. 

(d)  What  range  of  values  could  be  used 
as  the  slope  of  the  estimated  equation  line  90%  of  the 
time? 

(2)  Solution.  To  determine  the  estimated 
regression  equation,  we  use  nonparametric  linear 
regression . 

( 3 )  Workspace  decision  process. 

(a)  Microcomputer:  Choose  Nonparametric 
Simple  Linear  Regression  from  the  main  menu;  and,  once 
the  test  menu  is  displayed,  press  Enter. 

(b)  Mainfrajne:  Enter  NPSLR  at  the 
keyboard. 

(4)  Program  Interaction.   The  prompt  is: 
ENTER   X   DATA  (MORE  THAN  TWO   OBSERVATIONS 

ARE  REQUIRED) . 

Enter  the  X  data  separated  by  spaces.  The 
next  prompt  is: 

ENTER  Y  DATA  (NUMBER  OF  Y  ENTRIES  MUST 
EQUAL  NUMBER  OF  X  ENTRIES). 

Enter  the  Y  data.  The  following  is 
displayed. 
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THE   LEAST   SQUARES   ESTIMATED    REGRESSION 
EQUATION  IS: 

Y  =   -1 .265   +  .060668X. 

The  next  prompt  is: 

DO  YOU  WISH  TO  ENTER  SOME  X  VALUES  TO  GET 
THE  PREDICTED  Y ' S?  (Y/N). 

Enter  Y  ( If  N  is  entered,  the  program  skips 
to  hypothesis  testing  for  the  slope). 

ENTER  X  VALUES. 

Enter  61  52  46  36.   The  next  prompt  is: 

THE  PREDICTED  Y  VALUES  ARE:  2.44  1.89 
1  .53  .92. 

WOULD  YOU  LIKE  TO  RUN  SOME  MORE  X  VALUES? 

Enter  N.   The  next  prompt  is: 

WOULD  YOU  LIKE  TO  TEST  HYPOTHESIS  ON  B,  THE 
SLOPE  OF  THE  EQUATION?  (Y/N). 

Enter  Y  ( If  N  is  entered,  the  program  skips 
to  confidence  interval  estimation).  The  next  prompt 
is : 

ENTER  THE  HYPOTHESIZED  SLOPE. 

Enter  .05.   The  following  is  displayed. 

SPEARMAN'S  R  EQUALS:    .5756. 

THE  P-VALUE  FOR  H0 :  B  =  .05  versus  HI :  B  > 
.05  IS:   .0079. 

THE  P-VALUE  FOR  THE  TWO-SIDED  TEST  OF 
HYPOTHESIS  IS:   .0158. 
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Consider  a  significance  level  of  .05.  Since 
a  P-value  of. 0079  is  less  than  .05,  we  reject  the  null 
hypothesis  that  B  =  .05  in  favor  of  B  >  .05,  and 
conclude  the  slope  of  the  regression  line  is  greater 
than  .05.   The  next  prompt  is: 

WOULD  YOU  LIKE  A  CONFIDENCE  INTERVAL  FOR 
THE  SLOPE?  (Y/N). 

Enter  Y  ( If  N  is  entered,  the  mainframe 
program  ends;  or,  the  Nonparametric  Regression  test 
menu  reappears).   The  next  prompt  is: 

ENTER  THE  DESIRED  CONFIDENCE  COEFFICIENT; 
FOR  EXAMPLE:  ENTER  95,  FOR  A  95^  CONFIDENCE  INTERVAL. 

Enter  90.   The  following  is  displayed. 

A  90^  CONFIDENCE  INTERVAL  FOR  B,  THE  SLOPE 
OF  THE  ESTIMATED  REGRESSION  LINE,  IS: 

(  .05333  1  B  <_  .07  ). 

If  the  estimated  slope,  does  not  lie  within 
the  confidence  interval,  the  following  would  be 
displayed. 

THE  LEAST  SQUARES  ESTIMATOR  OF  B  LIES 
OUTSIDE  THE  CONFIDENCE  INTERVAL.  DISCARD  THE  LEAST 
SQUARES  EQUATION  AND  USE: 

Y  =   -1.4458  +  .060833X. 

THIS  EQUATION  IS  BASED  ON  THE  MEDIANS  OF 
THE  X  AND  Y  DATA,  AND  THE  MEDIAN  OF  THE  TWO-POINT 
SLOPES  CALCULATED  FOR  THE  CONFIDENCE  INTERVAL  ON  B. 
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The  next  prompt  is: 

DO  YOU  WISH  TO  ENTER  SOME  X  VALUES  TO  GET 
THE  PREDICTED  Y'S  FROM  THE  NEW  EQUATION?  (Y/N). 

To  compare  results,  let  us  input  the 
temperatures  in  the  new  equation.  Enter  Y  ( If  N  is 
entered,  the  mainframe  program  ends;  or,  the 
Nonparametr ic  Regression  test  menu  reappears).  The 
next  prompt  is: 

ENTER  X  VALUES. 

Enter   61   52   46   36.    The   following   is 


displayed. 


1 .35  .74. 


THE   PREDICTED   Y  VALUES  ARE:   2.265    1.72 


The  next  prompt  is: 

WOULD  YOU  LIKE  TO  RUN  SOME  MORE  X  VALUES? 

Enter  N.   The  next  prompt  is: 

WOULD  YOU  LIKE  TO  TEST  HYPOTHESIS  ON  B,  THE 
SLOPE  OF  THE  EQUATION?  (Y/N). 

To  compare  results  once  again,  we  enter  Y 
( If  N  is  entered,  the  mainframe  program  ends;  or,  the 
Nonparametr ic  Regression  test  menu  reappears).  The 
next  prompt  is: 

ENTER  THE  HYPOTHESIZED  SLOPE. 

Enter  .05.   The  following  is  displayed. 

SPEARMAN'S  R  EQUALS:    .8287. 
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THE  P-VALUE  FOR  H0 :  B  =  .05  versus  HI :  B  > 
.05  IS:   .0000. 

THE  P-VALUE  FOR  THE  TWO-SIDED  TEST  OF 
HYPOTHESIS  IS:   .0000. 

The  mainframe  program  ends.  The  microcom- 
puter program  pauses  for  input  from  the  keyboard  by 
prompting: 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  the  Nonpar  ajnetric 
regression  test  menu  reappears. 
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APPENDIX  D 
MAIN  PROGRAM  LISTINGS  FOR  MICROCOMPUTER  WORKSPACE 

7  KEN:AiAA;B'.BX-.BJ;CiCX;CI:D',DD;DX;DI:DXI;S;POS;NEGiXX;II;N',DEN:NN'.NUM',P; 

11  fl    TBIS   fOnCTION   computes   TBE   KENDALL   B   it ATI STIC  WHICH  IS  A   MEASURE 
2 J  a    OF  ASSOCIATION  BETWEEN  SAMPLES.    P-VALUES   ARE   GIVEN   FOR   TESTING  ONE 
3]   fl  AND   TWO-SIDED   HYPOTHESIS   FOR   NO   ASSOCIATION   VERSUS   ASSOCIATION. 

u1   «  SUBPROGRAMS   CALLED   BY   THIS   FUNCTION  INCLUDE:    TIES,    TIESK,    KEN D ALP , 
5]   n  INTERP,    INPUT  AND   NORMCDF . 

7]  p  DISPLAY  TEST  MENU  AND  INPUT  DATA 

8]  Nl:S-<-MENU   KEN03J 
9j    --(S  =  1J/B1 
10]  MENU  MAINOBJ 
11.   *0 

12  Bl:R<-INPUT  2 
13.   Ofl+i? 

m.   X-«-l  +  (Q+l)+5 
15]  Y-(Q+l)rR 
.16 j  0  ORDER   Y  IN  INCREASING  ORDER   OF  X 

;17,   4-t-ICAX: 

rial  fl  ORDER   X  IN  INCREASING  ORDER 

3-t-XZilXl 
p  COMPUTE  CURRENT  RANKING  OF  I 

C-*-iAA 
p  NOW  ORDER  Y  RANKS  IN  INCREASING  ORDER 

Di-AtAAl 
fl      IF  TIES  EXIST  IN  EITHER  X  OR   Y  RANKED    VECTOR   USE  MID -RANK  METHOD 
DD-^1   TIES   D 
XX-ol  TIES   3 
fl        FIND  ORIGINAL  RANKING  OF  Y  WITH  TIES  RESOLVED 
YY-DDLCl 
N*oX 
fl  COMPUTE  NUMBER  OF  DISTINGUISHABLE  PAIRS  V 

NN-^(.NxiN-m*2 
S-^pO 
AA-t-0 
p  POSITIVE  ONES  COME  FROM  A  RUNS  UP  CONDITION:    NEGATIVE   1  FROM  RUNS  DOWN 
p  ZERO  IS  SCORED  FOR  TIES.   MULTIPLY  THE  RESULTS  FOR  EVERY  ELEMENT  AND  SUM 
Ll:AA-^AA+l 


L19 

20 
21 
22 

23 
2«- 
125 
126 
.27 
.28 

r29 

.30 
31 
32 

33 

;34 

C35 

.36 

37 

38 

39 

.UO 

.U2 
U3 
.4U 
>5 
>5 

;u7 

.U8 
.49 
.50 

51 
.52 

53 

5'+ 

55 

56 

57 

58 

IS' 


BX'>-(,XX[,AA1>(,AA^XX}) 
CX-«-(XXCM3<UA+XX))x(  1) 
DX*BX+CX 
BY-*-(,YYZAAl>(AA^YY)) 

cy*IyyIaa2<Ua^yy)U(.   1) 

DY->-BY+CY 

DXYirDXy.DY 

POS-<r(,DXY>0} 

NEG'>-lDXY<oUri) 

S-^S.POS,NEG 

■>-UJi<(N-l))/Ll 
p  SUM  FINAL   VECTOR  TO   DETERMINE  S 

S-<-+/S 
p    OBTAIN  TEE  NUMBER  OF  TIES  IN  EACH  VECTOR  USING  THE  TIESK  FUNCTION 

U*TIESK  B 

Vi-TIESK  D 


sui-+/(2:u) 

SV->-+/l2'.V) 


p  CALCULATE  THE  B   STATISTIC  INCLUDING  THE  CORRECTION  FOR  TIES 
r+S*  ( (.NN-SU)^  (NN-SV)  )*0  .  5 
AT->r\T 

*(.N>13)/N0RM 
p         CALL   KENDALP  TO   CALCULATE  THE  RIGHT  TAIL  OF  THE  CDF  OF  B 

P*KENDALP  N 
fl  CALL   INTERP   TO   CALCULATE  P-VALUE   BY   INTERPOLATION 
PVAL->-AT   INTERP   ? 
,.^^      *{pVAL~~\)/L2 
C6Hj  PVAL*0.i 

!55l  *z;3 

'%%\  lrPA^^.J^k^'^,1  lJ.^^^  USING  NORMAL  APPROX. 

.6  71  NORM:NUM-(3xAT)x((i2xNN)i'0.5) 

.68]  DEN*-(2x(l2xN)  +  5))*0.5 

.6  9]  Z*NUM*DEN 

.70]  PVAL-^l- (NORMCDF   Z) 

.711  ^IFB   IS  POSITIVE  PRINT  OUT  DIRECT  ASSOCIATION, 

.72J  Lr3  r •*■  CT>0  J /L 5 

,73]  CBA->-<  INDIRECT^ 
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;7U]  -..^7 

7  =  1  Li-.CBA-^' DIRECT^ 

"A  L7  •.PV*2^PVAL 

79]  L8: ' KENDALL' 'S   B   EQUALS:     ' ,(^vT) .UTCNL 
"sol   <TBE  P -VALUE  FOR   BO:    NO  ASSOCIATION  EXISTS   VERSOS' 

81]   '  Hi:  <,(fCBA),<    ASSOCIATION-FXJSTS  IS:    '  ,  (UsPVAi )  .DTCiVL 

82]   'TEE  P-VALOE  FOR  TBE  TWO-SIDED  TEST  OP  BIPOTBESIS  IS:     ' ,l^9PV) ,QTCNL 
.33]  'PF.ESS  ENTER   WBEN  READI .  ' 

SuJ  WW->-^ 

!35j  */Vl 
7 


^ 


a 


7  KRVLiNUy!:DENOM;A:C;B',B',KiAA:BB:DD:E',F;NiOF',P:PVAL:R:SOFR:SR:TSOR:CBA:B 
>     TH^S   FUNCTION   COMPUTES   THE   KUSKAZ-WALLIS   TEST   STATISTIC   H   '^HlCd   IS 
p  A   MEASURE   OF  TSE  EQUALITY   OF   K   INDEPENDENT   SAMPLES. 

n  SUBPROGRAMS   CALLED   BJ   TBIS   FUNCTION  INCLUDE:    TIES,    TIESK,    INDEXPLS, 
q  FDISTN,    INTER?  AND   TBE   VARIABLES  PMAT KW 20 ,    PMATKW31 .    PMATKW3  3 ,    PMATKN3n 
;.  =  .   4  PMATKW^l.    PMATKWai,    AND   PMATKW>^3 . 
6 

17;   p    MENU  CBOICES  AND  ROUTE  TO  PROPER  STATEMET  FOR  ACTION. 
8    Nl:B*MENU  KRWL03J 

9]         -►(S=1)/S1  

^^--      liENU  MAINOBJ 


1 1  I 


2]  <2P'°*-^ 

:i3]  3l:' ENTER  TBS  NUMBER  OF  POPULATIONS   TO   BE   COMPARED    (.MUST  BE   GREATER   TBAN  T 
WO).  < 

;i5l  ^((jr<3)v(a  ix)«o))/ffi 

16]  *UpK)>l)/El 

17]  p               INITIALIZE  VECTORS  E  AND  F  AND   VARIABLE  C 

18]  E-^-F-SOFR-^-oO 

,19]  C->-Q 

20]  nTBIS   LOOP   FACILITATES  ENTERING  TBE  SAMPLE   VECTORS  AND   STORING  TEEM 

.21 1  Ca  A"- 'FIRST' 

^22  il:C*<:+l 

23  'ENTER  lOUR    ',(«CflA),'  SAMPLE.' 

2**,  D-Q 

25  *((po£5)sO)/iVEXr 

.26.  iJ-«-lpi? 

,27  p    CONCATENATE  SAMPLES  AS  TBE J  ARE  ENTERED  AND  STORE  TEEM  IN  VECTOR  E 

28  NEXT:E-<rE,D 

29,  P         RECORD   TBE  LENCTES  OF  TEE  SAMPLES  AS  TEE J  ARE  ENTERED 

30  F*F,^D 

^31  CEA-t-'NEXT' 

'^32  *(C<(K-1))/L1 

l33,  CEA-'LAST' 

2u  ■*(.C<K)  /  LI 

,35;  p       RECORD   SIZE  OF  ALL   SAMPLES  WEEN  COMBINED 

36  N--+/F 

3  7  p  ORDER   SAMPLE  SIZES  LARGEST  TO   SMALLEST 

38  0F*FC7F] 

39  P  ORDER   COMBINED  SAMPLE  VECTOR  TO  BE  USED  BI  TIES  FUNCTION 
>0  D-hECAT] 

>1  p    CALL   INDEXPLS  TO  INCREMENT  INDEXES  WBEN  TIES  OCCUR  WITEIN  ONE  SAMPLE 

>2  AA*F  INDEXPLS   E 

>3  p           CALL   TIES  TO   BREAK  TIES  BI  MI DRANK  METBOD 

uu  BB*1   TIES  D 

,^5,  C*0 

.^6.  a                           TBIS  LOOP  CALCULATES  TBE  E  STATISTIC 

.H7.  L2:C*C+1 

•!*f-  "oT,  >r,r>r  ,^r^^^  ,'^^ M^^S.JO^  ^^^^  SAMPLE  IS  CALCULATED 

,50,  ^ cRi\Wi¥^1lr^^^  °^  RANKS  SQUARED   DIVIDED   BI  TBE  INDIVIDUAL   SAMPLE  SIZE 

!52;  0                                 STORE  EACB   CALCULATION 

,53.  SOFR-*-SOFR,SR 

.5U.  '-(.C<K)/L2 

,55]  p             SUM  ACROSS  ALL   SAMPLES 

.56^  TSOR*-*/SOFR 

.5  7]  a             CALCULATE   FINAL   3   STATISTIC 

.38j  3*(TS0R^  ( 12'!-  (Nx  (N+1  .  }  }  )-'  3x(7/+l  ]  ; 

.59,  P  RECALCULATE  B  WITB   CORRECTION  FOR  TIES 

,60,  A-*-TIESK  E 

.61,  iVra-<-(+/(.4*3))-(+/A) 

62  DEN0M-<rN^((N*2)-l) 

.6  3,  B-'-Biil-CNUMiDENOM)) 

,6u.  p  S7STEM  OF  LOGICAL   STATEMENTS  ENSURE  PROPER  PROB .    IS  ACCESSED                       '  ' 

,65,  •*-<,OF[.l2<2) /OUTPUT 

,66.  *lK>*i)/FAPPROX 

,67.  ■>'Ik=3)/IF 

.6  8.  *i<<^/AOP=    3  3  3  3))v(0F{:11>3))/FAPPR0X 

.69.  -►(A/(OF=  2  1  1  1))/0UTPUT 

.70]  ♦((i:='+)A(C7FCiJ=2))/PU2 
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[71]       ■^{{^/(,0F-    3'1    1    l))v(A/(OF=    3    2    1    1 )  )v  (A/COFr    3    3    1    1))v(a/(0P=    3    3    2    1)) 

)/PUl 

-(J(:=u)/Pii3    ^ 
If : -»■  ( OF  C 1  ]  »+ ) /FAPPi?OX 

->-((a/(OF=    2    1    1))v(a/(0F=    3    1    1)))/0UTPUT 

*((OFCl]='i)A(CFC3]=l)}/P31 

*<(A/Mir=    3    2    1}}v(a/(0P=    3    3    1))v(a/(0F=    3    3    2))v(A/(aF=    3    3    3)))/P33 

*((a/(OF=    3    2    2))v(OFCl]=n))/P3U 
p  CALL  APPROPRIATE   VARIABLE  ACCESS  CDF 

P22:P'^PMATKW20L{N-^)'.  ;] 

-PAf 

P3 1 ;  P*-PMATKWZ  I  C  ( A?-  5  ) ;  :  ] 

-^PM 




P3  3  :P-<-PMiiri?W3  3  C  (A7-5  ) 

-»-PM 
P3  u  :  P-t-PAMriny3'+  C  (A?- 6  ) 

•*PM 
Pui  :P-<-PMAria'Ul  C  (iV-5  ) 

-»-p>f 

PU2  :  P*PMrXWii2  C  (W-  5  ) 
P^Z:P-<rPMATKWn2  [  (iV-7  ) ; 


_  .  _  _   ___._   _ ] 

p     CALL   INTER?   TO   CALCULATE  P -VALUE  31   INTERPOLATION 

'  ->-(?VAL='l)/ OUTPUT 

p  CALCULATE  P -VALUE  USING  TEE  F  VIST  W/ONE  LESS  D.F.IN  DENOM  APPROX. 
FAPPROX:F^((N-K)^H)T(iK-l)^((N-l)-H}) 
PVAL-<rl-(.(.lK-l),(.(.N-K)-l'))   FDISTN  F) 

*L5 

'  1 0  0  ]  OUTPUT :  PVAL-^ '  GREATER   THAN    .25' 
tlCi:  LB:  'THE   H   STATISTIC   EQUALS:     ■  .  (^'bH)  .UTCNL 

[102]   'T3E  P-VALUE   FOR   HO:    THE   POPULATION   MEDIANS   ARE   EQUAL    VERSUS    ^  ,3TCNL 
riOS]   '       31:  AT   LEAST  TWO   POPULATION   MEDIANS   ARE  NOT   EQUJiZ~~B~    '  ,  (usP'/^L) , 
'^'TCVL 

'''press  enter  when  READI.' 

AA-^rQ 

->-Nl 
El: 'ERROR:    ENTER  A   SINGLE  INTEGER   VALUE  GREATER  THAN   2;  TRI  AGAIN .' .DTCNL 

*S1 
7 


7  MANWiN',M:PV2:A:B',C;G;MM:NN:RX',U;Ny!l;PiNU',PVAL;NK',NUMZ;NUMZl:DEN:DENCiDE 
NCI :  TCiTCl;NUMiZ;Zl;  ALPHA;  CDF  iINDEX -.IPX  iCI;UALPHA;BB:CC;Ul;U2;PV;NNl;NN 
2iPVI;DIFFiAX:AY;AA;GG:AXl',PVpl;PV3:b;Q:R 
fi    THIS  FUNCTION   USES  THE  SUM  OF  RANKS  PROCEDURE  TO   CALCULATE  THE  MANN- 
P  WHITNEY   U  STATISTIC   WHICH   IS   USED   IN   COMPUTING  THE  P-VALUE  FOR  THE   TEST 
P  OF  LOCATION  AND   SCALE.      THE      C.I.    FOR    (MI-MX) ,    THE   SHIFT  IN  LOCATION   IS 
P  ALSO   COMPUTED.    SUBPROGRAMS   CALLED   BY   THIS   FUNCTION  INCLUDE:    TIES,    TIES2 
P  INDEXPLS,    VARMW,    MANWP ,    INPUT,    CONFMW ,    NORMCDF ,    AND   NORMPTB. 

■*-N3 
Nl:MENU  MANWHELP 
p  MENU  CHOICES'JND  ROUTE  TO  PROPER  STATEMENTS  FOR  ACTION, 


,91      N3:D-f-CH0ICEM  PACEDMENU  MANWOBJ 


.12 

13 

.l** 
15 
.15 
.17 
.18 

C20 
.21 
.22 
.23 
,2U 
.25 
.26 
.27 
28 

rso 

'31 
l32 
.33 
.3U 
35 
.36 
.37 
.38 
.39 
.40 

^2 

.U3 

.44 


£iJFF-«-0 

-^(Z3  =  2,3,4,5)/51 

*(£i  =  6,7.8)/S3 

-^(n  =  l)/A'l 

MENU  MAINOBJ 

*0 
B3:' ENTER  THE  DIFFERENCE  OF  THE  MEANS  OR  MEDIANS    (MX   -  MY).' 

DIFF-^n 

-►((pCJFF)>l)/E3 
p  ENTER   DATA    VECTORS 

Bl:' ENTER  X  DATA    {MORE  THAN  ONE  OBSERVATION  IS  REQUIRED).' 

N-<-a 

■>-(XpoN)  =  0)/El 
< ENTER   I   DATA . ' 

p  IF  CALCULATIONS  INVOLVE  VARIANCES  ADJUST  X  BY  TEE  DIFFERENCE  IN  MEANS 

N-^N-DIFF 
p         CONCATENATE  X  AND  Y  SAMPLE  VECTORS 

A*N,M 
g  DETERMINE   SIZE   OF    'A   AND   I   VECTORS   AND   ASSIGN   TO   NN   AND   MM 

m-^oN 

MM*oM 
p         COMPUTE  SIZE  LIMIT  OF  LEFT  TAIL   OF  NULL  DISTRIBUTION 

NM*(NNyMM)*2 

NMl^lNM 
P  ORDER  A   AND   ASSIGN  TO   B 

5*ACM] 

C-*- Inn, MM)    INDEXPLS  A 
p      CALL   TIES  FUNCTION  TO  BREAK  TIES  USING  MIDRANK  METHOD 

C'*'!  TIES   B 
p  IF  FALSE   CALCULATE  TEST  FOR   VARIANCES 

-»-(Ii  =  2,3,4,5)/S5 
p      CALL    VARMW  TO   GENERATE  RANKS  REQUIRED   FOR   VARIANCE  TEST 

GG-VARMWiNN+MM) 
a    CALL   TIES   TO  RECORD   TIES  IN  THE  DATA   AND   BREAK  TIES  IN  GG 


80 


>5 

y7 

>8 
H9 

50 
51 
52 
53 
5>+ 
!55 
,56 


Ct-CG  TIES   3 
p  CALCULATE  SUM  OF  X  RANKS  ■     " 

35:EX-<-+/(.CZ(.Ctl;(.\NN)2)l) 
a  CONVERT  TO  MANNWBIT   U  STATISTIC 

y-f-HX- (  (iViVx  (JVW+1) )  * 2  ) 

A  ^"^    IF  SIZE  OF  X  TIMES  SIZE  OF  I   >   80;    GO  TO  NORMAL  APPROX 

*((JVMx2)>80)/Ii2 

NNl*i/NN,MM 

NN2-*-\  /NN  ,MM 
aMANWP   FUNCTION  CALCULATES   LEFT  TAIL   CUMULATIVE  PROBS .    OF   U  STATISTIC 

P-^NNl    MANVP  NNl 


■§8'  fl      LOGICAL   STATEMENT  ENSURES   ONLI  LEFT  SIDE   OF  NULL   DIST  IS   USED 

"59"  •>-CU^NMl)/L3 

■50'  P       CONVERT   U  STAT  WHEN   GREATER   THAN  LEFT  TAIL    VALUES 

'61'  a->-(.NNxMM)-U 

■62  a  IF   Ul   IS  A   FRACTIONAL.    INTERPOLATE  P   VALUE 

'631  L3:*((.1\U)  =  0)/NON 

51^  y2-f-i  U 

[65]  *{U2>G)/?1 

r57l  *?3 

'-"^    P3:>7J*(i'CC/2  +  li+>C£/2  +  2j  )*2' 

■*CHECK 
NON:*(.U>Q')/GO 

P7+1 

*P2 
GOiPV-l-Pt^'i 
PliPVI^P^lU-^Dl 
CHECK :*iUliNMl)/L^ 

pv-^pvr 

PVI^PVl 

fL^ 
p  IF   CONFIDENCE  INTERVAL   DESIRED   GO   TO   LIO 
L2:->-(D=5  )/LlO 
a  COMPUTE   THE     NORMAL   APPROXIMATION  W/CORRECTION   FACTOR 

NUMZ*(U+0.5)-NM 

NUMZl-^iU-O.S'j-NM 

DEN-^'l  [MMxNNx(.MM+NN-^l)  )  +  12}*0.3 

Z'>-NUMZrDEN 

Z1->-NUMZ1tDEN 

NUMf-\  CU-NM) 

DENC-^li  (  InN+MM-1  )x  (DEN*2  )  )*(NN+MM-2  ))-(( (NUM-O  .  5  )*2  )+  (.NN+MM-2) )  )*0  .  5 

DENCl-f-ll  ( (.NN+MM-1  )x  (.DEN*2  ) )+  (.NN+MM-2  ) )-  U  (JVi/M+0  .  5  )*2  )*iNN+MM-2  )  )  )*0  .  5 

IC*{NUM-Q.5)+DENC 

TC1--(NUM+Q  .  5  )iDENCl 

*(UiNM)/ SECOND 

PVI*((.NORMCDF   Z)+(.  (NN+MM-2)   TDISTN  TC))*2 

PV-<-((l-(NORMCDF   Zl) )  +  (!-( (NN+MM-2)   TDISTN  TC1)))*2 

*Ln 
SECOND  :PVI-^((NORMCDF   Z)  +  (l- ((NN+MM-2)   TDISTN  TC)))*2 

PV-*-((\-(NORMCDF   Zl))  +  (  (NN+MM-2)   TDISTN  TC1))*2 

Ln:PV3-2x ( L / (PV,PVI ) ) 
*(PV3£1)/N5 
PV3*-1 

N5:PVM-*-(3.1)p  (PVI,PV,PV3  ) 
'THE  SUM  OF  THE  X  RANKS   IS:     ',(9RX),'.      THE   U  STATISTIC  EQUALS:     ',(«£/!), 

nrcNL 

LOGICAL   STATEMENT  FOR   VARIANCE  OUTPUT 


[105]    P 

106J      *(D  =  6^7_,8)/VAR 
!l07]       'THE  ^ 


___    P-i/ALUE  FOR  so:    MX   =   Ml  VERSUS  Hi:    MX    '  ,  (uLOCICCZJ-l ;  1]  ) ,  '    MY  IS:    ' 

(5  i*  9PVMID-1;12),DTCNZ  

•"LS 
VAR:PVM*(3 ,l)p (PV ,PVI ,PV3) 

'THE  P-VALUE  FOR   BO:    VX   =  V7   VERSUS  HI:    VX    ' , (9L0GICLD-5:ll ) , '    V7  IS:     ' 

(6  u  9PVM[D-5;ll),aTCNZ  

•  PRESS   ENTER  WHEN  READY . ' 

BS-f-D 

*N3 
L8: 'WOULD   lOU  LIKE  A   CONFIDENCE  INTERVAL   FOR  THE  SHIFT  IN  LOCATION    (Ml   - 

MX)?    (I/N).' ,DTCNL 

, ~'3B='N'  )/N3 

^117]  LlOiCC^iNPUT    3 

.113]     <;LPifA-' ioo-crH2  0o 

113]  a     ROUTE   TO   NORMAL   APPROX.    FOR   CONF .    INT.    OF  LARGER   SAMPLE   SIZES 

120]  ■*-((NM^2)>80)/L5 

121]  p  COMPUTING  CONFIDENCE  INTERVALS  BY   EXACT  P-VALUE 

122]  CDF*P 

123]  p        INDEX  POSITION  OF   VALUE  IN   CDF  i,   ALPHA 

124]  INDEX*- (  +  /(CDF<.ALPHA)) 

125]  ->-(INDEX>0)/L6 

125]  INDEX*-! 

127]   *L5 

.128]  p         COMPUTING   CONFIDENCE  INTERVALS   USING  NORMAL   APPROX.    W/C.F. 

.129]  L5:DEN*-((MMxNN^(MM+NN+l)  )*12)*0.5 

130]  UALPHA*(DENy<-(NORMPTH   ALPHA))+NM-O.S 

.131]  p  ROUND   U ALPHA   DOWN  AND   INCREMENT  BY  ONE 
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MX   5    ',(«CI[2]),'     )<  ,UTCNL 


;i3  2]       INDEX-*-lUALPHA*l 
;i33]    L&:IFX*-NN^INDEX 

ci*:px  confmw  a 

'  a  '  ,(.^c0),'     confidence  interval  for  the  shift  in   ' 

'  location  between  populations  x  and  i  is'.<.  qtcnl 

;i37l   '  (  ' ,{9CIZ12  ),'  ^   MY 

138]  ^ PRESS   ENTER  WEEN  READY.' 
139J  BB->-n  ^ 

;i4l1  El: 'ERROR:    SAMPLE  CONTAINS  LESS  THAN  TWO   ENTRIES'.    TRY  AGAIN .'  .DTCNL 

■iu33  E3:' ERROR:    YOU  RAVE  ENTERED  MORE  TBAN  ONE  VALUE.    TRY  AGAIN .' .UTCNL 
.144]   -»-33 
7 


V  NPLR;N;SUMX: SUMY ;XBAR: YEAR; SUMX2',SUMXY;B:A',WW:XX-,BB;U:D',ALPHA;PiCC;CDF: 

T ALPHA  :  NN :  Cl :  SLOPES ;  RR :  3R ;  I'J :  DENOM :  INDEX  iFF-.X:!  :Q:R;  CHA  ',E;?V 
fl      PROGRAM   CONDUCTS   NONPARAMETRIC   ZINEAR   REGRESSION.    IRE   ZEAST   SQUARES 
fi  ESTIMATED   REGRESSION   LINE   IS   COMPUTED   WITH   HYPOTHESIS   TESTING  AND 
fl  CONFIDENCE   INTERVAL   AVAILABLE   FOR   THE   SLOPE   3.    IF   3   DOES   NOT   LIE   IN   THE 
fl  C.I.    .-^N  .-ALTERNATE   .^.EGRESSION   ZINE   IS   PROPOSED.    SUBPROGRAMS   CALLED   ARE 
0  SPMANP,    KENDALP,    NORMPTH .    INPUT,    AND   CONFLR. 

[lPP-f-5 
p  DISPLAY   MENU  AND   INPUT  DATA. 
Nl:EirMENU  NPLROBJ 

*(E=1)/N2  

MENU   MAI NO 3 J 

->-Q 
N2:RirINPUT   2 

Qf-1*-R 

Z'^l^iQ  +  D'i'R 

J-»-(fl+l)+i? 


N'<rpX 

SUMXi-*/X 
SUMY*+/Y 

X3AR*SUMX*N 
YBAR-^SUMYiN 

SUMX2->-+/(.X*2) 


ASSIGN  TEE  SIZE  OF  X 
COMPUTE   THE   SUM  OF 


(.AND   Y)   TO  N 
X'S  AND   Y'S 


COMPUTE  THE  MEAN  OF   X   AND   Y 

COMPUTE  THE  SUM  OF  THE  X'S  SQUARED 
COMPUTE  TEE  SUM  OF  X  TIMES  I 


SUMXY-^+ZCX^Y) 
fl  COMPUTE    • S •  THE  SLOPE  OF  THE  ESTIMATED   LEAST  SQUARES  REGRESSION  LINE 

B-*-(.  (.N^SUMXY )-lSUMXxSUMY )  )*  ( (.NxSUMX2  )-(.SUMX*2)) 
fl  COMPUTE    <A' ,    TEE  Y -INTERCEPT 

A-*-YBAR-(.By(XBAR) 

FF* ' N ' 

'THE  LEAST  SQUARES   ESTIMATED  REGRESSION  EQUATION  IS'.'.UTCNL 

'7=  ',(.59A).'    +  •  (5Bfl).'X.  '  .DTCWL 

'DO   YOU  WISH  TO  ENTER  SOME  X   VALUES  TO  GET  THE  PREDICTED   7 "5?  (J/W).' 

wv->-a 

*(.WW='N'  )/Ll 
L2: 'ENTER  X   VALUES. • 

XX-^D 
p  CALCULATE  PREDICTED   Y'S 

YY^A-i-BxXX 

'THE  PREDICTED   Y   VALUES  ARE:     '  .(.'sYY). QTCNL 

'    WOULD   YOU  LIKE  TO  RUN  SOME  MORE  X   VALUES?    (.Y/N).' 

ww*-^ 

^(WW='Y' )/L2 
LI : 'WOULD   YOU  LIKE  TO   TEST  HYPOTHESIS   ON  B,    THE  SLOPE  OF  THE  EQUATION?    (.Y/ 
N).  ' 

ww-i-a 

^(.WW='N'  )/L3 

'ENTER  THE  HYPOTHESIZED  SLOPE. ^ 

p  COMPUTE  UI'S 

U--Y-(.BBy(X) 

CHA'>-<>    ' 
p  CALL   SPMANP  TO   COMPUTE  R  AND  ASSOCIATED  P-VALUES. 

D*:i   JPMANP    -J 

*(Dii:>0)/Lll 

CHA^'<    ' 
Lll:PV'>-2»DL2l 

■*(PV^1)/L19 

PV*1 
L19:  'SPEARMAN'  'S  R   EQUALS:     '  .(nvD^l"]  )  .DTCNL 

'THE  P-VALUE   FOR   HO:    B    =    ',(«SS),'  VERSUS' 

Hi:    B    ',(9CHAil    237.T5BB)  •  IS: 


(i^9DL2l  ).nTCNL 

-         A^9l2xDL2l)),DTCN 


'THE  P-VALUE  FOR  THE  TWO-SIDED   TEST  bP  HYPblHESIS   Is 

L 

p     IF   USING  THE  NEW  REGRESSION  EQUATION  BASED   ON  MEDIANS,    EXIT  HERE. 
L3:*(FF='N' )/LlS 

< PRESS   ENTER  WHEN  READY . ' 

ww-'-a 
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59]  -"Nl 

7  0   0  COMPUTE   CONFIDENCE  INTERVALS  ON  B  .  ,  , 

.71  LIS: 'WOULD   YOU  LIKE  A   CONFIDENCE  INTERVAL   FOR  TEE^  SLOPE?    (.I/N).' 

72  vw-i-a 

73  •*(WW='Nn/Nl 

7u  LlQ:CC-*-INPUT    z 

7  5  fi        CHANCE  ENTERED   VALUE  TO  ALPHA 

76    ALPBA-f-(100-CC)*200 

^77  fl    ROUTE  TO  NORMAL  APPROX.    FOR  CONF .    INT.    OF  LARGER  SAMPLE  SIZES 

78    *(.N>13)/L5 

7  9  fl  COMPUTING   CONFIDENCE  INTERVALS  BY  EXACT  P-VALUE 

.80]  P-^KENDALP  N 

81j  CDF*P12:: 

"8 2 J  n        INDEX  POSITION  OF  VALUE  IN  CDF   5  ALPHA 

8u1      *(JiVDrX>0)/I5 

85  J      INDEX->rl 

86  J       -►LS 

87]  p        COMPUTING  CONFIDENCE  INTERVALS  USING  NORMAL  APPROX.    W/C.F. 

88]  L5:DEN0M*({N^(N-1)^(  (2xA?)  +  5))t18)*0.5 

[8  9]  r^iSPa^-f-JiJiVOMx  (  |  iNORMPTH   ALPHA  }  ) 

l901  *L~ 

C  9 1-3  L6  :  r^LPff4-»-?  C  3  ;  INDEXl 

:92]  Z-AIPHA 

;93i  L^:CI^X   CONFLR   I 

91*  j  NN-I^CI 

9  5]  SZ;CPES-<-l+CT 

96]  RRi-i<.(.NN-TALPHA')*2) 

97]  AS 

98]  *(J?i?sO)/I20 
ragi  ^R->-^ 

f i 0 5 ]  i:2 0 ?si;-i-r (i + (  (.nn-^talpha  } ■;•  2 } ) 
[101I  Si? 

^102]  ^(SR^(zSLCPES))/L21 
"""'   5P-<-oSLC;p£'S 

i21:'A  ' , (SCO , '  CONFIDENCE  INTERVAL   FOR  B.    THE  SLOPE  OF    ' 

'  THE  ESTIMATED   REGRESSION  LINE ,    IS:' , DTCNL 

'  (  <(.5sSLOPESlRR2).'    S5^  ^  .  (S'sSLOPESZSRl) ,  '     ).',UTCNL 

' PRESS   ENTER  WHEN  READY . ' 

q  IF^3   OUTSIDE   THE  C.I.    CALCULATE  NEW  EQUATION   BASED   ON  MEDIANS 

-( iS^SLOPESZR?.:  ja(3^SL0PESZ3R2  )  )/Nl 
fl  ORDER   X  AND   I 

fl     CHECK  TO   SEE  IF  THE  SIZE  OF  SS  IS  EVEN  OR   ODD   FOR  FINDING  MEDIANS 

*(.(.2\NN)  =  0)/S1 
a  COMPUTE  MEDIAN  FOR  ODD   CASE 

B-^-SLOPESL  ( (JVW+l  )-i-2  )] 

*S2 
p  COMPUTE  MEDIAN  FOR  EVEN  CASE 

Sl:B-*-(.SLOPESZ(NNi2n+SLOPES[.(.(NN+2)T2)l  )t2 
p  DO  THE  SAME  FOR  THE  X  AND  Y  VECTORS 
S2:-((2Ta^)  =  0)/S3 


03 

lOH 

105 

106 

107, 

108; 

109 

llGj 

ll"*  " 

.112 

!ll3 

liu 

C115 

*^116 

117 

118 

119 

120 

121 

122 

;i23 

.12U 

125 

!l26 

:i27 

^128 

129 

130 

131 

132 

133 


YBAR'*-Yt(,(,N+l)i2U 
XBAR-rXlllN+l)*2n 


*OUT 
S2:YBAR*(Yi:(N*2  )]  +IC  ( (.N+2  )  +  2  )]  )  +  2 

XBA;?*(XC(iV*2)]+X[((iV+2)  +  2)]  )  +  2 
p  COMPUTE  NEW  INTERCEPT    'A' 

OaT:A*YBAR-{B^XBAR) 

'THE  LEAST   SQUARES   ESTIMATOR   OF  B   LIES   OUTSIDE  THE  CONFIDENCE  INTERVAL.' 

' DISCARD   THE  LEAST  SQUARES  EQUATION  AND   USE : ' , DTCNL 

'  Y   =  '  .(fA),<    +  ',(98), 'X    '.DTCNL 

'  THIS  EQUATION  IS  BASED   ON  THE  MEDIANS   OF  THE  X  AND   Y  DATA. AND   TH 

E' 
C13U]   'MEDIAN  OF  THE  TWO -POINT  SLOPES   CALCULATED   FOR  THE  CONFIDENCE  INTERVAL   0 

N  B.' .DTCNL 

ALLOW   USER  TO   DO   SOME  ANALYSIS   ON  NEW  EQUATION 

'  DO   YOU  WISH   TO   ENTER   SOME  X   VALUES  TO   GET  PREDICTED   Y"S    ' 

'FROM  THE  NEW  EQUATION?    (.Y/N).' 

PP-^Cl 

■>-(FF='Y'  )/L2 

pp^ ' Y ' 

*L1 


7  SIGN;AiC:B:DiPVAL:XiMO;N:CDF;ALPHA:CIiY',AA'.BBiCC;DD;PV'.PVI'.NNN;KPOS;ORD 
DiYY:0kDXikAtPHA:Z:Zl;QUA:WW;PVM'.Pir3ik;d 
P       THIS  FUNCTION   USES   THE  ORDINARY   SIGN   TEST  TO   CALCULATE  THE  K 
p  STATISTIC,    P-VALUE.    AND   CONFIDENCE  INTERVAL   AS  A   TEST   FOR   MEDIANS. 
p  THE  LAST  OPTION  WILL   DISPLAY   A   TABLE  OF   CONFIDENCE  INTERVALS   OF  ORDERED 
p  STATISTICS   WITH   CONFIDENCE   COEFFICIENTS. 

p  SUBPROGRAMS   CALLED   BY  THIS  FUNCTION  INCLUDE:    BINOM .    NORMCDF ,    NORMPTH , 
P  INPUT, AND   QUANC. 

♦iVu 
Nl:MENU  SIGNHELP 
a  MENUTWDICES  AND   ROUTE  FOR  PROPER  ACTIONS 


10]  N^:C-CHOICES  PACEDMENU  SIGNQBJ 
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*(C=2,3,u,5)/L8 

*(C=7,8.9,10)/i:9 

•*(C  =  11)/L20 
MENU  MAINOBJ 

fl"^        INPUT  DATA   FOR  SINGLE  SAMPLE  CASE 
L8:AAt-l 

X-^INPUT   1 

NNN'>-oX 

*(C=5)/L15 
MO-^INPUT  3 
D-^X-MO 

(T  PAIRED   SAMPLE  CASE 

R-rlNPUT   2 
Q*\t-R 

[3  2]  NNN-'-oDD 
[33]   -►(i:=10)/il6 
r3uT  W-f-INPUT   "• 
.35J   i/t-U-D-MO 

■3  6]  n  COMPRESS  D  TO  REMOVE  ZEROS 

373  Lll:Af(D:eO)/D 

p  RECORD   LENGTH   OF  A   AND   ASSIGN  TO   N 

N*pA 
a  KEEPING  TRACK  OF  POSITIVE  SIGNS 

KPOS**/{A>0) 
^(N>25)/NCRM 
PVALi^SINOM  N 
^<KPCS>Q)/P\ 
FVI~' 

PI :  PVI-*-!  -PVAL  ZKPOSl 
P2  :PV-rPVALL(.KPOS+l)l 

*L6 
fl  I?  N  IS   GREATER  THAN   3  0  USE  NORMAL  APPROX  V/   CONTINUITI   CORRECTION 
NORMlZl-  (  (.KPOS'O  .5)-(0.5x/V))^(C.5x  {N*0  .  5  )  ) 

Z-*-(  (aPOS  +  G.  5  )-(0.5xA')  )-i-(0.5x(/V*0.5)) 

PV^NORMCDF   Z 

PVIfl-iNORMCDF   Zl) 
p  IF  PAIRED   SAMPLE  TEST  GO  TO   LI 7  FOR  OUTPUT  STATEMENT 

I6:P73-«-2x(L/(Py,P7I)) 

-►(P73Sl)/iV5 

P73-<-l 
N5  :P7W-«-(3  ,1  )o  (.PV ,PVI .PVZ  ) 

'COMPUTATIONS  ARE  BASED   ON  A   SAMPLE  SIZE  OF:     ' .(fN),' . ' .QTCNL 

'THE  TOTAL  NUMBER  OF  POSITIVE  SIGNS  IS:     ' ,  i^KPOS) ,UTCNL 

-^(C—7    8  9)/L17 

'THE  P-i^ALUE  FOR   HO:    M   =    '  .(.vMO).'    VERSUS  HI:    M    '  .  (.fLOGICLiC-l)  :12  ) ,  '     ',( 

9M0),<    IS:     '  ,i^9PVM\:c-l;ll),DTCNt      

LIT: 'THE  P-VALUE  FOR   COMPARING   THE  MEDIAN  OF  THE  POPULATION  OF    ' 
' DIFFERENCES   TO   THE  HYPOTHESIZED   MEDIAN.     ' , UTCNL 
'HO:    M(.X-Y)    =    'AtiMQ),'    VERSUS  HI:    MiX-Y)    '  ,  (sLOCICC  (C-6  )  ;  1]  ) ,  •  '.(sMO),' 

,  IS:   •  (55Prac(c-5);i] )Tnrrro 

Lia: 'WOULD   YOU  LIKE  A   CONFIDENCE  INTERVAL  FOR  THE  MEDIAN?    il/N) .' ,UTCNL 

BB-^a 

^(BB='I' )/L16 

■^■QUANT 
R  INPUT  SIZE  OF  CONFIDENCE  INTERVAL 

L16:CC*INPUT   5 

ALPHAi-{100-CC)*200 

•>-iNNN>25')/NORMl 
p  COMPUTING  CONFIDENCE  INTERVALS   BY  EXACT  P-VALUE 

CDF-rBINOM  NNN 
p  INDEX  POSITION  OF   CDF  FOR  ALPHA    *    2 

Bir+ /  {CDF<.ALPHA) 

■*(.B>0)/SKIP 

fl-<-l 

•^SKIP 
P  COMPUTING   CONFIDENCE  INTERVALS   BY  NORMAL   APPROX. 
NORMl : KALPEA* f (0.5x( NNN* 0 . 5  n x ( NORMPTH   ALPHA ) )  +  C 0 . 5 x NNN ) - 0 . 5 
a  ROUND    -lALPHA    DOWN   10    '.lEAREST   INTEGER   AND    INCREMENT   3Y   ONE 

H-^^.KALPHA  +  l 
a  r?  SINGLE   SAMPLE   CASE   GO   TO   LI 

p   *   "  '  •  'calculate  and  print  out  conf.  int.  for  paired  sample  case 

L5:ORDD*DDZiDD2 

YY^oORDD 

CI-^-ORDDLBl  ,0RDD^(.YY-(.B-1))1 

'  A    '.(aCC),'    CONFIDENCE  INTERVAL   FOR  THE  MEDIAN  OF  THE    ' 

< POPULATION  OF   DIFFERENCES   IS:' , DTCNL 

'  (  '.(SCJCI]),'  ^   MEDIAN  U-Y)    S  ', (.901121) ,'     )'.QTCNL 

*QUANT 
P  CALCULATE  AND  PRINT  OUT   CONF.    INT.    FOR   ONE  SAMPLE  CASE 

L7:0RDX*XliXl 

YI*dORDX 
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110] 

lllJ 


ClOOl  CT-^ORDXZB'l  .ORDXKYI- (.3-1))'^ 

tlOl]  'A    '(.9CC).'       CONFIDENCE  INTERVAL   FOR   TBS  MEDIAN  OF  THE  POPULATION  IS:     < 

,  DTCNL 
[102]   '  (  '.(5CJC1]),'  i   MEDIAN   £    '(^CILll)'     )'.OTCNL       ^ 

L103J  QUANT: 'WOULD   lOU  LIKE   CONFIDENCE  INTERVALS   FOR  A   SPECIFIED   QUANTILE?    (.I/N 

) . ' .QTCNL 

■»"(VW='J' )/fll 

L2C:< ENTER   TEE  SIZE  OF  TEE  SAMPLE.' 

,  NNN->-C 

Cl09j  31: 'ENTER  DESIRED   QUANTILE;    FOR  EXAMPLE:    ENTER   20,  FOR  THE   20TH  QUANTILE. 

QUA-i-D 

->-(.(QUA^O  )v  (qi;a>ioo  )  )/El 
112]  NNN   QUANC   QUA 
113]   •  ' 

im]   I  *****  THIS   TABLE   GIVES   CONFIDENCE   COEFFICIENTS   FOR   VARIOUS   INTERVALS    ' 
115]  'WITH   ORDER   STATISTICS   AS   END  POINTS   FOR   THE    >  A^QUA) ,'TH   QUANTILE .' .QTC 
NL 
[lis:  'PRESS   ENTER   WHEN  READI .  ' 
[1171  33-^^ 
[118]   -A?u 

,113:  El:  -ERROR:    THE   QUANTILE   VALUE  MUST  LIE  BETWEEN   0  AND    100;  TRI  AGAIN.' 
Cl203  -31 
7 


7  S'^MAN '  X-"^  •  A-0'  ^'  '~HA  •  a  3  •  =  •  '^y 
a    *"  this' FUNCTION    :CMPUTES   THE   SPEARMAN  R   STATISTIC   VHICH   MEASURES 
a  THE   DEGREE   OF   CORRESPONDENCE   3ETVEEN   RANKINGS    OF   IWO    SAMPLES.    THE   ?- 
^VALJE   IS    JI7EN   FOR   TESTING    JNE   AND    r,^0- SIDED   JIF0THESI3    JF  ASSOCIATION 
0  SUBPROGRAMS   CALLED   31   THIS   FUNCTION   INCLUDE:    TIES,    IIESK.    SPEARP , 
R  INPUT,    SPAPROX,    INTERP,    AND   THE   VARIABLE  PMATSP. 

0         DISPLAY  MENU  AND   INPUT   DATA. 
Nl:3*MENU   SPMANOEJ 

-^(5  =  1)/S1      

MENU   MAIN03J 
*0 
Bl:R*INPUT  2 
Q-^l+i? 

X*1*(.Q+1)^R 
J*(fl+l)+i? 
0  CALL   SPMANP  TO   CALCULATE  TEE  STATISTIC  AND  ASSOCIATED  P-VALUES 
A*X   SPMANP   I 
-^(ACi]>o)/Li 
CHA-ri  INDIRECT^ 
■*L2 
Ll:CHA-^' DIRECT' 
22]  L2:P7-«-2xAC2] 
*"*      *{PV^1)/L2 
,.-  -  PV*1 

;25l  L3:' SPEARMAN  "S  R  EQUALS:    '  .(nvAUll)  .UTCNL 
26    'THE  P-VALUE  FOR  HO:    NO   ASSOCIATION   EXISTS   VERSUS' 
^""^   '  Hi:     '.(iCHA),'    ASSOCIATION-EXISTS   IS:     '  .(.nvAZll)  .DTCNL 

'THE  P-VALUE  FOR  THE  TWO-SIDED  TEST  OF  HIPOTHESIS  IS:    ',  N« (2x^12]  ))  ,DrCi7 

'PRESS  ENTER  WHEN  READI.' 

^Nl 

7 


27 
[28 


7  WISIG;A:B:D'.E;F:PV2'.Zl:Z;DEN;NUMZiNUMZl:PVALiX:M0',N',TPLUS:CDF:TALPHA',AL 
PBAiBiCI;IiAA;BB:CCiNN:DD;PViP0S;TP0SiNM;PVI;TP0SliNNN:C;PVMiPV3;R;Q;NU 

M ' Tc ' TCl ' ThAP ' DENT • DENTl 
1]   ft  THis   tUNCilON   OsES   hs  WILCOXON  SIGNED   RANK  TEST  TO   CALCULATE  THE  TPLUS 
"^   ft  STATISTIC.    P-VALUE,    AND   CONFIDENCE  INTERVAL   AS  A   TEST  FOR   MEDIANS. 
ft  SUBPROGRAMS   CALLED   31  THIS   FUNCTION   INCLUDE:    TIES,    WILP ,    NORMCDF , 
P  N0R11PTH,    CONFW,    INPUT  AND   THE   VARIABLE   PMATRIX. 
•;  5  J    '.'Vu 

•,a:       Nl:MENU   VIL3EZ? 

.7]       ft  ~FiE7}U   CHOICES  AND   ROUTE  FOR  PROPER  ACTIONS. 

8]      NuiC-f-CHOICEW  PAGEDMENU  WILOBJ 
9]         *(C=2,3,u,5)/L8 


2 

:3] 


♦(C=6.7,8,9)/i9 
->-(C=l)/A'l 
MENU  MAINOBJ 

ft  INPUT  DATA   FOR  SINGLE  SAMPLE   CASE 

L8:AA*1 

X-<rINPUT   1 

NNN-^oX 

*(C=5)/Z;i6 

HOi-INPUI   3 
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D*X-MO 

0  PAIRED   SAMPLE   CASE 

L9:AAf-2 

R-^INPUT   2 

Q-f-lfR 

DD->-X-J 

NNN-f-oDD 

*(C=9)/Z:i6 

m-^INPUT   4 

D-^iX-D-MO 
p  COMPRESS  D  TO  REMOVE  ZEROS 

Lll:A'^(D^O)/D 
R  RECORD   LENGTH   OP  A   AND   ASSIGN  TO   N 

N-^-qA 
fl  KEEPING  TRACK  OF  POSITIVE  SIGNS 

POS-^(.A>0) 
q  TAKE  TEE  ABSOLUTE   "ALUS   OF  A;    ASSIGN   TO   3   AND   ORDER   3 

8'>-\A 

a       REORDER   POSITIVE   SIGNS   TO   COINCIDE   WITH   PROPER   POSITIONS   IN   3 

POS-^-POSli^^.  \Ajj 
n  CALL  FUNCTION  TO  BREAK  TIES 

.   ,  E-^l   TIES  B 

46]  p      CALCULATE  TPOS  BI  ADDING  ACROSS  ALL  POSITIVE  VALUES  OF  E 
•■■-■*  TPOS->-+/(.POSxE) 
TPOSl-^TPOS 
fl         GIVES  SIZE  OF  LEFT  TAIL   OF  PROBABILITY  DISTRIBUTION 

iVM-«-(L((+/xiV)-^2)  )+l 
p  GO   TO   STATEMENTS   3ASED   ON   LENGTH   OF   VECTOR   E 

'32]  *(N>'5)/L3 

^53]  P  GENERATE  NULL   DISTRIBUTION  FOR   TPLUS 

[54]  F^WILP  N 

[5  5]  p  IF  TPOS  FALLS  IN  LEFT  HALF  OF  PROB   DIST  CALCULATE  PVALUE  AS  NORMAL 
[55]  Ll:^iTPOS<,(.NM-\))/TP.EG 

[57]  fl  OTHERWISE   USE   THE  NEGATIVE   T   STATISTIC 

[581  TPOS-^(.*/\N)-TPOS 

[59]  a  IF   TPOS   IS   FRACTIONAL    USE   BOTH   THE  INTEGER  ABOVE  AND   BELOW  AS   TPLUS 
l50j  TREG:^{  :i  '  TPOS'.  =  0  )/N0N 
61.  TPLUS*\.TPOS 

.62  P\:PV-*-l-{{FLTPLUS'\+F^TPLUS->rl'\  )*2) 
53    P7J-^-(FCrPLi/S+l]+FCrPLUS  +  2]  )  +  2 
54.  *CHECK 
'6  5  N0N'.-*{TP0S>Q)/GO 
55    PVf-l 
57.   -►P2 

58  GO'.PV^l-FL^TPOS)! 
.69   P2:P7J*FC(2'P05  +  1)] 

70  CHECK:*{TPOS'l<,(.NM-l))/L& 

71  Pl'2*P7 

72  PV-^PVI 

73  P7J-<-P72 
7H  ♦L6 

7  5   R  COMPUTE  NORMAL   APPROX.    W/CONTINUITI   CORRECTION  FACTOR 

76  L2:TRAP'^{N^{N+X))*<* 

77  NUMZ-<r{TP0S-^Q.5)-TRAP 

78  NUMZ1*-{TPOS-O.i)-TRAP 

^79.  riEA'*((A'x(iV+l)x((2xiV)+l))*24)*0.5 

80  Zi-NUMZ*DEN 

.8i:  Zl-^NUMZlrDEN 

82]  p  COMPUTE   STUDENT  T  APPROXIMATION  WITH   CONTINUITY   CORRECTION  FACTOR 

.83]  NUM*\(,TPOS-TRAP) 

8U  DENT-*-l{  (N^  (.DEN*2  ) )+  (A?-l ) )-  ( ( (NUM-O  .  5  )*2  )+  (JV-1 )  J  )*0  .  5 

8  5  P£'iVri-t-(((yVx(C£'W*2))r(iV-l))-(((yvra+0.5)*2)+(iV-l)))*0.5 
86  TCir{NUM-0.5)rDENT 

"87  rCl-<-(iVUM  +  0.5)*£i£:W2'l 

88  p  COMPUTE  AVERAGE  OF  TC  AND   ZC 

.89.  *(.TP0S<,(,(+/\N)*2))/SEC0ND 

90  PV-t-iCl-iNORMCDF   Zl) )  +  (.!-( (N-1)   TDISTN  TC1)))*2 

91  PVI*l(.NORMCDF   Z)  +  ((iV-l)    TDISTN  rC))  +  2 
1.92,  -►z;6 

.93.  SECOND  :PV*-(  (1-  (NQRMCDF   Z1))  +  ((N-1)    TDISTN  TCl))*2 
quj       PVIi-(.(NORMCDF   Z)  +  (l-((N-l)    TDISTN   TO))*? 

',9  5]  L6:P73*2xa/  (PV.pyj}  ) 
;96j       -|^?73S1)/;V5 

;98J  iV5:P7M-»-(3,l)p(P7J,P7,P73) 

,99]  'COMPUTATIONS   ARE   BASED   ON  A   SAMPLE   SIZE  OF    '     (vN)    ' .  <OTCNL 

.100]  'THE  TOTAL   SUM   OF  POSITIVE  RANKS   IS:     '(.vTPOSD.'.'.OTCNL 

.101.  R  ^         IF  PAIRED   SAMPLE  TEST   GO   TO   111    FOR   OUTPUT   STATEMENT 

,102.  ->-(AA  =  2)/Ll7 

.103]       'THE  P-VALUE  FOR   HO:    M   =    '(9M0),<    VERSUS   HI:    M    '  ,  (.fLOGIClC-l;ll  )  .  '     '  ,  (» 

^  MO).'  IS:'  .(^9PVMLC-l;l2).ilTCNL 

:i0U]  ♦LIS 

.105.  L17:'      THE  P-VALUE   FOR   COMPARING   THE   MEDIAN  OF  THE  POPULATION  OF    ' 
,106.   'DIFFERENCES   TO   THE   HYPOTHESIZED   MEDIAN.     ' .DTCNL 

Cl07]   '50:  M(.X-Y)    =  'AiMQ).'    VER^U^Hl:    MiX-Y)    '  ,  ('JiOC:JCCC-5 ;  1]  ) ,  '  '.(wMO),*, 
IS:'  .(^9PVMLC-i:l2).CiTCNZ 
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L18:' WOULD   lOU  LIKE  A   CONFIDENCE  INTERVAL   FOR  TEE  MEDIAN?    (.I/N).' 
*(SS='I' )/il5 

LlS-.CC-^INPUT    5 

ALPaA-^-UOO-COilOO 
fl         ROUTE  TO  NORMAL  APPROX.    FOE   CONF  INT  OF  LARGE  SAMPLE  SIZE 

■*(NNN>9)/Ln 

CDF-*rWILP  NNN 
a  INDEX  POSITION   OF   CDF  FOR  ALPEA    *    2 

CDF'^(.CDF^O)/CDF 

T  ALPS  A*  ( -r/  (CDF^ALPEA  ) ) 

*(.TALPSA::OWJUMP 

TALPHA-^1 

*JUMP 
a    COMPUTING   CONFIDENCE  INTERVALS   BI  NORMAL   APPROX.    W/C.F. 
Ln:TRAP'>-(NNNx(,NNN+l))*n 

DENt-iiNNNii  (.NNN+1  )x  (  (IxNNN'}*!  )  )*2U-  )*0  .  5 

T ALPEA* {DEN* (NORMPTH   ALPHA ) )+TRAP-0 . 5 

TALPHA 
a  ROUND   TALPHA   DOWN  TO   INTEGER   VALUE  AND   INCREMENT  31   ONE 

TALPHAf-lTALPHA  +  1 
a  IF   ONE   SAMPLE   CASE   GO   TO    LI 

JUMP:^{C=2  ,2  .'4  ,1) / LI 

fl  CALCULATE  AND  PRINT  OUT  CONF.    INT.    FOR  PAIRED   SAMPLE   CASE 

L5:CI*TALPBA   CONFW  DD 

'  A    >  .{900,'      CONFIDENCE  INTERVAL   FOR  THE  MEDIAN  OF  TEE    ' 

' POPULATION  OF  DIFFERENCES   IS:' , DTCNL 

•  (  '.(sCJCl]),'  i   MEDIAN (X-7)    S  ',(»CJC2]).'  )<  .UTCNL 
'PRESS   ENTER   WHEN  READI . ' 

3B->-a 
■*NU. 
a  ::ALCULATE  and   print  out   CONF.    int.    FOR   ONE   SAMPLE   CASE 

l7 -.01*1  alpha  confw   x 
'a  '  .i9cc),'     confidence  interval  for  the  median  of  tee  population  is:',. 
utcnI 

'  (    '  .(iJsCJCi:  ),'  £  MEDIAN  ^    ',(5CIC2]),'  V  .UTCNL 

•  PRESS  ENTER  WEEN  READ! . ' 

BB-^a 

■*NU 
7 
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APPENDIX  E 
MAIN  PROGRAM  LISTINGS  FOR  MAINFRAME  COMPUTER  WORKSPACE 


7  KENDALL: A: AAiB;BX;BI'.C:CXiC7;DiDD;DX:DIiDXIiSiP0S;NEG;XX;II:N '.DEN ;NN;NU 

M:PiPV:PVAL:SU'.SV',T:UiV;AT:Z:X:Y:CHA:Q;R 
1]   fl    is is   JUNCTION   COMPUTES   THE  KE^dAlL   B   STATISTIC  WHICH   IS   A   MEASURE 
2]   n  OF  ASSOCIATION   BETWEEN   SAMPLES.    P-VALUES   ARE   GIVEN   FOR   TESTING   ONE 
3I  a    AND    TWO -SIDED    HIPOTHSSIS    ?0R   NO   ASSOCIATION    VERSUS   ASSOCIATION. 
■+]   fl  SUBPROGRAMS   CALLED   3i    THIS   FUNCTION  INCLUDE:    TIES,    TIESK,    KEN D ALP , 
5]   n  INTERP   INPUT  AND   NORMCDF . 
5]    R--INPUI   2 
7J    Q-<rl^R 
8J    X-«-l+(Q+l)+iJ 
93    y*((3+l)+5 

10]  p  ORDER  I  IN  INCREASING  ORDER  OF  X 

AfJLiXl 

p  ORDER   X   IN  INCREASING  ORDER 

B^XLi^l 

fl  COMPUTE   CUFJiENT  RANKING   OF  I 

p  NOW  ORDER   I   RANKS   IN   INCREASING  ORDER 

D'>-AZU1 
fl      IF  TIES  EXIST  IN  EITHER  X  OR   I  RANKED   VECTOR   USE  MID -RANK  METHOD 

DD--1   TIES   D 

XX-^l   TIES   3 
R         FIND   ORIGINAL   RANKING  OF  I  WITH   TIES  RESOLVED 

Il^DDiCZ 

N*oX 
p  COMPUTE  NUMBER   OF  DISTINGUISHABLE  PAIRS 

NN-^iN^iN-l))*! 

S-t-pO 

AAi-0 
p  POSITIVE  ONES   COME  FROM  A   RUNS   UP   CONDITION:    NEGATIVE   1  FROM  RUNS   DOWN 
P  ZERO   IS   SCORED   FOR  TIES.    MULTIPLY  THE  RESULTS  FOR  EVERY  ELEMENT  AND   SUM 
Ll:AA-^AA+l 


.11 
12 

13 

.1^ 

.5 

.6 


18 
.19 
.20 
.21 
.22 
23 
r2U> 
"25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
=  9 
60 
51 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 


BX'^(XXl,AAl><,AA*XX}) 

CXirlxxUA2<UA*XX)Ui    1) 

DX-^BX+CX 

B7->- (,YICAA1>(,AA -¥77 )) 

CY-*-lYYLAAJi<lAA-^YY)U(~l') 

DYf-BY+CY 

DXY-rDX^DY 

P0S'^(,DX7>Q} 

NEG-*-lDX7<oUn) 

S-^S.POS.NEG 

■*UA<(.N-1))/L\ 
fl  SUM  FINAL   VECTOR  TO   DETERMINE  S 

S*+/S 
n     OBTAIN  THE  NUMBER  OF  TIES  IN  EACH   VECTOR   USING  THE  TIESK  FUNCTION 

U-^TIESK   B 

V-i-TIESK   D 

SUlr+/(2'.U) 

SVf-+/(2'.V) 
p  CALCULATE  THE  B   STATISTIC  INCLUDING  THE  CORRECTION  FOR  TIES 

Ti-S-^((NN-SU)x(.NN-SV))*0.5 

AT-t-lT 

■*{N>12)/N0RM 
n      CALL  KENDALP  TO   CALCULATE  THE  RIGHT  TAIL  OF  THE  CDF  OF  B 

P*KENDALP  N 
R  CALL   INTERP   TO   CALCULATE  P -VALUE  B7  INTERPOLATION 

PVAL-AT  INTERP  P 

•^(PVAL^~1)/L3 

PVAL*0 . 5 

0  CALCULATE   P   7ALUE   USING   NORMAL   APPROX . 
N0P.M:.-'IUM^\2*AT)xi  :2^NN  j*0  .  5  j 

D£'A'-<-(2x((2x/V)  +  5))*0.5 

Z*NUM*DEN 

PVAL-^l- (NORMCDF   Z) 
ft    IF   B   IS   POSITIVE  PRINT  OUT  DIRECT  ASSOCIATION. 
L3:*(T>0)/L5 

CHA*' INDIRECT' 

*L7 
L5:CHAi-<  DIRECT' 
L7:PV-<-2i(PVAL 

*(PV^1)/L8 

PV*1 
L8:<  KENDALL"  S  B   EQUALS    •  ,  (4ur) 
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75 
76 
77 
78 
79 


'TBE  P-VALOE  FOR   BC: 

'  Hi: 


NO   ASSOCIATION   EXISTS    VERSUS' 

•  ,  (.9CHA  ) ,  'ASSOCIATION  EXISTS  IS:     '  ,  i^vPVAL  ) 


'THE  P-VALUE  FOR  TEE  TWO-SIDED  TEST  OF  3IP0THESIS   IS:  '  ,  (U!isP7) 
I  I 


f|3 

H 

u 

5 

15 

15 

17 

18 

19 

20 

21 

22 

23 

24] 

25] 

26] 

28] 
.29^ 
.3C 

31 

32 

33 

34 

35 

36 

37 
,38 

39 
>0 

>2 
43 

4U 

>7^ 
48 
49 

50 

.51 
52 

53 

54 

55 

56 

57 

58 
.59 
,60 
.61 
.62 
,63 
L54 


•20 

67 

68 
.69,, 
.70 

71 

72 

73 

74 

75 

76 

77 

78] 


7  KE0SKAL:N0M:DEN0yi'.A:C',E;D:K',AA-,3B;DD;E:FiN;0F',P;PVAL',RiS0FR',SR',TS0E;CBA 

fl     2*515  FUNCTION   COMPUTES   THE   KUSKAL-VALLIs   TEST   STATISTIC   3   WEICS    IS 

fl  A    MEASURE   OF  TBE   EQUALITY   OF   K   INDEPENDENT   SAMPLES. 

B  SUBPR0ORAM3   CALLED   EI   TBIS   FUNCTION   INCLUDE:    TIES.    TIESK,    INDEXPLS . 

a   FDISTN,    INTERP  AND  TBE   VARIABLES  PMATKW20 ,    PMATKWil ,    PMATKW3  3 , 

0  PMATXWi^,    PMATKWni,    PMATKWni ,    AND   PMATKM^i . 

UPP-^5 
Bl:' ENTER  TBE  NUMBER   OF  POPULATIONS  TO   BE  COMPARED    (.MUST  BE  GREATER  THAN  T 

WO).' 

-(Ti?<3)v((l|Z)=:0))/?l 

*IU&)>1)/E1 
a  INITIALIZE   VECTORS   E  AND   F  AND   VARIABLE   C 

E-'-F'^SCFR-^o  0 

rTEIS  loop  FACILITATES  ENTERING  TBE  SAMPLE   VECTORS  AND   STORING  TBEM 

CBA-*- '  FIRST ' 
Ll:C:*C+l 

' ENTER   JOUR    ' , («CHA ) , '  SAMPLE. ' 

D-^Q 

D-lcD 
a         CONCATENATE  SAMPLES  AS  TBEI  ARE  ENTERED   AND   STORE  TBEM  IN   VECTOR   E 
NEXT:I^E.J 
a  RECORD  TEE  LENGT3S   OF  TBE   SAMPLES  AS  TBEI  ARE  ENTERED 

F-'tF.qD 

CBA*'NEXT' 

*(.C<(.K-1))/L1 

C3A->-'LAST' 

*(,C<K)/L1 
9  RECORD   SIZE  OF  ALL   SAMPLES   WEEN   COMBINED 

N'>-'/F 
a    ORDER   SAMPLE  SIZES  LARGEST  TO   SMALLEST 

fl  ORDER   COMBINED   SAMPLE   VECTOR  TO   BE   USED   BI  TIES   FUNCTION 

D*ELtE1 
fl    CALL   INDEXPLS  TO   INCREMENT  INDEXES  WEEN  TIES  OCCUR  WITBIN  ONE  SAMPLE 

AA-rF  INDEXPLS   E 
p  CALL  TIES  TO   BREAK  TIES  BI  MIDRANK  METEOD 

BB*1   TIES  D 

C-*-0 
B  TBIS  LOOP  CALCULATES  TEE  B  STATISTIC 

L2:C-*-C+l 
B  SUM  OF  RANKS  FOR   EACB   SAMPLE  IS   CALCULATED 

Si?-«-+/SSC(FCC]  +  (MCC;]  ))] 
fl  CALCULATE   SUM  OF  RANKS  SQUARED   DIVIDED   BI  TBE  INDIVIDUAL   SAMPLE  SIZE 

SR*(.SR*2)-tFLC:i 
B  STORE  EACB   CALCULATION 

SOFR->-SOFR,SR 

*{C<K)/L2 
fl       ,      SUM  ACROSS  ALL  SAMPLES 

TSOR-*-*/SOFR 

B-<- (TSORy:  (12*  (.Nx(.N+l))))- (3^  (N+D) 
fl  RECALCULATE  B   WITB   CORRECTION  FOR  TIES 
Af-TIESK  E 

NUM'<-(+/(Ai'3))-(  +  /A) 
DEN0M-<rN*((N*2)-ll 
H-<rBi(l-(NUM*DENOM)) 
fl  SYSTEM  OF  LOGICAL   STATEMENTS  ENSURE  PROPER  PROS.    IS  ACCESSED 
[OF lll<2) /OUTPUT 
\K>n-)/FAPPROX 
*{k=3)/IF 

(A/_(OF-    3    3    3    3))w(oni'2>3))/FAPPROX 

2    1    1)}v(a/(0F=    3    3    1    l))v(A/(Oi?=    3    3    2    1)) 


*(A/(OF=    2    11    1})/0UTPUT 
♦  f  fK=^  }a  (OF''''  ■'='>;  ^/DiiT 
*('a/(0='=    3'l'i    1;)v(a/(CJF=    3 


)/DU-t 

*CZ=i)/?u3 
IF:*(0F[.12>^)/FAPPR0X 
■*-((^/(0F=    2    1    1))v(a/(0F=    3    1    1))) /OUTPUT 
-►((OF[1]=4)a(OFC3  3=1))/P31 

l}r^j)°nl^    I    I    l\V>'Xi>°A^.\^'^,T,V'^''^^^^°^^    3    3    2))v(A/(0F=    3    3    3)))/P33 
-►((a/(OF=    3    2    2  ))v(0Fl1]=4)  )/P34 

fl  CALL  APPROPRIATE   VARIABLE  ACCESS   CDF 

P23:Pi-PMATKW20L(N-ii);;l 

■*-PM 
P31:Pf-PMATKW3lL(N-5)i;l 

■*-PM 
P33:P*PMATKW33L(.N-5);  ;] 

•>-PM 


'.  •  •'■ 
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P3  U  :  P->.PMATKWZ  U  C  (.V-  6  )  ;  ;  ] 
-►PM 

■*-PM 
P^7:P*PMATKW^2L{N-z); ;] 

*PW 
P43  . p^PMATKWn 3C(N-7)',:1 

a  CALL   JNTERP  TO   CALCULATE  P-VALOE  37  INTERPOLATION 

PM:PVAL*U  INTERP  P 

•^(.PVAL=    1)/ OUTPUT 

■^L5 
a      CALCULATE  P -VALUE   USING  TEE   F   DIST  W/ONE  LESS  D.F.IN  DENOM  APPROX. 
FAPPROX:F^{  <,N-K)*-H)*  (.  U-1}^  (  (.N-D-H  )) 

PVAL^l-iilK-D.iiN-K')-!)^   FDISTN  F) 

■*-L5 
OUTPUT :  PVALi- '  GREATER   TEAN    .25' 
L5:'THE  H   STATISTIC  EQUALS:     • , (4«ff ) 


THE  P-ViALUE  FOR   BO:    THE  POPULATION  MEDIANS  ARE  EQUAL   VERSUS    • 


L  3  3  J  >J  J 

Cioo]   '  ' 

[lOlJ   -^0 


:iC2J  El: 'ERROR:    IGU  MUST  ENTER  A   SINGLE  INTEGER   VALUE   GREATER   IRAN   2;    TRI   AGAI 
N.  " 


[103]   •  • 
[104]   *51 
V 


1 
2 

3 
U 

5^ 
5 

[13 

.10 

11 

12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
C22 

;23 

,24 

25 

.26 

127 

.28 
.29 
!30 

:3i 

.32 

C34] 


7  MANNWHIT:N:M'.PV2-,A:3:C:G:MM'NN'.RX:U:NM1;?-NU:PVAL',NM:NUMZ:NUMZ1:DENC1:Z 
;  Z 1 ;  ALPHA  -.CDF;  INDEa  :  IPX :  CI ;  u  ALPHA  ;  33 -.CC ;  Jl:U2  :  PV  :NN1:  NN2  :  PVI ;  DIFF  -AAiGG 
iPVM :PV3  : D :  Q -, R  •. DEN -.DENG iTC  ■, TCI ; NUM 
p     THis   fOnctIoN   USks   the  sum   OF  RANKS  PROCEDURE  TO   CALCULATE   TEE 
fl  MANN-WEITNEY   U   STATISTIC   USED   IN   COMPUTING   THE   P-VALUE   FOR   THE  TEST 
B  OF   LOCATION   AND    SCALE.       THE   C.I.     FOR   M{I)-M{X),    THE   SHIFT   IN 
p  LOCATION,    IS   ALSO    COMPUTED.    SU3PR0GRAMS   CALLED   INCLUDE:    TIES,    TIES2 
fl  INDEXPLS,    VARMW,    MANWP ,    INPUT.    CONFMW ,    NORMCDF ,    AND  NORMPTH . 
DIFF-fO 
I  I 

'DO   700  VISE  TO   COMPARE  THE  MEDIANS   OR   VARIANCES   OF  TEE  POPULATIONS?' 
t  t 

32:'         ENTER:  1   TO   COMPARE  MEDIANS;  2  TO   COMPARE   VARIANCES.' 

AA-t-a 

*{UA^1)A(AA*2))/E2 

■>'Iaa=i)/ni 

'      the  test  to  compare  variances  requires  the  two  population  means' 
'or  medians  to  3e  equal.     if    the7  differ    37  a  known  amount.   ' 
'the  data  can  be  adjusted  before  appl7ing  tee  test.' 

t  I 

53:'     ENTER  TEE  DIFFERENCE  OF  MEDIANS    (M(X)  -  M(I))  OR   900  TO  QUIT).^ 

DIFF-^0 


■((.pDIFF)>l)/E3 

■Idiff=';     ■  ■ 


:900)/0 

'THE  NULL  B7P0TBESIS  STATES   -   TBE  POPULATION  VARIANCES  ARE  EQUAL;    7(X)  = 

va).' 

'  WBICH  ALTERNATIVE  DO  70U  WISE  TO  TEST?' 

I  I 

N2:'ENTER:      1  FOR   HI:  V{X)   <    7(7);   2  FOR  Bl:    7(X)  >  7(7);   3  FOR  Bl:    7(X) 
*  7(7). • 

-►(  (£1*1  )a  (P*2  )a(D*3  )  )/ffU 
-31 
Nl: '      TBE  NULL  B7P0TBESIS  STATES   -  TBE  MEDIANS  OF  X  AND  7  ARE  EQUAL;    MiX) 
=   M(7). ' 
I  I 

•  NHICE  ALTERNATIVE  DO   70U  WISE   TO   TEST?' 

I  I 

B6:'ENTER:      1  FOR  El:    M(X)  <  «(7):   2  FOR  Bl:    MiX)    >   M(7);   3  FOR  Bl:   M(X) 
=  M(7).  ' 

[36]   -►T(Ii3:l)A(D22)A(ZJ*3))/ff4 

[37]  fl  FNTER   DATA    VECTORS 

v38l  31: 'ENTER   X   DATA    (MORE  IB AN   ONE   OBSERVATION   IS   REQUIRED).' 
.3  9  J   ;V*a 
40j   *(.  (op/V)  =  0)/El 

41.  'ENTER   7  DATA. ' 

42.  M*Q 

43.  a    IF   CALCULATIONS  INVOLVE   VARIANCES  ADJUST  X   37  TBE  DIFFERENCE  IN  MEANS 
4  4    N*N-DIFF 

4  5.  p  CONCATENATE  X  AND   7   SAMPLE   VECTORS 

46    A-*-N,M 

47.  p  DETERMINE  SIZE  OF  X  AND   7   VECTORS  AND  ASSIGN  TO   NN  AND  MM 

^.48.  NN-^pN 
.49.  MMi-pM 

50.  p  COMPUTE  SIZE  LIMIT  OF  LEFT  TAIL   OF  NULL   DISTRIBUTION 

51.  NMi-(NNxMM)*2 
C52J  NM1*INM 
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.53 

5" 
.55 

,^59 
60 

.61 
62 

!63 

.66 
67 
68 
69 
70 
71 

>72 

[74 

[75 

^75 

^77 

78 

79 

80 

1^81 

'83 
[su 
J85 

l87 

88 

89 
.90 
,31 

;li 

.98 

99_ 

100 

101 

;i02 

103 

,10U 

105 

.106 

107 

108 

.109 

.110 

.111 

Cll2 

;il3 

.IIU 

.115 

.115 

,117 

.118 

119 

120 

,121 

,122 

:i23 


;i2U 
125 

126] 
127 


fl  ORDER  A   AND  ASSIGN  TO   B 

B*-AiU.l 

C-<-(.NN.MM)   INDEXPLS  A 
fl       CALL   TIES  FUNCTION  TO  BREAK  TIES  USING  MIDRANK  METHOD 

G-1   TIES   3 
a    IF  FALSE   CALCULATE  TEST  FOR   VARIANCES 

*(Ail  =  l)/S5 
fl      CALL   7ARMW  TO   GENERATE  RANKS  REQUIRED   FOR   VARIANCE  TEST 

GGf-VARMWCNN+MM) 
a    CALL   TIES   TO  RECORD  TIES  IN  THE  DATA   AND   BREAK  TIES  IN  GO 

q''   "  "  ,  CALCULATE  SUM  OF  X  RANKS 

S5:.W^+/(CC(CC1;  (i/ViV)]  )]  ) 

0      , ,  CONVERT  TO  MANNWHIT  U  STATISTIC 

U-^-RX-iiNNyiNN+l))*!) 

Ul-^U 
n  ,,„   IF  SIZE  OF  X  TIMES  SIZE  OF  I   >   80;    GO  TO  NORMAL  APPROX 

*((WMx2)>80)/L2 

NNl-<-\/NN.MM 

NN2'-'  /NN.MM 
QMAmP   FUNCTION  CALCULATES  LEFT  TAIL   CUMULATIVE  PROBS .    OF   U  STATISTIC 

P->-NN2    MANWP   NNl 
fl  LOGICAL   STATEMENT  ENSURES   ONLI  LEFT  SIDE  OF  NULL   DIST  IS   USED 

*(.U-S.NM1)/L2 
fl  CONVERT  U  STAT  WHEN  GREATER  THAN  LEFT  TAIL   VALUES 

U-*-(NN^MM)-U 
p  IFUIISA   FRACTIONAL,    INTERPOLATE  P   VALUE 

L3:->-(.(.l\U)  =  0)/NON 

U2*IU 

*(U2>0)/P1 

?V-t-l-(Cpr:;2+i]  )42) 

-?3 


P1:PV'<-1-  ( (.P':U22+P[,U2+11  )*2) 
i»3:?vr-^(?[!72  +  i:+PLi72  +  2J       ' 


_    +2 

■*CEECK 
N0N:*(.U>QWGO 

PV->-l 

♦  P2 
GO':PV->-l-?rui 
?2:P7I-^?C(i/  +  l)1 
CHECK :  *  ( Ul^NMl  ViU 

PV-^PVI 
PVI*PV2 

p  COMPUTE  THE  ■  NORMAL  APPROXIMATION  W/ CORRECTION  FACTOR 
L2:NUMZ*(,U-i-Q.5)-NM 

NUMZl-t-lU-0.5)-NM 

DEN->-(.(MMxNN*(,MM+NN+l))*12)i'0.5 

Zi-NUMZ*DEN 

Z1*NUMZ1*DEN 

NUM^\ (U-NM) 

DENCA  (  (  (NN+MM-1  )x  (DENi'2  )  )*  (NN+MM-2  ) )-  (  (  (NUM-0  .  5  )*2  )*(NN+MM-2  )  )  )*0  .  5 

DENCl*i  ( ( (.NN+MM-1  )x  (,DEN*2  )  )♦  (iVA'+WM-2  3  )-  (  (  iNUM+0  .  5  )*2  )*iNN+MM-2  )  )  )*0  .  5 

TC-*-(.NUM-Q.5)-tDENC 

TC1*{NUM+0.5)*DENC1 

-(U^NM)/ SECOND 

PVI-^UNORMCDF   Z)  +  ((.NN+MM-2)   TDISTN  TC))*2 

PV*Ul-(.NORMCDF   Zl))  +  (l-((iVW+MM-2)  TDISTN  TC1)))*2 

*LU 
SECOND :PVIf(.(NORMCDF   Z)  +  (.l-((.NN+MM-2)   TDISTN  TC)))*2 

PVf((.l-(NORMCDF   Zl})+(.(.NN*MM-2)   TDISTN   rcl))  +  2 

^(.PV3^1)/N5 
P73+1 
N5 :P7M+(3 ,1 )p (PVI.PV.PV3  ) 
i^THE  SUM  OF  THE  X  RANKS  IS:    ',(«5X),'.  THE  U  STATISTIC  EQUALS'.    \  (.9U1) 

P         LOGICAL  STATEMENT  FOR   VARIANCE  OUTPUT 
^(AA=2)/VAR 

'    THE  P-VALUE  FOR  HO:    M(X)  =  M(7)  VERSUS  HI:    MiX)    '  ,  (vLOGICZDiU  )  ,^    M(.I) 
IS:    '  .C+jfpyMCZ?:!]  ) 
I  I 

■*-L8 
VAR:PVM*(3  .Dp  (PV  ,PVI  .PV3) 
'THE   P-VALUE   FOR   HO:    7 {X)    =  7(1)  VERSUS   31:    7(X)  ',  (.'sLOGICZD ;  12  ).  '    7(7) 
IS:     '  ,  (.'^^PVMZD-.lZ  ) 


',1281 

l129]  *0 

Cl3  0]  L8:' WOULD   lOU  LIKE  A   CONFIDENCE  INTERVAL   FOR  THE  SHIFT  IN  LOCATION {Mil) 
M(X))?      a/N).' 

131]  flS-^Q 

132.  -t-iBB^'N'  )/0 

^133.  LlO:CC-<rINPUT   5 

.13U,  ALPffi9-»-(100-CC)  +  200 

.13  5.  p  ,      ROUTE  TO   NORMAL  APPROX.    FOR   CONF .    INT.    OF  LARGER  SAMPLE  SIZES 

136,  *((WMx2)>80)/Z;5 

^137,  p  COMPUTING  CONFIDENCE  INTERVALS  BI  EXACT  P-VALUE 

138.  CDF-*-P 

.13  9.  p  INDEX  POSITION  OF   VALUE  IN  CDF  S   ALPHA 

.11+0,  INDEXi^(.*/  (.CDF^ALPHA)) 

lUl.  ■*-{INDEX>0)/L& 

.142]  INDEX-*-! 

~9T 


■*-L6 
Q  COMPUTING  CONFIDENCE  INTERVALS   USING  NORMAL   APPROX.    W/C.F. 

L5  :  UALPHA-^ (.DENOMZ^  (NORMPTE  ALPHA  )  )+NM-0  .  5 
B  ROUND   UALPHA   DOWN  AND   INCREMENT  31   ONE 

INDEXt-lUALPBA  +  1 
L&iIPXi-NN.  INDEX 

CIi-IPX   CONFMW  A 

•  A    •  («CC),'  PERCENT  CONFIDENCE  INTERVAL   FOR  TEE  SHIFT  IN    • 

'  LOCATION   BETWEEN  POPULATIONS   X  AND   I  IS:' 

I  I 

[  ^  C  '  .(iPCICl]  ),'  S  Ma)    -   MU)   S  '  ,(5>CJ[2]),'  )• 

El: 'ERROR:    THE  SIZE  OF  JOUR  SAMPLE  IS  LESS  THAN  TWO;    TRI  AGAIN.' 

*B1 
EZ:' ERROR:    lOU  HAVE  ENTERED  MORE  THAN  ONE  VALUE.    TRI  AGAIN.' 
I  t 

-►S3 
E2:' ERROR:    lOU  HAVE  NOT  ENTERED  A   VALUE  OF  1   OR   2;    TRI  AGAIN.' 

I    T 

-52 
E^:' ERROR:    lOU  HAVE  NOT  ENTERED  A    VALUE  OF  1 ,    2 ,    OR   3 -,    TRI  AGAIN.' 

*(AA  =  2WN2 

V 


7  NPSLR:  N ;  SUMX  iSUMJ'.XBAR ;  IBAR ;  SUMX2  '.SUMXY  :3  ;A:WW:XX;3B:U:D:  ALPHA  :P:CC;  CDF 
; TALPUa  ;NN:CI:  SLOPES  \RR  :  SR-,TI  -.DENOM :  INDEX ;  FF  -.XiI-.QiE:  CHA  xPV 
p    PROGRAM   CONDUCTS  NONPARAMETRIC   LINEAR   REGRESSION.    THE   LEAST  SQUARES 
q  ESTIMATED   REGRESSION  LINE   IS   COMPUTED   WITH   HIPOTHESIS   TESTING  AND 
fl  CONFIDENCE  INTERVAL   AVAILABLE  FOR   THE  SLOPE  3 .    IF  3   DOES   NOT  LIZ   IN 
fl  THE   C.I.    AN  ALTERNATE  REGRESSION  LINE  IS  PROPOSED.    SUBPROGRAMS   CALLED 
fl  ARE:    SPMANP,    KEN D ALP ,    NORMPTH ,    INPUT,    AND   CONFLR. 

n  INPUT  DATA 

R-^INPUT   2 


'srax-f+zx 

XBAR'^SUMX*N 
IBAR'^-SUMI*N 

'sraX2-<-+/(X*2) 


ASSIGN  THE  SIZE  OP  X    {AND  7)   TO  N 
COMPUTE  TEE  SUM  OF  X'S  AND  I'S 

COMPUTE  TBS  MEAN  OF  X  AND  I 

COMPUTE  THE  SUM  OF  THE  X'S  SQUARED 
COMPUTE  THE  SUM  OF  X  TIMES  J 


SUMXI*-+/(X^7) 
p  COMPUTE    • S • ,  THE  SLOPE  OF  THE  ESTIMATED  LEAST  SQUARES  REGRESSION  LINE 

Bfi(.Ny.SUMXI')-{SUMX^SUMI))*(.(NxSUMX2)-(.SUMX*2)) 
P  COMPUTE    'A',    THE  I -INTERCEPT 

A->-YBAR-(.B^XBAR) 

FF*-  <  N  < 

'    THE  LEAST  SQUARES  ESTIMATED  REGRESSION  EQUATION  IS:    ' 

I  I 

J  ^  I  =    '  .(fA),'    +  '  .(.93), 'X.' 

'DO  lOU  WISH  TO  ENTER  SOME  X  VALUES  TO   GET  THE  PREDICTED  I"S?    il/N).' 

ww->-a 

*{WW='N' )/Ll  — 

L2: 'ENTER  X   VALUES. ' 

XXfQ 
p  CALCULATE  PREDICTED  I'S 

H*A  +S  X  XX 

'THE  PREDICTED   7  VALUES  ARE:     ',(911) 
I  I 

'  WOULD   lOU  LIKE  TO  RUN  SOME  MORE  X  VALUES?    (T/N).' 

^<.WW=  '7'  )/L2 
Ll: 'WOULD   lOU  LIKE  TO   TEST  EIPOTHESIS  ON  3,    THE  SLOPS  OF  THE 

,V ) .  ■ 

WW-'-a 

*{WW='N'  )/L3 

'ENTER   THE  HIPOTHESIZED   SLOPE.' 

BB-t-a 
n  COMPUTE  UI'S 

U*7-(BB^X) 

CHA*'>    ' 
p      CALL   SPMANP  TO   COMPUTE  RHO  AND  ASSOCIATED  P -VALUES. 

D*X   SPMANP   U 

*(DZ11>Q)/L11 

CHA- > <    ' 
Lll:PV-2^DL2l 


EQUATION?    (7/ 
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CS8 
59 

.60 
SI 

.52 

63 

6U 
65 
66 
67 

;63 

70 
71 
72 
.73 
C7U 
75 
75: 


PV-l 

L19:'  SPEAEMAN'  <S  R  EQUALS:    '.(HsCCl]) 
I  I 

<TSE  P-7ALUS  FOR  30:    3   =    <.(.is33),'    VSRSaS  31:    3    '.(.mCEACl    2J),(«55),'  IS: 
I  f 

'TBE  P-VALUE  FOR  TEE  TWO   SIDED  TEST  OF  EIPOTBESIS  IS:     '  ,  (UisPV) 
t  I 

Q     IF   USING  T3E  NEW  REGRESSION  EQUATION  BASED   ON  MEDIANS.    EXIT  3ERE. 
L3:*(.FF='I'  }/0 

fl  COMPUTE   CONFIDENCE   INTERVALS   ON   B 

'WOULD   lOU  LIKE  A   CONFIDENCE  INTERVAL   FOR  TEE  SLOPE?    (.I/N).' 

ww*a 

*(.WW='N'   )/0 
L1Q:CC'>-INPUT   5 
0        C3ANGE  ENTEEED   VALUE  TO   ALP3A 

ALP3A-<-(lQ0-CC)*200 
Q     ROUTE   TO   NORMAL   APPEOX.    FOR   CONF .    INT.    OF  LARGER   SAMPLE  SIZES 

*(N>12:/L5 


71 'a    p  COMPUTING   CONFIDENCE  INTERVALS  31   EXACT  P-VALUE 


C77, 

[78]  P-'rKENdALP  N 

■■-3  CDF*?Z1-: 

80  0  INDEX  POSITION  OF    VALUE  IN   CDF  £   ALPHA 

8 1  INDEX*-  ( +  /  ( CDFS.ALPEA  )  ) 
■82  ■*-(INDEX>0)/L6 
^83  INDEX-*-! 

8"*  -t-Le 

3  5  0         COMPUTING  CONFIDENCE  INTERVALS   USING  NOEMAL   APPEOX.    W/C.F. 

C35j  Zz:DENOM->-(  (<Vx(iV-l)x(  v2x.V)  +  5  )  )-!-13}*0.5 

, 8  7 1  TALPSA^DENOMx ( |  iMORMPTB  ALPHA  ) ) 

L88J  ■^L'i 

.89]  Zo:TALPHA*?Zl:INDEXZ 

;90j  L^:CI->rX   CONFLR   I 

[91]  NN-\^CI 

[92 J  SL0PES-1*CI 

[93]  RE->^l(.(.NN-TALP3A)*2') 

[9'a]  *(RRxO)/L2Q 

[9  5  J  RR'>-1 

[96 j  L2Q:SR-»-r  (1  +  ((NN+TALP3A)*2)) 

'97;  ^iSR^:2SL0PES)j/L2-L 

[98]  SR->-o  SLOPES 

99]  L2l:'    A    '(9CC)<    PERCENT  CONFIDENCE  INTERVAL   FOE  B,    THE  SLOPE  OF    ' 

100]   '  TEE  ESTIMATED   REGRESSION  LINE,    IS:< 

101  '  ' 

102]   J  ^  (  '  ,(9SL0PESLRR1).'    <   B   <    '  ,  (.vSLOPESLSRl) ,  '     ).' 

a   IF  B   OUTSIDE  THE  C.I.    CALCULATE  NEW  EQUATION  BASED   ON  MEDIANS 
a-{(.B^SLOPESiRR:i  )a  (B^SLOPESlSR:i  )  )/0 
p  ORDER  X  AND   I 

X-<-XCAX] 

i-fiun 

A     CHECK  TO   SEE  IF  TEE  SIZE  OF  SS  IS  EVEN  OR  ODD  FOR  FINDING  MEDIANS 

-<-((.2]NN)  =  0)/Sl 
p  COMPUTE  MEDIAN  FOR  ODD   CASE 

B-i-SLOPESZ(.(.NN+l)*2)2 

+52 
P  COMPUTE  MEDIAN  FOR   EVEN  CASE 

SI :  S-<- (SLOPES  C  (WW*  2  )] +SL(3PES  C  ( (A'iV+2  )  + 2  )])  + 2 
p  DO  THE  SAME  FOR  THE  X  AND  I  VECTORS 
S2:*((2TiV)  =  0)/S3 

IBAR->-lt{(N+l)*2U 

XBAR-^XtllN+l)*2n 

■*OUT 
S3:75Ai?-^(7[(W+2):+7C({W+2)*2):  )  +  2 

X3ili?-<-(XC(;i/+2)]+XC((A'+2)  +  2)]  )  +  2 
p  COMPUTE  NEW  INTERCEPT    'A' 

OUT :  A-^IBAR-  (B^XBAR ) 

' THE  LEAST  SQUARES   ESTIMATOR  OF  3   LIES  OUTSIDE  TEE  CONFIDENCE  INTEEVAL . « 

'DISCARD   TEE  LEAST  SQUAEES   EQUATION  AND   USE:' 

I  I 

I  ^  I  =  '  ,i9A),'    +  '  ,(.9B),'X    • 

'     THIS  EQUATION  IS  BASED   ON  THE  MEDIANS  OF  TEE  X  AND   7   DATA   AND    ' 

'THE  MEDIAN  OF   THE  TWO-POINT   SLOPES   CALCULATED    FOR   THE   CONFIDENCE   INTERV 

AL    JN   3. ■ 
^132:   '  ' 

,13  2:  p     ALLOW    USSR    TO    DO    SOME   ANALYSIS   ON   NEW   EQUATION 
Cl34J   'DO   lOU  WISH  TO   ENTER  SOME  X    VALUES   TO   GET  PREDICTED   I"S  FROM  TEE  NEW  E 

QUATION?         {J/N).' 

*(PP='J' )/L2 

FF*<I' 

*L1  -     .  .. 
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7  SIGN '.AiC:3;D:PVAL',X;M0;N '.CDF: ALPHA  :CI:I:AA:BB:CC',DD',P7;PVI',NNN',KFOSiORD 
D ;  11  lORDX  iKALPaA  ;  Z  :  Zl ;  QUA  :  -t^WiPVMiPVZ  ;RiQ 
ri]   p       i'HIS   JUNCTION   USES   THE  ORDINARY   SIGN  TEST  TO   CALCULATE  TEE  K 
^2'   fl  STATISTIC,    P-VALUE,    AND   CONFIDENCE  INTERVAL   AS   A   TEST   FOR  MEDIANS. 
'3"   fl  THE  LAST   OPTION   WILL   DISPLAY   A    TABLE   OF   CONFIDENCE  INTERVALS   OF 
'u'   0  ORDERED   STATISTICS   WITH   CONFIDENCE   COEFFICIENTS. 

'5'  a    SUBPROGRAMS   CALLED   BI  THIS  FUNCTION  INCLUDE:    BINOM,    NORMCDF ,    NORMPTH , 
^6'       fl  INPUT, AND   QUANC. 
'j'  II 

"8"    •        DID   lOU  ENTER   THIS   PROGRAM  FOR   TEE  SOLE  PURPOSE  OF   GENERATING  CON 
FIDENCS  INTERVALS   FOR  A   SPECIFIED   SAMPLE  SIZE  AND   QUANTILE?    il/N).' 

■g]  ww^a 

"lli   '        THE  NULL   HYPOTHESIS   STATES   -   THE  POPULATION  MEDIAN    (.M)   IS   EQUAL   T 
;      0   TEE        HYPOTHESIZED   MEDIAN    (.MO);    HO:    M   =   MO.     ' 

"13 ■   •  WHICH  ALTERNATIVE  DO   YOU  WISH  TO   TEST?' 

lu"   '  ' 
15"  53:'  ENTER:         1  FOR  HI:    M  <  MO;         2   FOR  HI:    M  >  MO;  3  FOR   ffl :  M   s  MO.' 

15'   C-f-Q 

1?:   -»-((Cs1)a(Cs2)a(C=:3  )  )/£:3 

FOR   SINGLE-SAMPLE  PROBLEM;         2   FOR  PAIRED-SAMPLE  PROBLEM. ■ 


isl  32: 'ENTER:      1 
"igl  AA-i-D 
20]  *{(AA  =  1)'-^(AA' 


)a(AA=2))/E2 

'21]  ■>'(AA  =  2')/L9 

'2  2]  fl  INPUT  DATA   FOR   SINGLE  SAMPLE   CASS 

'2  3]  X-*rINPUT   1 

24]  NNN-^oX 

25]  MO-^INPUT   3 

25]  D^X-MO 

^271  *L11 

r23j  fl  PAIRED   SAMPLE   CASE 

l29l  L9:R-(-INPUT   2 

l3  0]  5-t-lti? 

[31]     :i*-i-^iQ-n)fS 

C32]  J-<-(Q+l)+i? 

[33]  DD-rX-Y 

[3u1  NNNiroDD 

r3  5]  MO-^INPUT   u. 

[35]  Z}t-(X-y)-MO 

^37^  fl                      COMPRESS   D   TO   REMOVE   ZEROS 

38]  ill:^-f-(Z3xO)/D 

.3  9.  fl                  RECORD  LENGTH  OF  A  AND  ASSIGN  TO  N 

M-O.  N*pA 

yi,  P                                                 KEEPING  TRACK  OF  POSITIVE  SIGNS 

>2  iiTOS-«-+/(i4>0) 

>3,  fCN^SOl/NORM 

>u  PVAL->-BINQM  N 

.H5.  *(XP0S>0)/P1 

>6  P7I-<-l 

>7  ■*P2 

>8,  Pl:P7J-<-l-P7^LCZP0S] 

.U9.  P2:P7-<-P7iilLC(XP0S  +  l)] 

,50,  *L& 

.51.  a    IF  N  IS   GREATER  THAN   30  J7S^  NORMAL   APPRQX  W/   CONTINUITY   CORRECTION 

.52.  WOPM:Zl-f-((XPOS-0.5)-(0,5xiV))+f0.5x(iV*0.5)) 

.53,  Z-<-((i!:POS+0.5)-(0.5xA^))+(0.5x(JV*0.5)) 

,54,  PV-f-NORMCDF   Z 

5  5  PVIf-1- (.NORMCDF   Zl) 

,56.  p  JP  PAIRED  SAMPLE  TEST  GO  TO   LI 7  FOR  OUTPUT  STATEMENT 

.57,  Z;6:P73-t-2x(L/(P7,P7J)) 

,58.  -►(P7351)/iV5 

,59  P73-<-l 

.50.  iV5:P7W-f-(3,l)p(P7,P7J,P73) 

,61,   'COMPUTATIONS  ARE  BASED   ON  A   SAMPLE  SIZE   OP:  •  ,  (»iV) 

,62,    •  ' 

,53,   'THE  TOTAL   NUMBER  OF  POSITIVE  SIGNS  IS:     ',(.9KP0S') 

.6U  I     ■ 

.55  ^(AA=2)/L17 

.66]       'THE  P-VALUE  FOR  BO:    H  =    '.(suMO),'    VERSUS  El:    M    '  ,  (vLOGICZCiH) ,'    '.(sMO) 

,'  IS: ',(u»P7M[C:l]) 
[67]   f  ' 

6  8  "^L 1 8 

;69;  L17: 'THE  P-VALUE  FOR   COMPARING  THE  MEDIAN  OF  TEE  POPULATION  OF    < 
.70.   'DIFFERENCES   TO   THE  HYPOTHESIZED   MEDIAN,' 
71    '  ' 
C72]   'ffO:  M(X-I)    =  '.(aMO),'  VERSUS  31:    M(.X-I)    '  .  i^LOGICZC;lj  )  ,  '     ',(«M0).'.  IS 

:  '  ,  >«P7MCC;1J  ) 
'73  J    '  ' 

;7u;  LIS: 'WOULD   YOU  LIKE  A   CONFIDENCE  INTERVAL   FOR  THE  MEDIAN?    (Y/N).' 
.75    flB-«-a 
.76.  ^(.BB='Y')/L16 
,77,  EQUANT 

.7  8.  p  INPUT  SIZE  OF  CONFIDENCE  INTERVAL 

.79.  L16:CCi-INPUT    5 
,80.   >lLPff>?-<-(100-CC)t200 
.81.   f<.NNNi30)/NOHMl 

.8  2.  p  COMPUTING  CONFIDENCE  INTERVALS  BY  EXACT  P-VALUE 

.83.   CDF-'rBINOM  NNN 

.84.  p  INDEX  POSITION  OF  CDF  FOR  ALPHA    *    2 

.8  5    B*->r/(.CDFiALPBA) 
:86]  ^{B>0)/SKIP 
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5-f-l 

*SKIP 
p  COMPUTING  CONFIDENCE  INTERVALS  BI  NORMAL  APPROX . 
NORMl : KALPHA->-(  (.0  .  5^  (NNN*0  .  5  )  )x(NORMPTH   ALPHA  )  )  +  (0  .  5xNNN)-0  .  5 
q  ROUND   KALP3A   DOWN  TO  NEAREST  INTEGER  AND  INCREMENT  31   ONE 

B-^LKALPBA  +  l 
0    _     _  IF  SINGLE  SAMPLE  CASE   GO   TO   LI 

p    ■     ~      CALCULATE  AND   PRINT  OUT  CONF .    INT.    FOR  PAIRED   SAMPLE   CASE 
L5:0RDD-^DD£&DD1 

•J'S-^-oORDD 

CI'^0RDDi31.0RDDr(.I7-(.B-l'))l    '  ^„  „„„  . 

'A    •  (5CC).'  PERCENT   CONFIDENCE  INTERVAL   FOR  TEE  MEDIAN  OF  TEE    ' 
I   ' POPULATION  OF  DIFFERENCES  IS : ' 

'  <  ',(»CJC1]),'  S  MEDIAN(.X-I)   S  '  ,(»CIC2]  ),'  )' 

-^UANT 
q  CALCULATE  AND  PRINT  OUT  CONF.    INT.    FOR   ONE   SAMPLE  CASE 

L7:0RDXi-XZAXl 
:i07]  IJi-paRDX 
108]  CI*0RDXL31,0RDXZ(I7-(.B-1))1 
10  9]   '.;  '  ,:sCC),'    PERCENT  CONFIDENCE  INTERVAL   FOR   TEE  MEDIAN  OF  TEE  POPULATIO 

N  IS:     ' 
'110]   '  ■ 

llli;   •  (  '.(iBCJCl]  ),'  ^   MEDIAN  ^    •  ,(!iiCJ[2]),'  )' 

112    '  ' 
!ll3]  QUANT'.'WOULD   lOU  LIKE   CONFIDENCE  INTERVALS   FOR  A   SPECIFIED   QUANTILE?    (.I/N 

;iiiAj  VA^^'T! 

115j  ■*Chrn/='I'  )/31 

115j  *0 

!ll7]  3»-:' ENTER   DESIRED   SAMPLE  SIZE    (.SINGLE  INTEGER   VALUE).' 

11 8  J  NNN*-: 


■*-((QNNN)>l)/En 
♦Ul ""■  '  ^  ■" 


.  .i.\NNN)vO)/E^ 
31 -.'ENTER  DESIRED   QUANTILE;    FOR   EXAMPLE:    ENTER    20,  FOR   TEE   20TE   QUANTILE. 


C122]  QUA'>-a 

[I22I  *((QUAiO)'^(gUA>lQO))/El 

"l2'i;  NNN    iUANC   QUA 

LI253   '  ' 

[126 J   '  *****  TEIS  TABLE  GIVES   CONFIDENCE  COEFFICIENTS  FOR   VARIOUS  INTERVALS  W 

ITE    ' 

'ORDER  STATISTICS  AS  TEE  END  POINTS  FOR  TEE    ' .("QUA^.'TE   QUANTILE.^ 

I  I 

El: 'ERROR:    TEE  QUANTILE  VALUE  MUST  LIE  BETWEEN   0  AND   100;  TRI  AGAIN. ^ 
I  I 

*B1 
E2:' ERROR:    700  HAVE  NOT  ENTERED  A   VALUE  OF  1   OR   2;    TRI  AGAIN.' 
t  t 

■*-B2 
E3i 'ERROR:    lOU  EAVE  NOT  ENTERED  A    VALUE  OF   1,    2,  OR    3;  TRI  AGAIN.' 
I  I 

♦53 
E'*:' ERROR:    lOU  EAVE  NOT  ENTERED  A   SINGLE,    INTEGER   VALUE:    TRI  AGAIN.' 

♦flu 
7 


7  SPEARMAN '  X '  I '  A'  Q'  R'  CEA  •  PV 
Cl]   R       this' FONCTlof)   COMPETES  TEE  SPEARMAN  R   STATISTIC  WEICE  MEASURES 
^2    P  TEE   DEGREE   OF   CORRESPONDENCE   BETWEEN   RANKINGS   OF   TWO    SAMPLES.    THE  P- 

3    p  VALUE  IS   GIVEN  FOR   TESTING  ONE  AND   TWO-SIDED   EIPOTEESIS   OF  ASSOC. 

U    fl  SUBPROGRAMS   CALLED   31   TEIS   FUNCTION  INCLUDE:    TIES.    TIESK,    SPEARP , 

5  fi  INPUT,    SPAPROX,    INTERP ,    AND   TEE   VARIABLE   PMATSP. 

6  R*INPUT   2 
"•^         Q-^-l^R 

X*l*(Q-i-l)fR 


[127] 

128 

129, 

130 

131 

132 

.133. 

13U 

135 

136 

137 

138 

139 

,140, 

.mi. 

i¥ 


r-<-(o+i)+i? 

'ILL    


,_0]  q  CALL    SPMANP   TO   CALCULATE   TEE   STATISTIC  AND   ASSOCIATED   P-VALUES 
1.1 1 2  A-f-X    SPMANP    I 

Ci2j  *<^ci:>o}/z;i 

^ !  C3A* ■ INDIRECT ' 

■>'L2 
LI :CHA*' DIRECT' 
L1:PV^2-^A121 

*(P751)/L3 

P7-<-l 
LZ:  'SPEARMAN'  'S  R  EQUALS    '  ,  (45/1  Cl]  ) 


14 
15 
15 
17 
18 
19 
20 
21 
[22] 
23 
24 
25] 


TEE  P-VALUE  FOR  EQ:    NO   ASSOCIATION   EXISTS   VERSUS' 

Hi:     ',(.9CHA),'    ASSOCIATION  EXISTS  IS:' ,(^9AZ22) 

TEE  P-VALUE  FOR  TEE  TWO-SIDED  TEST  OF  EIPOTEESIS  IS:'  A^vPV) 
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7  h7ILCOX;AiB-,DiE-,F:PV2'.Zl',ZiDEN',NUMZiNU}JZl',PVAL:X:MO-,N',TPLUS:CDF:TALPHA',A 

LPHA:H:CIiIiAAiBB',CC;NN;DD-,PV;P0S;TP0S;NM;FVI',TF0Sl;NNN',C;PVM;PV3;H;Q-,I 

RAP : NUM : DENT: DENTl : TC ; TC\ 
fl  TEIS  PROGRAM   USES   THE  WILCOXON   SIGNED   RANK  TEST  TO   CALCULATE  TEE      T 
fl  STATISTIC.    P-VALUE,    AND    CONFIDENCE   INTERVAL   AS   A    TEST   FOR   MEDIANS. 
p  SUBPROGRAMS   CALLED   BY   TEIS   FUNCTION  INCLUDE:    TIES,    WILP ,    NORMCDF , 
p  NORMPTE,    CONFW ,    AND   INPUT. 
I  I 

•        TEE  NULL   E7P0TEESIS   STATES    -    TEE  POPULATION  MEDIAN    (3?)  IS  EQUAL   T 
0   TEE        E1P0TEESIZED  MEDIAN    (MO)-.    EQ:    M   =   MO.' 


3]  '              WEICE  ALTERNATIVE  DO   lOU  VISE  TO   TEST?' 

9]  •  ' 

;iO]  B3:'  ENTER:         1   FOR  El:    M  <   MO;         2   FOR  El:    M   >  MO;           3  FOR  HI:    M   i  MO.' 

11  C-*-Q 

;i2;  *((Csl)A(C*2)A(C3:3))/ff3 

13]  B2: 'ENTER:         1    FOR   SINGLE-SAMPLE  PROBLEM;         2   FOR  PAIRED-SAMPLE  PROBLEM.' 

■"^^       •>-((AA£l)^(.AAx2))/E2 
^(AA  =  2WL<5 
p  INPUT  DATA   FOR  SINGLE  SAMPLE   CASE 


15 
J.16. 

£l7. 

.18; 
19; 
2o; 

21' 

22 

233 

25^ 
26 


X^INPUT   1 
NNN^oX 
MO^INPUT   3 
D^X-MO 
->-Lll 
fi  PAIRED   SAMPLE  CASE 


L<5:R-<rINPUT 
Q-e-l+J? 
X-(-l  +  ((3+l)+i? 

•27]  Y^'kQ^D^R 

[28]  DD^X-'I 
NNN-'-qDD 
MO-^INPUT   U 
D*{X-I)-M0 


,29 
30 

,31 
32 
33 
3U 

;35 

,35 

!37] 
38" 
39 
^■0 

.42 
.43 

45 

!47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 
,60 
.61 

:62 

63 
64 
55 
66 
67 
58 
59 
70 
71 
72 


:z^i 

•M: 

77 

78; 

.79. 

:   • 

84' 

;85: 

.85, 
.87. 

COMPRESS  D   TO  REMOVE   ZEROS 
Lll:A-^(.DxQ')/D 
P  RECORD   LENGTE   OF  A   AND  ASSIGN  TO  N 

Nt-pA 
p  KEEPING  TRACK  OF  POSITIVE  SIGNS 

POS-fU>0) 
p  TAKE  TEE  ABSOLUTE    VALUE   OF  A;    ASSIGN  TO    3   AND   ORDER   3 

0       REORDER  POSITIVE  SIGNS  TO   COINCIDE  WITE  PROPER  POSITIONS  IN  B 

pos'>rPosuuAn 

fl  CALL   FUNCTION  TO   BREAK  TIES 

Ef-1   TIES  B 
p       CALCULATE  TPOS   BI  ADDING  ACROSS  ALL  POSITIVE  VALUES  OF  E 

TPOS^+/{POSy-E) 

TPOSl-^TPOS 
fi  GIVES  SIZE  OF  LEFT  TAIL   OF  PROBABILITI  DISTRIBUTION 

WM-«-(L((+/iiV)T2))  +  l 
fi  GO   TO   STATEMENTS  BASED   ON  LENGTE  OF   VECTOR   E 

-»-(W>16)/L3 
fi  GENERATE  NULL   DISTRIBUTION  FOR   TPLUS 

FirWILP  N 
fl      IF  TPOS   PALLS   IN  LEFT  EALF  OF  PROB   DIST  CALCULATE  PVALUE  AS  NORMAL 

^(.TPOS^iNM-l  )  )/TREG 
p  OTEERWISE   USE  TEE  NEGATIVE  T  STATISTIC 

TPOS'^(  +  /\N)-TPOS 
p  IF  TPOS   IS   FRACTIONAL   USE  BOTE  TEE  INTEGER  ABOVE  AND   BELOW  AS  TPLUS 
TREG:-*((1\TP0S)  =  Q)/N0N 

TPLUSf-lTPOS 
Pl:PV'>-l-{iFLTPLUSl+FLTPLUS+ll  )  +  2) 

PVI--(.FZTPLUS  +  11+FLTPLUS  +  21')*2 

■*CEECK 
NON:-*-(TPOS>0)/GO 

PV->-l 

■*P2 
GO:PV-*-l-FZ(.TPOS)'\ 
P2:PVI->rFllTP0S  +  \)l 
CEECK:*(.TP0S1^(.NM-1))/L6 

PV2->-PV 

PV-^PVI 

PVI*PV2 

p  °OMPnTE   NORMAL    APPRO X.    W/CONTINUITI   CORRECTION   FACTOR 
L2  -.TRAP-^kN^  (:V+1  ;  J-i-y. 

NUMZ-^iTPOS  +  O  .  5  )-TRAP 

NUMZl-r{TPOS-0  .5)-TRAP 

Z)£'A'-<-(  (A'x  (A?  +  l  )x  (  (2xA^)  +  l  )  )  +  24  )*0  .  5 

Z^NUMZ^'DEN 

Zl^NUMZHDEN 
fi  COMPUTE   STUDENT  T  APPROXIMATION  WITS   CONTINUITY   CORRECTION  FACTOR 

NUmA  (TPOS -TRAP) 

DENT^U  (A'x  (DEN*2  )  )*(N-1  ))-((  (NUM-0  .  5  )*2  )*(N-1)  )  )*0  .  5 

DENTl^ll  (A'x  (DEN->'2  ))*(N-l))-((  (NUM+0  .  5  )*2  )+  (iV-1  )  )  )*0.  5 

TC*(NUM-0.5)tDENT 

TCl-^(NUM+0  .  5  )*DENT1 
0  COMPUTE  AVERAGE   OF  TC  AND    ZC 
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[88 
89 
90 
91 
92 
93 
94 
95 
96 
97 

99] 

.loo; 

101. 
102. 

103. 

!iou. 


*{T?0S^U+/\N')-i-2))/ SECOND 

PV*C(>1-(N0RMCDF   Zl  ))+(!-( (iV-1)  TDISTN   rci)))*2 

PVI-^-UnORMCDF   Z)  +  iiN-l)    TDISTN  TC))*2 

■>-L6 
SECOND :PV->-(Cl-(NORMCDF   21))  +  ((iV-l)  TDISTN  TC1))*2 

?VI*{{NORMCD?   Z)  +  (l-{(iV-l)  TDISTN  TO))*! 
L6:PV3-<-2xa/(P7,P7J)) 

*(P73Sl)/iV5 

P73-*-l 

iV5  :P7M-<-(3  ,1  )o  (.PVI  J>V .PV2  ) 

'  COMPUTATIONS   ARE   BASED   ON  A   SAMPLE  SIZE  OF:     '  ,  (asA^) 
t     I 

'TEE  TOTAL   SUM  OF  POSITIVE  RANKS  IS:     ',(5rP051) 

fl  IF  PAIRED   SAMPLE  TEST   CO   TO   L17   FOR   OUTPUT  STATEMENT 

-■(.AA  =  2)/L17 

'TEE  P-VALUE   FOR  EO:    M   =    '.(sMO),'    VERSUS  HI:    M    '  ,  CfLOGICLC-.ll  ) .  '     <  .  (.tiMO 
),;    IS:  <  .(u«P7MCC:l] ) 

-519 
L17  '.  -THE 


■VALUE  FOR   COMPARING   TEE  MEDIAN   OF  TEE  POPULATION  OF 


'DIFFERENCES  TO   THE  EIPOTHESIZED   MEDIAN,' 


S: ' ,(u»P7MCC;lJ ) 


VERSUS   31:    !^CA-Z)    '  , '.vLOGICZC;!:  ) , 


.  i^MO), 


Lis: 'WOULD   lOU  LIKE  A   CONFIDENCE  INTERVAL   FOR  TEE  MEDIAN?    (.I/N).' 

BB*a 

*(.3B='I'  )/Lie 


;l. 


Z-a-.CC^iy'U'"    5 

ALPHA-^-^lOO-CC^^lQQ 
q         ROUTE   TO   NORMAL   APPRO X. 

-(jVW7/>15  )/Z~ 
a  COMPUTING  CONFIDENCE  INTERVALS   31  EXACT  P- 

CDFi-f/ILP  NNN 
0  INDEX  POSITION  OF  CDF  FOR  ALPHA    *    2 

TALPHA'>-  ( +/  ( CDF <. ALPHA  ) ) 

■*{TALPHA>Q)/JUMJ> 

TALPHA-^1 

■'JUMP 
a  COMPUTING   CONFIDENCE  INTERVALS   31   NORMAL 

Lii:DEN0MZ-^((2^NNNx(NNN+l  )x  (.i2^NNN  )  +  l)  )^3  )*0.  5 

TALPHA-t-C  ( (DENOMZx  (NORMPTB   ALPHA  )  )  +  (.NNNx  (NNN+l )  )-2  )xU  ) 
0         ROUND   TALPBA   DOWN  TO   INTEGER    VALUE  AND   INCREMENT  BI  ONE 

TALPEA-t-l  T ALPHA  + 1 
p  IF  ONE  SAMPLE  CASE  GO  TO   L7 

JUMP:*UA  =  1)/L7 

fl  CALCULATE  AND  PRINT  OUT  CONF .    INT.    FOR  PAIRED   SAMPLE  CASE 

Ll:CI-^TALPHA   CONFW  DD 

'A    'A^CC).<    PERCENT   CONFIDENCE  INTERVAL  FOR  TEE  MEDIAN   OF  THE 

'POPULATION  OF  DIFFERENCES  IS:' 


FOR   CONF  INT  OF   LARGE  SAMPLE  SIZE 
■VALUE 


APPROX.    ?//C.F. 


I 


(  '  ,(5CIC1]),'  £  MEDIAN(.X-I)   ^    •,(5CJC2]),'  )• 


1U5 

r   1 

l'*6 

1 

1U7 

I   1 

1U8 

*0 

1U9 

E2:  ' 

ISO 

1     1 

151 

*B2 

152 

E3:  ' 

153 

1  t 

15"+ 

-►S3 

p  CALCULATE  AND  PRINT  OUT  CONF.    INT.    FOR  ONE  SAMPLE  CASE 

L7:CI-*-TALPHA   CONFW   X 

'A    ' ,{9CC),'    PERCENT   CONFIDENCE  INTERVAL  FOR   TEE  MEDIAN  OF  TEE  POPULATIO 

N  IS:' 


(  ' ,(«CICl]),'  ^  MEDIAN  S 
ERROR:  10 U  HAVE  NOT  ENTERED  A  VALUE  OF 
ERROR:    lOU  BAVE  NOT  ENTERED  A   VALUE  OF  1,    2,    OR   3 i    TRY  AGAIN.' 


',(«>CJC2]),'  )« 

1   OR   2:    TRY  AGAIN. • 
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APPENDIX  F 
LISTINGS  OF  SUBPROGRAMS  BASIC  TO  BOTH  WORKSPACES 


V  CBN-^BINOM  N:NiP:XiK:CDF 
p  THIS   FUNCTION   iS   A   SUBPROGRAM  OF  TEE  SIGN  TEST    (.SIGN).      IT   CALCULATES 
fl  THE   CDF  OF  TEE  BINOMAL   WHEN  PROBABILITY   -    .5.  N=SAMPLE  SIZE. 

Jf-f-o  ,  \N 

CBN-^-^\(.KlN)x(,Pi,K)^((.l-?)*N-K) 

7 


7  UfL   CBIN  R 

2  P    THIS   FUNCTION  IS  A   SUBPROGRAM  OF   CONFIDENCE  INTERVAL   GENERATOR  FOR 

3  p  THE   QTE   QUANTILE.      IT  RETURNS   THE   VALUE  OF  THE  BINOMIAL   CDF  AT  R, 
Ifl   Q  WITH  N ,P-L   WHERE  N=SAMPLE   SIZE  AND   P-PR0BABILIT1 . 

7 


7  CON->rXX  CONFLR   7l;3B  iSS',AA:A:XR:IR:B ',C;S 

;i]  p       THIS   FUNCTION   IS   A   SUBPROGRAM   OF  NONPARAMETRIC  LINEAR   REGRESSION 

2 J  P  (NPLR).      IT   CALCULATES   THE   TWO-POINT   SLOPE   FOR    EACH   PAIR   OF   POINTS 

.3  J  P  (XZ.II)   AND    ^XJ.IJ.j    ALL    K-J   AND   XI-XJ .    ALL   SLOPES   ARE   ORDERED   AND 

.u1  p  USED   TO   FIND   THE   CONFIDENCE^ INTERVAL   FOR   3.    THE   SLOPS   OF  ESTIMATED 

5]  p  EQUATION.      XX=   X  DATA  SAMPLE  AND  11=   I  DATA  SAMPLE. 

6J 

7]  p  RECORD  TEE  SIZE  OF  XX  AND  INITIALIZE  VARIABLES 

8]  BB^pXX 

9]  SSf-pO                                                  o 

.10]  AA*0 

11  p  THIS  LOOP  COMPRESSES  XX  AND  II  DOWN  TO  WEERE  TEE  XX  <  ALL  OTEER  XX 'S 

12.  Ll'.AA-^AA+l 

13  A-r(.XXUA2<XX'i 

lU  XR-^A/XX 

15  IRlrA/Il 

16  B^pXR 

17  -►(S  =  0)/L3 
.18,  C-f-0 

19.  p  THIS  LOOP  CALCULATES  TEE  SLOPE  OF  FACE  PAIR  OF  PAIRED   DATA. 

.20.  Ll:C-<-C+l 

.2l1  S-^illZAAl  -IRLCD*  (XXLAAl  -XRLCl  ) 

,22]  SS-^SS.S 

.23.  ->-(.C<b)/L1 

^2U  L3:->-(.AA<BB)/L2 

.25]  CON->-ipSS'),(.SSLiSSl) 
7 


7  CONFU-'rAA   CONFMW  BB-.A-.B  ;C:D:E:F:C:H 

[11  P       THIS   FUNCTION   IS  A   SUBPROGRAM  OF   THE  MANN-WHITNEI  TEST    (MANW). 

.2]  ft    IT   COMPUTES   CONFIDENCE  INTERVAL   ENPOINTS   FOR   THETA,    THE   SHIFT  IN 

.3 J  P  LOCATION   BETWEEN   X  AND   I.    AA=INDEX  POSITION   OF   C.I.    ENPOINT . 

,u1  p  3B=C0MBINED   DATA   SAMPLES. 

,5]  r    assign  size  of  x  vector  to  a;  index  position  for  conf  int  to  b. 

,=>]  a--aa[,'l:i 

C7]  3^AAi2l 


10 
.11 
.12 
.13 
.l** 

15 
.16 

17 

18 
.19 
,20 
.21 


p  ASSIGN   X   VECTOR   TO   C;    I   VECTOR   TO    D 

C-^A-f'BB 

D^-A^BB 
ft  REORDER  X  AND   I   VECTOR   VALUES  TO  ASCENDING  ORDER 

C+CCAC] 

D-^DLlDl 
ft  INITIALIZE   VECTOR   E  AND   VARIABLE  F 

E-*-pO 

F*l 
ftlNNER  AND   OUTER   LOOPS   CALCULATE  ALL   POSSIBLE  DIFFERENCES:    EVERY   I 
ft    ELEMENT  MINUS   EVERY   X   ELEMENT.    H    VECTOR   STORES   THESE  DIFFERENCES. 
L2:G*1 
Ll:£:-DCf]-CCC] 

H'^H,E 
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[22]  G*G*1 

25  -^(i'5(pU))/I2 

26.  B         REORDER  E   VECTOR   VALUES  TO  ASCENDING  ORDER 

27  H-rHLlBl 

28  fl  JiViJFX  APPROPRIATE  CONF  INT   V ALOES  FROM  B 
:29]  COiVFM"«-i?CS],flC((pff)-(S-l))] 

7 


V  CONFt-A   CONFW  B;C;D;E:C;H :NiF',A:B 
ft      TBIS  FUNCTION  IS  A   SUBPROGRAM  OF  THE  WILCOXON  SIGNED  RANK  TEST 
p  iWISIG).    IT  PROVIDES   CONFIDENCE   INTERVAL   END   POINTS   BASED   ON   THE 
n  AVERAGES   OF  ALL  PAIRS   OF   DATA    SUCB   THAT  ALL   XI   S  XK.    A=   INDEX 
n  POSITION  OF   C.    I.    END  POINT  AND   B=   INPUT  DATA' 

C-^-BL&Bl 
a  START  E  ACCUMULATION   VECTOR   OFF  WITB   ORGINAL   VECTOR      VALUES 

if-^2  0 

H*-B,C 

D-^Q 
a  OUTSIDE  LOOP   INCREMENTS   D   AND   RESETS   E  TO   D 
L2  :D->^D  +  1 

E*-D 
fl  INSIDE  LOOP  GENERATES  NEXT  SET  OF  AVERAGES  AND   CONCATENATES  TO   ORIGINAL 

VECTOR 
Ll:E*E+l 

F-^(CZD1+CZE1  )*2 

a    "■       CONTINUE  INNER  LOOP   UNTIL   E  EQUALS  TEE  SIZE  OF   C 

*{E<(.pC))/Ll 
a        CONTINUE   OUTER   LOOP   UNTIL   D   EQUALS   TEE  SIZE  OF  C  LESS   ONE 

-^(D<U?C)-l)j/L2 
fl  ORDER  FINAL  ACCUMULATED    VECTOR  B 

B^EihEl 
fl  INDEX   CONF  INT   VALUES  OUT  OF  B 

CONF*BZAZ,Bli{pB)-U-l)n 

y 


7  Px-DF  FDISTN   X;A',M:N:RM:RN',LM',LNiSM:SN:M2',N2 

11  P     TBIS   FUNCTION  IS   SUBPROGRAM  OF   THE  KRUSKAL-WALLIS   TEST    (KRWL)    AND 

2  P  THE   STUDENT  T  DIST ATDISTN ) .    IT  CALCULATES  APPROX.    CUMULATIVE  PROBS. 

3  a    OF  X   USING  TEE  F  DISTRIBUTION  W/   DF=(.M,N)   DEGREES  OF  FREEDOM. 

5!  P    TAKES   VECTOR  ARG. 

6  A-f-.M'iXiN+X^Mi-lfDF.N't-l^DF 

7,         -L^\(M>2)^N>2 

8  p   TREAT  TEE   1ST.    FOUR  SPECIAL   CASES. 
,9]         ■*s'L<  .    1    0    9M.N 

,10]  Lll:P-*-{    10il*O.S)  +  o  +  2 
^11.       ■►O 

12  L12:P*il*0.S 
13J   -►O 

^lu]  L21:P-«-l-(l-il)*0.5 

15  ♦O 

16  i22:P*A 

17  *0 

18  p  BEGIN  THE  GENERAL   CASES.    INITIALIZE  TEE  QUANTITIES. 


19  L:RM->-(l\Mi2).0/RN->-l\N^2 

20  SM*(RM=0)xRMiLMt-lM2-*-0.'. 

21  SN'>-lRN=0)xRN^LN'>-lN2'>-0.'. 


xRN^LN--lN2'^0.5^N-2 

22]  P  TREAT  SN .CASE  N-EVEN  FIRST. 

23.  -^NODD^W^iAN 

2U  ■*MEVEN^\Q=N2 

,2  5,  SA'-«-l+  +  /x\(l-A)o.xi  +  (-l+0.5xM)  +  iiV2 

'26  -^MEVEN 

.27  N0DD:->-ME7EN^\l=N 

28  SN'>-(1-A)*Q.5 

29  p    TREAT  TBE  PORTION  OF  SN  TEAT  DOESN^T  CEANGE  W/  M  ODD   OR  EVEN. 
,30,  *(2+l+DLC)x\0=LN 

31]  SiV-oSiVxi-u^/xN  (l-.4)o.xl+2xM2*l+2xiLi7 

.3  2"  a    IREAT   THE  M-EVEN   SUBCASE 

'33]  *<■  iM=2  )  ,l  =  2\M)/ME7EN,3+lfaLC 

3  5l  ■'■MEVEN 

36.  p    TREAT  TBE  M-ODD   SUBCASE 

,37,  5W-«-SWx2  +  ol 

3  8  ♦FWZJxiirM 

,39,  SiV-^Sil/xx/i++    l+2xt0.5+M2 

UO  MEVEN:*M0DD^\1=2\M 

,U1,  -<-ENDx\2=M 

,U2,  SM-»-l  +  +  /x\>lo.xl-  +  2xiM2 

>3  *END 

>4.  p    NEXT  TREAT  TBE  SPECIAL   CASES   FOR   ODD-M. 

,U5,  MODC:SM-<-(OxM=l)  +  ((A*o.5)  +  o  +  2)>:3=M 

>6,  *EiVDxi(i=M)v3=M 

,47,  SM*(A*0.5)x(♦OT2)xl++/x\Ao.xl-  +  l+2x^i;?f 

.U8]  END:P*(.SN^A*M*2)  +  i>*xRNxRN-RM')+C(.RMxRNx    lOil*0.  5  )  +  o  +  8  )-2x^iVxSMx  (l-A  )*0  .  5 
7 
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7  XX-^N  INDEXPLS   3 ;AA:3B ;CC:DD',X1-,NN',NC;N'.B :CiD',E:F 
n      THIS   FUNCTION   IS  A   SUBPROGRAM  OP  WRUSKAL-WALLIS    iKRWL) ,    AND   MANN- 
Q  WHITNEY    (MANW).      IT  COMPUTES   A   MATRIX   OF  RANKS   OF   THE   ORIGINAL   DATA 
a   WHERE  THE   FIRST  ROW   OF  THE  MATRIX   IS   THE  FIRST  SAMPLE.      IF  TIES   OCCUR 
n  IN   THE  RANKING   WITHIN  A   SAMPLE ,    SUCCEEDING   TIES  ARE   INCREMENTED   31 
p  ONE.    N=    SIZES  OF  THE  SAMPLES   BEING  PASSED   IN  3 .    B=   ALL   DATA   SAMPLES 
fl  COMBINED. 

fl        RECORD  SIZE  OF  N  TO  DETERMINE  NUMBER  OF  SAMPLES  IN  B 

AA-*-pN 
fl  LOCATE  LARGEST  SAMPLE  SIZE 

BB^f/N 
p  ORDER   3   SMALLEST  TO   LARGEST 

DD^BZiBl 
p        SET   UP  MATRIX   OF  SIZE  REQUIRED   TO   STORE  SAMPLE  RANKING. 

XX-'^UA.BB^pO 
p        CONCATENATE   0  AND  N 

NN-^O.N 
p      FIND   CUMULATIVE  SUMS   OF  SAMPLE   SIZES. 

NC^^\NN 

CC^O 
p  THIS  LOOP  INDEXES  OUT  ORGINAL   SAMPLE  VALUES  FOR  FURTHER  CALCULATIONS. 
Ll:CC^CC+l 
a         LOCATE   CC'TH   ORIGINAL   SAMPLE   VALUES   IN   3 

Xl-c5C(iVCCCC]+\iVCCC]  )]  ^  .„ 

P         RECORD  POSITIONS  OF  ELEMENTS   OF  X   VECTOR  IN  B  AND  ASSIGN  TO   C 

C-^DD\X1 

p  THIS  LOOP  DETECTS  TIED  INDEXED  POSITIONS  AND  INCREMENTS  TEE  INDEXING  OF 

EACH   SUCCESSIVE   TIED   POSITION   37  ONE 
p  TEST   DTH   ELEMENT   OF   C  AGAINST  REST   OF   C  FOR  TIES 

L2'.E*(.CID1=D-¥C) 
a.  SET  F  EQUAL   TO   THE  APPROPRIATE  SIZE    (DTH  SIZE)    ZERO   VECTOR 

p  CONCATENATE  F  AND   VECTOR  OF  O'S  AND    I'S  (I'S  APPEAR  WHEN  TIE  OCCURED) 

F^F  ,E 
p  ADD  RESULTANT  F  VECTOR  TO   C 

C-^C*F 

D-f-D+l 
R  CONTINUE   FOR   ENTIRE   C   VECTOR 

^(D^(pC))/L2 
XXLCC\1*C, (BB-NiCCl )pO 
*iCC<AA)/tl 
V 


7  IN-^INPUT  A',B'.X;7 

111      fl  THIS  FUNCTION  DOES  MOST  OF  THE  INPUT  PROMPTING  AND  ERROR  CHECKING. 
2    p  A=TBE  TIPE  OF  PROMPT  DESIRED. 

.3,  fi   IT  IS  A   SUBPROGRAM  OF  SIGN,    WISIG,    MANW,    KEN,    SPMAN ,   AND  NPLR. 
>.    •*(L1.L2,L3,LU,L5)CA] 

5  LI -.'ENTER  THE  DATA    (MORE  THAN  TWO  OBSERVATIONS  ARE  REQUIRED).    ' 

6  JiV-oO 

7  ->-((pQlN)  =  q)/El 

8  *l(pIN)=2)/El 
9]    *0 

.10]  L2: 'ENTER  X  DATA    (MORE  TEAN  TWO  OBSERVATIONS  ARE  REQUIRED).* 
11    X-«" 


,12    ^appX)  =  0)/El 

.lul  'ENTER  I  DATA    (NUMBER  OF  I  ENTRIES  MUST  EQUAL  NUMBER  OF  X  ENTRIES).' 

15]  i^q 

.15.   -►(TppJ)  =  0)/£'2 
.17.   *((pX)*(pr))/E2 
.18,  IN-^(pXX;X,I 

;2  0;  L3: 'ENTER  TEE  EIPOTHESIZED  MEDIAN.* 

21     INlrQ 

,22,  -^((pIN)>l)/EZ 

,23,   -►O 

.2u1  Ln: 'ENTER  THE  EIPOTHESIZED  MEDIAN  FOR  TEE  DIFFERENCES  OF  TEE  PAIRED  DATA.* 

,25]  IN->ra 

,26    ■*((pIN)>l)/E3 

'^:     '0 

;28l  L5:'  ENTER   THE   DESIRED   CONFIDENCE   COEFFICIENT:     ' 

,23.   •  • 

,30.   '  FOR  EXAMPLE:    ENTER   95,  FOR  A    95   PERCENT   CONFIDENCE  INTERVAL.' 

31    IN->-0 

.32.  *((,IN^0)^(IN>100))/En 
.33.  *((1\InUq)/E5 
.3U.   ->-0 

.3  5,  El: 'ERROR:    TEE  SIZE  OF  JOUR   SAMPLE  IS  LESS  TEAN  TERES;    TRI  AGAIN.* 
,36,   •  • 
.37.  *(A  =  1)/L1 
,38,   *L2 

.3  9.  E2:< ERROR:    SAMPLE  SIZES  ARE  NOT  EQUAL:    WANT  TO   TRI  AGAIN?    (J/N).< 

4  0    B--0 
>i;  *(B='I'  )/L2 
Cu2J  'ENTER  RIGHT  ARROW      -      TO   QUIT.* 
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E3:'^ ERROR:    TEE  HIPOTEESIZED   MEDIAN  MUST  BE  A   SINGLE   VALUE;    TRI  AGAIN.' 
*(A=3)/L3 

*L'4 

E'-i:'^ ERROR:    THIS   VALUE  MUST  LIE  BETWEEN   0  AND    100;  TRI  AGAIN.' 

E5:<^ERR0R:   TBIS  VALUE  MUST  BE  AN  INTEGER;    TRI  AGAIN.' 

*Ls 
7 

7  INTERt-A   INTER?  B;C:D:E;F:G:FF:GG;PL 
fl      TEIS   FUNCTION  IS  A   SUBPROGRAM  OF  THE  KRUSKAL-WALLIS    (KRWL) , 
p  KENDALL'S   B    UEN),    AND   SPEARMAN'S  R    (SPMAN  AND   SPMANl) .    WHEN 
fi  PASSED   THE  TEST  STATISTIC  AS   THE  LEFT  ARGUMENT,    IT   CALCULATES  THE 
a   ASSOCIATED   P-7ALUE   BY   INTERPOLATION  OF  THE   TABLE   VALUES   OF  STATS. 
a    AND    CDF'S   WHICH   ARE   PASSED   AS   THE   RIGHT  ARGUMENT. 

a      SEPARATE  THE  CDF  TABLE  AND   STATS.    INTO   SINGLE   VARIABLES. 

^F-irB  ''!•' 

GG-<raZ2-A 
n     FIND   WHERE  A   FIRST  EXCEEDS  OR  EQUALS  ONE  OF  THE  TABLE   VALUES. 

C'*-<\(A^FF) 
n    IF  THE  STAT.    DOES  NOT  EQUAL  AN!  TABLE   VALUES  SET  INTER   -    "1. 

♦((+/C)=0)/L2 
fl       INDEX   LOCATION  OF  FIRST  OCCURENCE  OF  MATCH 

Df^T.1 

_.,  fl         IF  INDEXED   POSITION  EQUALS  ONE  INDEX  P-VALUE  OUT  OF  GG. 

i7j  ♦(j>i)/i:i 

is:  INTER-^GGLDI 
19]   *0 

B         OTHERWISE   CONDUCT  INTERPOLATION  TO   GET  PROPER  P-VALUE. 
LltF^FFCD-lj-FfCP] 
F-^FFLD-lZ-A 
G*GG[.D1-GGiD-ll 
PL*{FxG)iE 


20 
21 
22 
23 
24 


,  2  5  T  INTER-rGG  [D  - 1  ]  +PL  . 
l25j   -*0 

[27]  L2 : INTER*   1 
7 


7  KENP^KENDALP  NiA',B;C:D-.EiNN:X;XX\AA;F;TiCiP:TPL 
B        THIS   FUNCTION  IS   A   SUBPrOGRXm  OF  KENDALL'S   B    (.KEN)   AND  NON- 
fl  PARAMETRIC  LINEAR  REGRESSION    (.NPLR) .      IT  CALCULATES  THE   CUMULATIVE  "^ 

B  DISTRIBUTION  OF  3   FOR  A   SAMPLE  SIZE  N. 

X-<-  1  1 
fl     INITIALIZE  FREQUENCIES  FOR  X  FOR  SAMPLE  SIZE  NN. 

NN*2 
p        DETERMINE  SIZE  OF  RIGHT  PROB .    TAIL   VECTOR 

G*(lUNy-(N-l))tn))-i-i 
p  OUTER   LOOP  INCREMENTS  THROUGH  THE  N  SAMPLE  SIZES  TILL   THE  DESIRED   ONE  IS 

GENERATED 
Ll:Dt-pQ 
CfO 
F-oX 
NN->-NN+l 

B-UNNx^NN-l))*!)*! 
A*B 
B  INNER  LOOP  GENERATES  NN+1   FREQUENCIES  FROM  THE  VECTOR  OF  NN  FREQUENCIES 
L2:A*-A-1 
C^C+XCCS-A)] 
Di-D.C 
AAt-Q 
p      IF  SIZE  OF  D   EQUALS  NN  AND   INDEXES  OF  X   STILL  REMAIN   GO   TO   Ln 

*(.(.(pD)=NN)A((B-A)<F))/L^ 
fl        OTHERWISE   CONTINUE  TO   INCREMENT  THRU  OLD   FREQS 

*((B-A)<F)/L2 
p       WHEN   LEFT  HALF  OF  NN+1    VECTOR  IS   COMPLETE   GO   TO   L5 
■*L5 
,^.,  a       THIS  LOOP  ALLOWS   ONLI   NN  TERMS  TO   3E   USED   TO   GENERATE  NEW   VECTOR 
;283  L-'-.A^A-l 
.29"  AA-^AA^l 
.30l  C*C+XL(.B-A)1-XLAA1 
^31.  D*D,C 
32    *((L((S-1)*2)  +  1)>(5->1))/Z;u 

33.  p  INVERT   VECTOR  D  AND  ASSIGN  TO   E 

.3u]  L5:E-*-eD 

;3  5]  fl  FOR  NN  OF  APPROPRIATE  SIZE  EITHER  DROP  FIRST   VALUE  OFF  E  OR  NOT 
^36.   ♦(A?W=3,6,7,10.11,m,15,18,19)/L3 

;33;  fl  COMPLETE  NEW  VECTOR  OF  FREQS  X  Bl  CONCATENATING  D  WITH  E 

•iO;  fl  *    '     CONTINUE   UNTIL   SIZE  OF  SAMPLE  N  IS  REACHED 

•+1,  *(.NN<N)/L1 

4  2]  fl     GENERATE   VECTOR  OF   CORRESPONDING  P  STATS  OF  PROPER  SIZE  G 
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P-t-O.iCC-l) 
s      CALCULATE  B   STATS   FROM  TEE  P   VECTOR 

T-^l  (((uxP)^(iVx  (//-I )))-!) 
p     TAKE  ONLl   G   ENTRIES   FROM  THE  FREQUENCY   VECTOR 

TPL-^Tx(N}((.N-l))i2 

fl   "  CBANGE  FREQS  TO   CDF   VALUES  AND   OUTPUT  B   STATS  W/ APPRO.    CDF   VALUES 
KENP-^(.2,C)pT,iT  +  \XX)*(.lN)).TPL 
V 


7  MAN-^N  MANWP  M;FiQ\P',S',T',U xV iB'.NN-.N-.MiUU-.MM 
B        THIS   FUNCTION  IS   A    SUBPROGRAM   OF  THE  MANN-WHITNEI  TEST    {MANW). 
fi  IT  GENERATES  THE   CUMULATIVE   DIST.    FOR   THE   U  STATISTIC.    N   =   SIZE  OF 
fl  LARGER   SAMPLE;    M   =   SIZE  OF  OTHER  SAMPLE. 

aCOMPUTE  NUMBER  OF  TERMS  TO  BE  INCLUDED   IN  LEFT  TAIL  DISTRIBUTION  LESS   1 

MM-^iliiN^MJ^Z))*! 
a  SET   F   VECTOR   EQUAL   TO    1  CONCATENATED   WITH      MM   ZERO'S 

Fi-l,MMpO 
q  SET  ?   EQUAL   TO   THE  MINIMUM   OF  N+M  OR  MM 

P-^i/dN+M)  .MM) 
fl  SET   Q   EQUAL   TO   THE  MINIMUM  OF  M  OR  MM 

Q*l/M,MM  ^  ,-. 

fi     GO   TO   LINE  DENOM  IF  MM  IS   LESS   THAN  N+liSIZE  OF   X+l') 

■*DENOM>-\{MM<N+\) 
p  IF  MM^N+1    GENERATE  FIRST  BLOCK  OF  RECURSIVE  RESULTS   USING  NUM  LOOP 
NUM'.T-^-N+l 
L2:U->r.m-^l 
n  PRIMARY   FORMULA    USED   IN   GENERATION  OF  FIRST  BLOCK  OF  RECURSIVE  RESULTS 

a    assigns' NEW   DECREMENTED    VALUE  TO   U  AND   TESTS  IF  T   <   THIS  NEW   U 

->-L1^\IT<U-^U-1) 

*L2xi  (psr-f-r+i ; 
p  GENERATE  FINAL   RECURSIVE  RESULTS   USING  DENOM  LOOP 

DENOM:  S-f-l 

p     PRIMARY   FORMULA    USED   IN  GENERATION  OF  FINAL  RECURSIVE  RESULTS 

L3:FZV^.-^F''V2  ^F''7-S1 
->-L3^\  (  (MM+l)S7-«-7+l) 

p        CONVERT  FREQUENCY  TABLE  TO   CDF   VALUES  FOR  FINAL   OUTPUT 
MAN-*-(.+\F)*(.N'.(.N+M)) 
7 


7 
8 

9] 
10 
11 
12 
13 


V    Z-^NORMCDF  XiA-.B-.C-.D 
a    EVALUATES  NORMAL   CDF  AT   VECTOR   X.    FOR   IX<4 ,    26.2.11    IN  ABRAMOWITZ   AND 
p    STEUU7r,~P .~^^2T~I3~USED  .'FUR'ZARGE  XT    THE   CONTINUED   TRACTJVTrJirWAZZ~    ■ 
PP. "3377    9!2.11,    ls~VSEV-IT  DFPTH-IST   ApPE7iRS~W~m7E~AT-ZE7!.ST~l^-SIVN7- 
aFICANT  FIGURES  .'PORTED   TV  miNFRAME   lyi^T'IUW  I3J~ff3S   BFTETT  CHAN'GED'TO 
p"37DTD   VNDERTZOW  PROBLEMS  WITH   x\ 

->-((p.4-«-(|X<4)/Z-«-X-«-,X)  =  p,X)/3+DLC 

S-f-nxSU)/X 

->-(.ip,B)  =  pX)/BIG  ^  ,  n 

ilf-0,5+(^i4x((*A*2)xo2)*0.5)x+/x\CA*2)o.*    l+2x  i  (C-«-10r  nOxf/ I^  )+l 

ZC(    X<u)/ipX3-^il 

-►ap,>l)  =  pX)/0 
BIGlC-*-    16589790  56295540  52050600  1993U6U0  3680160  3U1952  15232  256 

£1+  2027025  32U32400  75675600  60540480  21621600  3843840  349440  15360  256 

fl-<-l-(fl  +  2x((o2)x*fl*2)*0.5)x(+/((0.5x,S*2)o.*0,i7)x((p,S),8)pC)  ++/((0.5x,5* 


^   ^   2)0  *0,i8)x((p  S),g)pO 
[15]   ZC(U24)/ipX]-^5 


7    Zf-NORMPTH  P-.A-.B tCjD-.qiT-.S'.RiF 
p    IMPLEMENTS  ALGORITHM  AS   111    BY   BEASLEY        SPRINGER,    APPLIED   STAT,    1977 
p    FUR   3   7TCT0R~ITrPUT~Z}F  FRACTIONS , 'RETURNS   CORRESPO NDIWT NORM AZ"0 (JAN - 
nTIZES~CZIIMFD   TlZrCVRATTI  TS~BETTER  TB'AT}-^-l>x'i-^*-^--¥DR'T)RE]i.Ttm"ACtUR'ACY , 
<^ESPYCi:^ZZ7~  :'DR   'XTP.'EFiE   ?    '/kZITEir^   aDZ'ONE   OR   tlORE  NEJ7TDT!~?.2:?'HSDN~ZuDPS. 

-~^/T?23  jT(?'5I}"J7?32 

-(v/c  .  \Q^,P-Q.i )50.42) )/3+QLC 

S*-Z^,Q 

■*EXT 

r-"- (0.422  I  0)/Z-^,Q 

5-^(0. 42<|«)/i 

A-^    2.50662823884   18.61500062  52  41.39119773  534 

S+   8.4735109309  23.08336743743 


25.44106049637 
21.06224101826  3.13082909833 


r-»-rx(  ((2'*2)».*o,i3)+.xA)  +  i+(  (r*2  )o.*i4)+.xfl 
zr(o.422|g)/ip,ii]*r 

♦(F=1)Z0 
EXT-.C-^      2.78718931138    "2.29796479134    4.85014127135-2.32121276858 
D*-    3.54388924752    1.63706781897 
S*(xS)x((i?o.*o,i3  )  +  .xC)*l  +  (((i?*(l«0.5-|S)*0.5)<».*    1    2)+.xZ)) 
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l\i 


ZCC0.u2<l(3J/ipC3*S 
♦0 
ERR: 'ONE  OR  MORE  P   VALUES  ARE  OUT  OF  RANGE. 


V   PirPERVl  N',X;Y;Z 
fl  TEIS   FUNCTION  IS   CALLED   31  SPEARP.    IT   GENERATES  ALL   POSSIBLE 
a    PERMUTATIONS   OF  N  RANKS.    N= SAMPLE  SIZE. 
■*-0>^  \N~P-^    1  1  pi 
/■i]    Z->-PERM   N-1 
■53    P-«-iX*0 

7I    !-*-(-(  iiV)eX)\Z 
/si    IZiXl-'-N 

9]    P-'-UXxlJV-l),iV)p(,P),.7 
:iO]  *L1 

7 


7  M-'-N   QUANC   Q;I ;J :K:L;M:U 

[1]  P  TBIS   FUNCTION   GIVES  A    CHOICE  OF  NONPARAMETRIC   CONFIDENCE  INTERVALS 

2:  0  FOR   THE   QIH   QUANTILE  OF  A    CONTINUOUS   POPULATION.      N= SAMPLE  SIZE, 

3]  B  Q=QUANTILE.    IT   CALLS  THE  SUBPROGRAM   CBIN . 

Is]  a  CENTERS   CHOICES  AT  -THE  -ORDER   STATISTIC  NEAREST  ESTIMATE . 

LOJ     I->-'.0  .  z-rNxQ^-_20 

rr]  a    J   IS    I. z    STANDERD   DEVIATIONS {APPROX . ) 

CSj     J-«-L0.  5-rl.  3x  (Jxl-iJ-rlOC  )*0  .  5 
l9J      I^I'r.Jx      \      -■_ 

;10j  fl  RUN   OUT  K  BOTH   WAIS. 

fill   Ji:-^  l  +  xl  +  L/(IC2]  ,/V-JCl]  ) 

[12 J  fl  CONFIDENCE   COEFFICIENTS 

[I3j       L^-/<^{2.qK)o{N,Q.0\^Q)    CBIN  ((.liXli+K-DAI  111 -K)) 

r  iu]      M^ 1  ORDER   STATISTICS  I  COEFFICIENTS ' 

LIS]  a    P^ijr   c??   FORMATTING  OUTPUT 

[15:   3J*.w,;0.1]  290  (  (ISO '-'),'''  ) 

'::•:  7-  *  :  6  c  59(2.pa)o (i-r( ti':2j-x),i-ci]+ii:-i) ) 

ris]   £/*£7,((d/?).4)o'    tl       .   > 
Cl93  M-U, ill    U,    12  5  »((p^),l)p£ 


7  SPW*iV  SPAPROX  X-J 
fl    THIS  FUNCTION  IS  A   SUBPROGRAM  OF  SPEARMAN'S  R    (SPMANP) 
fl  IT  APPROXIMATES   THE   CUMMULATIVE  PROB   FOR   R   WHEN  PASSED   THE   SAMPLE 
a    SIZE  IN  THE  LEFT  ARGUMENT  AND   THE  ABSOLUTE   VALUE  OF  R  IN  THE  RIGHT 
p  ARGUMENT.      SUBPROGRAMS   OF  THIS  FUNCTION  INCLUDE:    TDISTN 

a    CALCULATE  TEE   CONTINUITY   CORRECTION 

J-<-5tNx~1+/V*2 
p  TRANSFORM  THE  STATISTIC  R  INTO   ONE  THAT  CAN  BE   USED   WITH  THE  STUDENT 
P  T   DISTRIBUTION 

X-<-(X-J)x((iV-2)  +  l-(X-7)*2)*0.5 
p  CALL   THE  T  DIST   FUNCTION  TO   CALCULATE  THE  P -VALUE 

SPN*l-iN-2)   TDISTN  X 

7 


7  SPEAR-rSPEARP  N',Cl;A',B]CiD:E:M:N',Dl',D2\LIM-,R:CDF 

;i]  P    THIS   FUNCTION   IS  A    sbBPROGRAM   OF   SPEARMAN'S  R    (SPMANP).    IT 

2 J  R  CALCULATES   THE  EXACT   CUMULATIVE  DIST.    FOR  R   FOR   THE   SAMPLE  SIZE 

3]  p  PASSED   AS   THE  RIGHT  ARGUMENT.    BECAUSE  OF  THE  LARGE   COMPUTER   MEMORY 

>]  p  REQUIREMENTS.    N   IS   LIMITED   TO   SIX   ON  THE  PC  AND    7    ON  THE  MAINFRAME. 

53  P  SUBPROGRAMS   CALLED   BI  THIS  FUNCTION  INCLUDE:    PERM 

5J    A*fl-«-0 

7]  p  INITIALIZE   VARIABLES.    VECTORS,    AND  MATRICES. 
^91    Cl-<-Dl— 30 
^3  J    M-«-(iV,.7)oO 
UC:  C^O.-AN-l) 

11'  s.    T?:iS   LOOP    GENERATES   AN   .VxA7  ARRA7    OF   THE   POSSIBLE    VALUES   OF   DIFFERENCES 

;i2  3  0  BETWEEN  ANY  TWO  PAIRED  RANKS   BETWEEN  SAMPLES. 

,13,  Ll:S+3+l 

11*;  MLB;1'>-C-A 

15    A*A+1 

15    *(B<N)/L1 

17  p  NOW  CALCULATE  THE  SQUARES  OF  ALL  POSSIBLE  DIFFERENCES. 

13    Mi-M*2 

19  P  CALL   PERM  TO   LIST  ALL  POSSIBLE  PERMUTATIONS   OF  N  NUMBERS 

20,  D-rPERM  N 

21  p  CALCULATE  SIZE  LIMIT  OF  FINAL   VECTOR   OF  R  STATS., 
J22    LIM'>-1-^(({N*3)-N)*12) 

23]  p  INITIALIZE   VALUES   BEFORE  INDEXING  OUT  COMBINATIONS   OF  ALL  POSSIBLE 
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;2U]    fl    SQUARED    VALUES. 


'.291      EtA;:\'^MU;DL:A21 


30 
31 
32 

33] 

■311,- 

35 
36 
37 
33] 

^h 

40. 
Ul. 


i?-<-0 

fl  raJS  LOOP   CALCULATES  ALL   POSSIBLE   COMBINATIONS   OF  THE   SQUARED    VALUES. 
L2:A->-A  +  l 


lA<N)/L2 
«  ADD   DOWN  ALL   ROWS  FOR  EACH   COLUMN  TO   SUM  UP   SQUARES   COMBINATIONS. 

D2-^+/E 
a   ADD   UP   NUMBER   OF  DUPLICATED   SUMS   OF   D- SQUARED    VALUES  AND   COMPRESS 
fl  VECTOR   DOWN   TO    UNIQUE    VALUES. 


L3:Cl-«-Cl,  (L/1'2  ) 
Dl^Dl.l  +  /l(.l/D  _. 
D2'<-i(l/D2)xD2)/D2 


Dl->-Dl,(+/l(l/D2)=D2)) 


ri3 

15 

18 
19 
20 
21 
22 
23 
2H 

.26 

■^27 

28 

29 

30 
31 
32 
33 
3U 
35 
36 
37 
38 
39 
UO 

m 

U2 
43 
44 
45 
46 
U7 


'40  j 

■  30' 

C51 

[52 


■((pCl)<LIM)/L3 
a    TRANSFORM   SUM  OF  SQUARES   VALUES   TO   SPEARMAN'S  R   STATISTIC. 

R->-l-(.6xcl)i{Nx(.(N*2)-l)) 
p  CALCULATE   CDF   VALUES   ASSOCIATED   WITB   THE  R   STATISTIC. 

CDF*{->r\Dl  )r  WV 

fl  FORM  TWO   ROW   MATRIX   FOR   OUTPUT   OF  R   STATS  AND   CDF   VALUES. 
SPEAR* ( 2  ,  (pCl  )  )p  (i? , CDF ) 
7 


7  SPMt-X  SPMANP  7;C:D;DD',DliD2;N;DEN0MR;XX;Il;NS',NUMR;P;PViPVALiSU;SV;RH0; 
U ;V ;ARHO ; Ul ;V1 iX;Y iWW 
[1]   fl       THIS   FUNCTION   Is   A    SUBPROGRAM   OF  NONPAR   LINEAR   REGRESSION    (NPLR) 

l2J  a    AND   SPEARMAN'S   R    (SPMAN).       IT   COMPUTES   THE   SPEARMAN   R   STATISTIC 
[3]   a  AND   ASSOCIATED   P-VALUES .    THE   LEFT  ARGUMENT  THAT   IS   PASSED   IS   THE   X 
M       fl  SAMPLE;    THE  RIGHT  ARGUMENT  IS   THE   I   SAMPLE. 
C5J   a  SUBPROGRAMS   OF   THIS   FUNCTION   INCLUDE:    TIES,    TIESK,    SPEARP .    SPAPROX, 

5]   fl  INTERP,    AND   THE   VARIABLE  PMATSP . 

7] 

8]   fl  ORDER   I  IN  INCREASING  ORDER   OF  X 

9]    V-*-IZUl 
"a  ORDER   X  IN  INCREASING   ORDER 

P  COMPUTE   CURRENT  RANKING   OF   I 

C-lrltAV 

a  NOW  ORDER   I  RANKS  IN  INCREASING  ORDER 

Dl-*-VUVl 
P      IF  TIES   EXIST  IN  EITHER  X  OR   I  RANKED   VECTOR   USE  MID -RANK  METHOD 

DD-^1   TIES   Dl 

XX-t-1   TIES   U 
P         FIND   ORIGINAL  RANKING  OF  I  WITB  TIES  RESOLVED 

Il-^DDLCl 
p  RECORD   SIZE  OF  INPUT   VECTOR 

N-<rpX 
p  CALCULATE  DIFFERENCES  BETWEEN  RANKS  OF  X  AND   7  VECTORS 

D'>-XX-Y1 
p  DETERMINE  TEE  SUN  OF  SQUARES   OF  TEE  DIFFERENCES 

Z32-<-+/(D*2  ) 
p     OBTAIN  THE  NUMBER  OF  TIES  IN  EACH   VECTOR   USING  TEE  TIESK  FUNCTION 

U\-rTIESK   U 

Vl-t-TIESK  Dl 


SU-<-((+/(C/l*3))-(+/£/l))  +  12 
SK-<- ((+/( 71*3  ))-(+/7l)  5  +  1: 


_   .  ,  .  . ^2 

NS*N^(.(N*2)-l) 

p  CALCULATE  THE  R   STATISTIC  INCLUDING  THE   CORRECTION  FOR  TIES 
NUMRir(NS)  +  (r6)xD2)  +  i(.    6  )^(SU+SV)) 
DENOMRir{  (A?S-  (I2xsa$  )*0.  5  )x  (  (iVS-(12xS7)  )*0  .  5  ) 

RHO-^NUMRiDENOMR 
ARHO->-\RHO 

p         CALL  SPEARP  TO  CALCULATE  THE  RIGHT  TAIL  OF  THE  CDF  OF  R 

PirSPEARP  N 

*L2 
Ll:^(N>10)/L3 

P-f-PMATSPLCN-S)-,-,! 
p  CHANGE  SIZE  OF  P  TO   AN  M^N  MATRIX 

P-^PCl;;] 
p    CALL   INTERP  TO   CALCULATE  P-VALUE   BI  INTERPOLATION 
L2:PVALfARH0    INTERP   ? 


■'433      ^iPVAL='l,/L7' 


PVAL^Q 


p    CALCULATE  P   VALUE   USING  STUDENT  T  APPROX. 
...    L3:PVAL-^N   SPAPROX  ARHO 
53]    L7:SPMir(RH0),PVAL 
7 
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7  Pf-K  TDISTN  X;7 
"  IM^^l^^ZTM^J^  ^  ^^^SP-'^OGRAM  OF  THE   CUMULATIVE  PROBABILITY   GENERATOR 
n  FOR   SPEARMAN'S  R    iSPEARP).      IT   CALCULATES   THE   CDF  AT  X   USING 
p  THE  STUDENT'S   T  DIST  WITH   K  DEGREES   OF   FREEDOM. 
^7^\xi  If'^^-''^^^   CALLS   ON  THE      <F'    DISTRIBUTION  FUNCTION    iFDISTN). 

P-«-0.5i(l,X)  FDISTN  X*2 

PC7/ipX3-«-0.5+7/P 
*0xii=A/7 

pc(~/)/ip;^:*o.5-(~v)/? 


7  TI--3B   TIES  BiC;DiI;N;T',I:Z;K:M;L;PP;NRiMM 
ri]   R     THIS  FUNCTION  IS  A   SUBPROGRAM  OF  KENDALL'S   B    (KEN).    SPEARMAN'S 
•,2]  a    R    :SP>yiNP).    KRUSKAL-WALLIS    (KRWL)  .    MANN-WHITNEI    (MANW) 
'3]   a  .iND    WILCOXON    ^WISIG) .       IT   CHECKS   THE   RIGHT  ARG .    VECTOR    FOR   TIES   AND 
[ul   n  CHANGES   THE  TIED  POSITIONS   OF  THE  LEFT  ARG.    BI   THE  MIDRANK  METHOD, 

Is:      .Vf-cs 

7]  fl  IF  NO    VECTOR  OF  RANKS  IS  PASSED;    GENERATE   ONE 

8]    *((poSfl)*0)/Z,5 
9]    BB'>-\N 

10]  L6:J-<-pO 

11  L-<-A?pO 

12  T-*-l 

.13]  q  CHECKING  FOR   TIES   31  INCREMENTING  THRU   THE   VECTOR 

ll^\    L3:C-*-(.T*3)=3rT: 

Cl5j  0  COUNT  NUMBER   OF  TIES;    IF  NO   TIES   GO   TO   L2 

Cl7J  *(D=0)/L2 

ClS]  0  RECORD   WHERE  TIES   STARTED  AND   HOW  MAN!  RANKS  INVOLVED 
[191  I*I,TAD  +  1) 

[201  a  INCREMENT  NEXT  T  31  THE  NUMBER   OF  TIES  ENCOUNTERED  PLUS   1 

[21]  L2:T^T*-(D-^l) 

C22I  qJF  T  LESS   THAN  SIZE  OF  ORIGINAL   VECTOR   GO   TO   L3    AND   START  AGAIN  AT  NEW  T 
23,  •'(.T<N)/L3 
I*zl 
fl  ASSIGN  TEE  RANKS   OF  TEE  LEFT  ARG.    TO   TI 

TIi-BB 
fl  IF  NO   TIES  FOUND   QUIT 

*(J=0)/0 

p  LOCATE  THE  INDEXED  POSITIONS   OF  TIED  RANKS 
L5:PP-^((I[1  +  Z]  )-l)  +  i((JCl  +  Z]+J[2+Z]  )-l) 
P  FIND   THE  MIDRANK   VALUE  OF  THESE  RANKS 

NR^(+/TIIPP1  )+JC2+Z] 
p  SET  UP  VECTOR  WITH  ZEROS  AND  ONES;    ONES  WHERE  TIE  RANKS  INVOLVED 

K-^0 

MM-t-NoO 
L^:K*K+1 

M'<-iTI=TILPPZK12') 

MM-'-MM+M 

*(X<IC2+Z3 )/Lu 
fl  SET   UP  A    VECTOR  WITE   ZEROS  WEERE  TIED  RANKS   OCCUR 

L"-~MM 
p  TRANSFORM  ONES   OF  MM   VECTOR  TO  MIDRANK   VALUE 

MM-'-NRxMM 
p  TRANSFORM   ONES   OF  L   VECTOR  TO  REMAINING   UNCEANGED  RANKING   VECTOR 

L'*-TI*L 
n  FILL   IN  MIDRANK   VALUES 

TI->-L+MM 

Z-<-Z  +  2 
p    DO   TEE  SAME  FOR  ANI   OTHER  TIES  INVOLVED   BUT  WITE  NEWLY   COMPUTED  TI 

*(iY>Z)/L5 

7 


2u 

^26 
27 
28 

;29 

30 
31 
*^32 
33 
34 
35 
36 
37 
33 
39 

;uo 

.U2 
^U3 
nu 
U5 
46 
47 
U8 
49 

.51 


7  TIE-^TIESK  AAxAA',B;C;D;I-.NiT 

Li:   5      THIS   FUNCTION   .S    ".  CUBPROGRAM   OF  KENDALL'S   3  '.ZEN),    SPEARMAN'S 
[22   a  .?  {SPMAN)    AND    -.SPMANl  )  .    AND    KRUSKAL-WALLIS    kKRWL).       IT   CHECKS    THE 
■,3]   a  P.IGHT   ARGUMENT   FOR    TIES   AND    RECORDS   THE   NUMBER    OF   OCCURENCES   OF   EACH 
'       a  TIE  AND   THE   TOTAL   NUMBER   OF  TIES   IN   THE   VECTOR. 


p      ASSIGN  ORDERED      VECTOR     TO  B  AND  INITIALIZE  VALUES 

NfoAA 

B*-AALMA1 
9]    TIEi-pO 
10]  Tt-1 
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:ilj  0  CHECKING  FOR   TIES   31  INCREMENTING  THRU  THE   VECTOR 

"12]  L3:C-<riT*B)=BCT2 

13]  fl  COUNT  NUMBER  OF  TIES'.    IF  NO  TIES  GO  TO  L2 

,iu1  D-^+/C 
•^15]   •^(iJ  =  C)/L2 

16J  TIE*TIE ,  (.D  +  1) 

17]  a  INCREMENT  NEXT  T  BI  TEE  NUMBER   OF  TIES  ENCOUNTERED  PLUS   1 

18J  L2:T-^T-*-(D  +  l) 

J 19]  RIF  T  LESS  THAN  SIZE  OF  ORIGINAL   VECTOR  GO  TO  L3  AND  START  AGAIN  AT  NEW  T 
:20]   *(r<(iV-l))/L3 
7 


7  VAR-^VARVm  B:C:D',E;D1;S1 
11]   B       THIS   FUNCTION   IS   A   SUBPROGRAM  OF  THE  MANN -WHITNEY   TEST    (.MANW). 
2]  0   IT  GENERATES   THE  RANKING   SCHEME   USED   IN   CALCULATING  THE 
3]   0  DIFFERENCES   IN   SCALE    (1  ASSIGNED   SMALLEST,    2   ASSIGNED   LARGEST.    3 

u]   a  ASSIGNED   NEXT   LARGEST.    ^  SECOND   SMALLEST .ETC .    3Z   TWOS   TILL   PROPER 
5 J   PI  SAI^LE   SIZE   IS   REACHED).    SAMPLE  SIZE  IS  PASSED   IN  THE  RIGHT  ARC. 

7J    D^O 

8  J  a  FIND   FLOOR   OF  MIDPOINT  OF   VECTOR  AND   ASSIGN   TO   C 

9]    C-^l(Bi2) 

10]  p  LOOPS  GENERATE  RANKING  VALUES  LEFT  HALF  FIRST 
11   L2:Z)-eU+l 
,12,  E-^E,D 
13    ^(ioE)=C)/L3 
,1«-,  D-^D  +  2 
!  1 5  J  S-^S ,  D 
ISj   *((p£')<C:)/L2 

:i7_  q  NOW   GENERATE  RIGHT  HALF 

L2:D1-1 
D1*D1+1 
Sli-El.Dl 
*i(pEl)=CWL5 
L6:DI*-D1+1 
El^Sl.Dl 
*(.(oEi)  =  C)/L5 
Jl-s-Dl  +  3 
El-^E1,D1 
■>-UpEl)<C)/Le 
fl  IF   SIZE  OF   VECTOR  IS -ODD   VALUE  CONCATENATE  MIDDLE  RANK  IN  BETWEEN  HALFS 
L5:-*((2Tfl)*0)/L7 

fl  IF  SIZE  IS  EVEN  CONCATENATE  LEFT  HALF  WITH  THE  REVERSE  OF  THE  RIGHT 
VAR-^-E ,  (QEl) 
->-Q 
L7:VAR*E.B,(&F1) 
7 


7  WIL^WILP  NNiN',A;PiTiNN:WiPP',NM 

;i]  fl     THIS  FUNCTION  IS   A   SUBPROGRAM  OF  THE  WILCOXON  SIGNED  RANK  TEST 

.2,  fl  iWISIG).    IT   GENERATES   THE   CUMULATIVE  DIST.    FOR   THE  TEST  STATISTIC 

.3.  fl  THE   GENERATOR   USES  A   RECURSIVE  FORMULA,    NN= SAMPLE  SIZE. 

.4,  WM-e-(L((+/iiV/V)  +  2))  +  l 

5  Nf-2 

,6,  0         SET  P  EQUAL   TO  PHOB .    DIST.    WHEN  N  EQUALS   2. 

7  P-fUpi 

,8,  LS'.N-^N+l 

,9]  A-cl 

,10]  0  SET  T  VECTOR  TO  PROPER  SIZE  OF  ZEROS. 

,11.  r-^(+/(iiV))pO 

,12,  fl  IF  A£N   USE  TRUNCATED   RELATION  TO   COMPUTE  OCCURRENCES. 

,13,  L^'.-^U^N)/L1 

,1^.  a  IF  A>N  USE  FULL  FORMULA   TO   COMPUTE  OCCURRENCES. 

.15.  *(.A>N)/L2 

„16,  fl  WHILE    (A-N)   IS  NEGATIVE  TRUNCATE  FORMULA      TO   AVOID   A   NEGATIVE  INDEX. 

Cl7]  LliTUl-'-PU} 

■;i8]  A^A  +  1 

,19.  +2:u 

20  fl        IF  A   IS   LARGER   THAN  THE  LENGTH   OF  P   GO   TO   LB . 

l21j  L2:*{A>iaP))/La 

r22]  a     ONCE    iA-i^)    BECOMES   POSITIVE:    THE  RECURSIVE   FORMULA    CAN   3E   USED. 

'23]  TZAI^PiU-i^Dl^PCA: 

2Uj  -t-L? 

12 5"  fl      ONCE  A   IS   LARGER   THAN  THE  LENGTH   OF  P  TRUNCATE  FUNCTION  AGAIN, 

26  Z,6:rCil]-^P[(il-iV)] 

,27.  L7:Ai-A+l 

28  fl     IF  A   AS   AN   INDEX   HAS  NOT  EXCEEDED   N(.N+l)/2    CO  AGAIN. 

29  ■*(A^(+/(\N)))/L^ 

.30  fl  CONVERT  T  INTO   P  AND   CONCATENATE   1  FOR   USE  IN  NEXT  ITERATION  OR  OUTPUT. 

.31.  L5:P'^((  +  /(.\N))pT).l 

,32.  PP-<rNM^P 

,33.  WIL->-(+\PP)*(2*NN) 

'3U  p    CHECK   IF  LENGTH   OP  INPUT   VECTOR  EXCEEDS  NUMBER  OF  N'S   GENERATED. 

C3  53  ^iNN>N)/L3 
7 
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APPENDIX  G 

LISTINGS  OF  PROGRAMS  USED  TO  GENERATE  C.D.F.  COMPARISON 

TABLES 


7  KENTEST'.NiALPHAiB;A;C',TAOiP'.NUM;DENiZiZZ;ERRZZ;M'.ZZCiDiAA',PP;NUMC;I;H;F 
lFS;S;J',EkRZZC',iC:KK 

PROGRAM   GENERATES   TABLES   OF   C.D.F.    COMPARISONS   FOR   KENDALL'S   3 

SET  SAMPLE  SIZE  AND   ALP3A    VALUES. 

Ni-ll 

ALPHA*    0.010.C2J.:3C.:5  0.03a.l33.1S  i 

n  TSIS  LOOP  INCREMENTS  SAMPLE   SIZE. 

Ll:PP*pO 

Mf-pO 

XJf-«-pO 

FS-«-oO 

22-00 

2ZJ-00 

D*-Q 

fl  COMPOTE   CUMULATIVE  DIST  AND   ASSOCIATED   STATS. 

P*KENDALP   N 
p  T3IS  LOOP   CALCULATES  ALPEA    VALUES  AND   APPORXIMATIONS . 

L2: 5-^9-^1. 

C'^ALfPHrL  l3  a 

A->-(  +  /(Pt2:2^C)) 

TAU-PI'.-.At 

KK*KK,TAU 

AA->-e*N^~l-i-N*2 

PP-*-PP,PC2;A1 
p  COMPUTE  NORMAL   APPROXIMATION. 

NUMf3x(TAU)x(.(.N^(.N-l)')*0.5) 

DEN->-(.2x((,2xN)  +  5))*0.5 

Z-^NUMiDEN 

ZZ->rZZ,  (.1-NORMCDF   Z) 
p  COMPUTE  NORMAL  APPROXIMATION  WITB  CONTINUITY  CORRECTION. 

NUMC*3x(TAU-AA)^(  (.N^  (N-1 )  )*0  .  5  ) 

ZCirNUMC*DEN 

ZZC*ZZC Al-NORMCDF  ZC) 

■*iB<l)/L2 
p  COMPUTE  ERROR  DIFFERENCES . 

ERRZZ*PP-ZZ 

ERRZZC-^PP-ZZC 
B  PRINT  OUT  TABLE  OF  VALUES. 

IfO.l 

5^.  t  XEST  STAT      VALUE      ' 

«■<-«,'        *  PROBLB    t   S]  :  FOR  SAMPLE  SIZE  EQUAL   TO    ',(2  0  9N),» 

*L3 
L5:J*19n 
L3:M'>-M.tn    [U7C(18pl97),.7,61p(8pl97),J': 

■»-(£i  =  u)/L4 

J?-<-  1  18  pS 

*(Z3xO)/Z;5 

F-<-'XUC<T  9.9999  >'  DFMT(.l    7    pKK) 

F*-    1  62  pF 

FS-^E,F 

■*L1 
LB:FS*pQ 

H->-17oS 

FS*FS,3 

C-^-Q 

/-<-'  I  '  (7  5  9PPIC2  ) 
FS-'^FS.F 
*(.C<7)/L9 
L7:M-<-M,Cl]  FS 
D-^D+1 
1*1 


*<iNl,N2.N3.L5)[.D2 
NltS*^ EXACT  C.D.F. 

.7-«-198 

*L3 
N2:S*' ERROR;    NORMAL 

PP*ERRZZ 

*L3 


1  07 


30" 

31 

32 

33 

3U 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 
51 
52 
53 


N3 IS*' ERROR;    NORM.    W/CC    ' 

PP*-SRRZZC 

■*L3 
LH:M 

->(iV<14)/Ll 

7 


7    KWTE5T ;A-,B',C; ?P : PC : NN: ALPHA  \N',K',P'.F',E', PVALUE ; PVAL ; PFl ; PF ; PVF ; PVFl ; D ; iTHJ? 
fl  ;  £:i?iJ  ? ;  Mi??  1 ;  CD  F ;  CC ;  ii: ;  I ;  J :  Ff :  S ;  M 

R  THIS  PROGRAM   GENERATES  TABLES   OF  C .D .F .    COMPARISONS  FOR  THE  KRUSKAL- 

fl  WALLIS  TEST. 

i?-*-    4     4    4 
CCf-2 

ALPUA-<r    0.01    0.02    0.03    0.05    0.08    0.13    0.18 

^L2 
Ll:At-NN 

p  CALL   KRUWALP  TO    GENERATE  EXACT  DISTRIBUTION  FOR   SELECTED   SAMPLE  SIZE. 
L2:?*KRUWAZP   A 

B-^0 

N-^+/A 

CC-^CC+l 

K*-pA 

Mf-pO 

PP-f-D  0 

KK-^oO 

?C-*-oO 

PF^oO 

PPl-oO 

CDF^PL2;1 
L4:P-*-5+l 
R  LOCATE  POSITION   OF   EXACT   CDF   VALUE   5  ALPHA. 

D-^{+/(,CDF^ALPHAL31  )) 

-^(D=0)/L3 

D*-l 
a  DETERMINE   CORRESPONDING  TEST  STATISTIC   VALUE. 
L3  :£(^P[1;l)] 

KK'*-KK   H 
fi    RECOkD   EXACT   VALUE  OF   CDF. 

PVALUE^Pll-.Dl 

PPt-PP. PVALUE 
fl  COMPUTE   CORRESPONDING  P-VALUE   USING  CEISQ  APPROX. 

PVALirU-l)    CHISQ   B 

PCf-PCAl-PVAL) 
aCOMPUTE   CORRESPONDING  P-VALUE  USING      F  APPROXIMATION 

P-^((iV-X)xff)+((ii:-l)x((A'-i  )-£?)) 

PVF'>-(.(.K-1).(N-K))    FDISTN  F 

PF-^rPF,  (1-PVF) 
fl  COMPUTE   CORRESPONDING  P-VALUE   USING  F  APPROX  W/    1  LESS  D.F.    IN  DENOM 

PVFl-'-UK-DAiN-K)-!))   FDISTN   F 

PPl-ePPl,  (i-pypi) 

*(B<7)/i4 

D-f-O 

ERRH-f-PP-PC 

ERRF-^PP-PF 

ERRF1-^PP-PF% 
fl  PRINT  OUT  TABLE  OF   VALUES. 

I-*--0,l 

S-'-'TEST  STAT.    VALUE    ' 

M->-M.<  PROBIH    i    '  ,DA7C46]  , '];  FOR  A    GROUP   OF   3  SAMPLES   CONSISTING  OF   4 


4,    AND    <  AfCC).  '    0§S 
54]      ^L8 

55]  L5:J'i-<-< 

56]  LS-.M-^M.Zn  (17p'-'),J,61p(8p'-'),J' 
57]       ->-(D=5)/L9 
58]       Ff-pO 
59]       flf-iepS 
60]       F-^F,H 
61]       C--0 
,62]       *(DtO)/L10 

L63]  ::ii:c-f-(;-i-i 

[54]   -p-i   ' ' , (7  4  s^ZC^i ) 
'5  31  F*F.FF 
,66.   -v(C<7)/Lll 

FF*'  T'  .(7  5  e/fifCC]) 
F^F.FF 
■*(.C<7  )/L10 
Ll2:M-^M,ll2    F 
D-*-D  +  l 
1*1 


67 
.68 

.69 

.70 

71 

.72 

73 

74 

t76 
'77 
[78 


-(.Nl  ,N2  ,N3  ,Nn  ,L5)ZD1 
Nl:S-r^ ALPHA    VALUE 


KK^PP 
J-H'  +  i 
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79]  --L8 


80 
82 


83]  N2:  Sir 'ERROR:    F  DIST 


84 
85 
86 
87 
88 
,89 
90 


iV2  :  S-<-  •  ERROR :  CEI SQUARE ' 

KK*ERRH 
L8 


KK-rERRF 
■*-L8 
Nn:S->-' ERROR;    F  W/~l   DF' 

KK->-ERRFl 
■*-L8 

-►CCC<3  )/Ll 
7 


V   K»TESTSK:RiR2xFC \Nl^ ;ALPEA;Ii ;K;J ',W ',B -,11  ',PVAL;FF1',EE ;MM;M:X\ERRC iERRFl;F ', 

PV;B'.KK;D:PVF1;Q;FF;U'.C;I;J;S 

a  THIS  PROGRAM  GENERATES  C.D.F.  COMPARISONS  FOR  THE  KROSKAL-WALLIS  TEST 
fi  BASED  ON  3  0000  RANDOM  VECTORS  OF  ^0  RANKINGS  WHICH  SIMULATE  COMPARING 
p  5  SAMPLES   OF   8  OBSERVATIONS   EACH. 

W-f-10000 

X-^5 

ffB-«-pO 

/V-^-UO 

JV/V-s-0 
Ll:NN*-NN-^l 

5-*-8 

MM-^0 
a  GENERATE  A   RANDOM   VECTOR   OF  RANKINGS. 

X-«-u0?U0 

IJ^dO 
n  SUM  TEE  RANKS   OF  EACH   SAMPLE. 
L2:y-^+/(MM+CC+X)) 

Q^Q+8 ' 

MM-^MM-^B 

■*iQ<.»'Z)/L2 

52-^+/((JJ*2)  +  8) 
R  CALCULATE  THE   H   STATISTIC  FOR   EACH   VECTOR   OF  RANKINGS. 

H-*-((3  +  '+10)xj?2)-123 

HH-*-HH   B 

♦  ((A'A'  +  10000)sl,2,3  ,H)/m 
p  CHANGE  RANDOM   SEED   EVERY   10000  ITERATIONS. 

URL*   0  60  60  1000  I'u+DTS 
Ln:*(.NN<W)/Ll 
0  ORDER  ALL   ff'S  FROM  LARGEST  TO   SMALLEST. 

BHi-HBL'VHEl 

B-f-O 

DfO 

M*oO 

ALPHA*   0.01    0.02    0.03    0.05    0.08    0.13    0.18 

KK*pO 

PFl*oO 

PC-op  0 
L3:B*B+1 
p    DETERMINE  DESIRED   B   VALUE  CORRESPONDING  TO  ALPBA    VALUE. 

B*HH [ {NN^ALPHA [S] ) 3 

KK*KK  H 

p  compOte  corresponding  p-value  using  ceisq  approx. 

PVAL-i-U-l)    CHISQ   E 

PC*PC.  il-PVAL) 
p  COMPUTE   CORRESPONDING  P-VALUE   USING  F  APPROX  W/    1  LESS  D.F.    IN  DENOM 

P*UN-K)xB)*(,^x(.39-H)) 

PVFl*((.K-l)A(.N-K)-in   FDISTN  F 

PFl-'rPFlAl-PVFl') 

*(B<7)/L3 

ERRC--ALPHA-PC 

ERRF\*ALPBA-PF1 
p  PRINT  OUT  TABLE  OF  VALUES. 

1*0.1 

S*'TEST  STAT.    VALUE    ' 

M*>^,'         PROBZH    >    ^  ,3i27^U6]  , ':  :  BASED   ON    ',(5  0  «(/),'  GENERATED   5''S  FOR    5 
SAkPLES   OF   8  OBS.    EACH.  ' 

*LS 
LS'.J*'-' 
L8: «■<-«,  [J]  (17p'-'  ),.7,61p(8p'-'),J 

*{D=^)7L% 

F*pO 

B*16dS 

F*F,B 

C*Q 

♦(Z?*0)/L10 
IllrC+C+l 

FF*'    1' ,(7  3  9KKLC2) 

F*F.FF 
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:723   •*(C<7Vill 
73]  -L\2 

7u1  iiO:C-<-<7+l 
75]   fF-«-'  r  ,  (7  5  VKKICZ  ) 
76]  F->rF,FF 
77]   *(C<7)/L10 

78]  L12:M-^M,C1]  F 
79]   Z)-*-D  +  l 
80J   J-«-l 

[3  2 J  Nl:S*^ ALPHA    VALUE 
8  3 ]  KK^ALPEA 
8u]   J-(-'+' 
85]  ■*L9, 

86]  N2:S-^' ERROR:    CEISQOARE' 
87]  KK-^ERRC 
8  8  J   ■^'Zffi 

89]  N3:S*' ERROR;    F  W/~l    DF' 
90]  KK-^ERRFl 

[32]  L3:;^ 

£/-«-  20000  30000  1 

A;-<-y[A] 

t96J   ^U<4)/L1 
V 


r32] 

[QU] 


V  MANTEST;U:NN:MMiALPHA:3;A;C;T:P:NUMiDEN:Z;ZZ:ERRZZ;M;ZZC;D:PP-NUMC:I;Ei 
F;FS:S;J  -^TT :  f ;  ITC ;  TC ; MUMC ;  iVt/MT ;  L)£iVOM ;  DENOMC ;  ERRA  VE ;  ERRTT ;  ERRTC  :  ERRZZC ;  i 

2]   «  THIS  PROGRAM   GENERATES  TABLES   OF  C.D.F.    COMPARISONS   FOR  THE  MANN- 
3]   fi  WEITNEI   TEST. 

^5]   0  SET  SAMPLE  SIZE  AND   ALPHA    VALUES. 

re]    NN^7 
■7]    MMi-7 


[8]    ALPff^-*-  0.01  0.02  0.03  0.05  0.08  0.13  0.18 
■g]   Q  THIS   LOOP   INCREMENTS   SAMPLE  SIZE. 

.10]  Ll:PP-^pO 
*^11    W-<-pO 
12"  KKfpQ 
13.   fS-t-oO 
^14.  ZZ-t-pO 
,15.   ZZC-^pO 

ff-^pO 

TTC-i-oO 

NNf-NN+1 

MM-^MM+1 

D-f-O 

S-f-0 


£15 

^17 

18 

19 
20 

;2i 

^22 
23 
24 
[25 
,26 
,27 
,28 
,29 
.30^ 

•U^ 
32 

33 


„  „. COMPUTE  CUMULATIVE  VIST  AND  ASSOCIATED  STATS. 

P-i-NN  MANWP  MM 
?-,.!:  o  .  ^^^^   ^°°^   CALCULATES  ALPHA    VALUES  AND  APPORXIMATIONS . 

C-<rALPHALBl 

>!■«-( +/(PSC)) 

t/-*-i?-l 

KK-<rKK,U 

PP-^PP,P\.A1 
".„„.  „  ..  COMPUTE  NORMAL  APPROXIMATION. 

NUM->-U-((NN^MM)t2) 
,   D£'N-»-((iVyVxMMx(w;V+MW+l))  +  12)*0.5 
3U]  Z^NUMiDEN 
3  5]  ZZ-^ZZ  .NORMCDF   Z 

'\%\     ^r,T,y,n     f  n^r.^°A^H^rlJW^k  s^^^^°^^^^^°^    ^^^^    CONTINUITY    CORRECTION. 

37]  NUMC-t-(U+0 .5)-i(NNxMM)i2) 

3  81  ZC-^NUMC*DEN 

3  9]  ZZC-^ZZC.NORMCDF    ZC 

■.■^O]  R  ^COMPUTE  STUDENT  T  APPROXIMATION 

,41]  DENOMir{{  (.NN+MM-1)-{Z*2  )  )*{NN+MM-2  )  )*0  .  5 

'■i;2]  T-^Z-i-DENOM 

43]  ZT-^rZ' ,  (NN+MM-"^  "^    '^DISTN    "" 

"^--i  \,rr..~J^^^^'^~^   STUDENT   I  APPROXIMATION  WITH   CONTINUITI   CORRECTION. 
45  J  NUMT-^i  NUM  )-<-0  .  }: 

hu^i      ^P9M?,t'-kiL'^'^.'^'^'^ '''■'■"'  ^^"^*2  ;  U(.NN+MM-2  ) )-  ( ( (A7ra+0  .  5  )'»2  )riNN+MM-2  ) )  )*0  .  5 


.47 


!48;  TTC->-TTC,CiNN+MM-2')   TDISTN  TC) 


49 
-50 
,51]  ERRZZ'>-PP-ZZ 


,52 

,54 
55 


TC*NUMT*DENOMC 


■{B<7)/t2 
P„         COMPUTE  ERROR  DIFFERENCES. 


ERRZZC-PP-ZZC 
ERRTT*PP-TT 
ERRTC-i-PP-TTC 
ERRAVE-'-PP-  UZZ+TT)*2) 
ERRTCZC-^PP-  (  (.ZZC+TTC ) ♦  2  ) 

0  PHJwr  oc/r  fiisLE  OF  values, 

I--0.1 


1  10 


>•  DFMTil    7    pKK) 


5^ I T£ST  STAT  '    VALU^      ' 

M->-m''  '      PROBZn   ^    Ul;    FOR  SAMPLE  SIZES  N  EQUAL   TO    ',(2  0  9NN),<    AND 

M  EQUAL   TO    ',(2  0  9MM)  .'^ .  ' 

LS-.M-i-M.Lll    UAVZil8pl97),J,61p(.apl97),J2 
*(Z)  =  8)/LU 
S-<-  1  18  oS 

F->-'BG<\       ZZ'S 
F*  1  5  2  0? 
FS*S,F 
■*L7 
L6:FS-<-pO 
ff-«-17DS 

L9:C-^C+1 

F*-'    I'  ,(7  5  5PPCC]  ) 

FS-^FS.F 

-(.C<7  )/L^ 
L7:A}fM,Cl]    FS 

D-f-D  +  1 


•(iVl  ,A72  ,A'3  ,Nn,N5  ,N&  ,N7  ,L5  )  IDl 
:S->-^EXACT  C.D.F.  ' 


JVl: 

•7-^198 

->-L3 
Nl-.S-^' ERROR;    NORMAL 

FP-ERRZZ 


N2:S*'E?.R0R\ 


NORM.    W/CC 


■'\^ 


,103 

,10U 

105 

.106. 


N^:S-<-< ERROR;    T  DIST 

PP->-ERRTT 
*L3 
NS-.S"-' ERROR;    T  W/CC 

♦L3 
N6:S-^- ERROR;    AVE   T/Z 
■  PP-<rERRAVE 

*L2 
N7  IS-*-' ERROR;    AVE  TC/ZC    ' 

PP*ERRTCZC 

■*L2 
LniM 

*(NN<9')/L1 
7 


7  SIGNTEST;N;ALPBA;B;A;C;KiP;Z;ZZ;ERRZZ;M;ZZC;DiPP;I;B;F;FS;S;J;ERRZZC;ZC 
fl  TBIS  PROGRAM  GENERATES  TABLES  OF  C.D.F.    COMPARISONS  FOR  TBE  SIGN  TEST. 


N-^22 
ALPBA*   0, 


SET  SAMPLE  SIZE  AND   ALPBA    VALUES, 


COMPUTE  CUMULATIVE  DIST  AND  ASSOCIATED   STATS. 


,01  0.03  0.06  0.12  0.22  0.35  0.5 
TBIS  LOOP  INCREMENTS  SAMPLE  SIZE. 
LliPP-f-pO 
M*0  0 
KK->- pO 
ZZ-«-p0 
ZZC-t-pO 
Ni-N+1 
Df-0 
St-0 

P-^BINOM  N 

a  TBIS  LOOP  CALCULATES  ALPBA    VALUES  AND   APPORXIMATIONS . 

L2:3fB+l 

c-ALPSA rsi 
A-*-K^/  {p^c: ) 
y.^A  - 1 

pp-^ppIpca] 
n  compute  normal  approximation. 

Z-«-(^- (0  .  5x/V)  )*  (0  .  5x  (;V*0  .  5  ) ) 

ZZ-^ZZ,  (.NORMCDF   Z) 
fl         COMPUTE  NORMAL   APPROXIMATION  WITB   CONTINUITI   CORRECTION, 

ZC:*((X+0.5)-(0.5x/V)U(0.5x(iV*0.5)) 

ZZC<-ZZC ,  {NORMCDF   Zc) 

*(S<7)/i2 
0  COMPUTE  ERROR  DIFFERENCES, 

ERRZZ*PP-ZZ 

ERRZZC^PP-ZZC 
P  PRINT  OUT  TABLE  OF   VALUES. 

1*0.1 


1  1  1 


*(N1,N2.N3,L5)LD2 
NllS*' EXACT  C.D.F. 


S-*-'TEST  STAT.    VALUE      • 

M*M,'  ^PROBIK   i  Kl',    FOR  SAMPLE  SIZE  EQUAL  TO    ',(2  0  «/V),» 

*L3 
L5:J-^19n 
L3:M^M.Cn    CU7C(18pl97  )  ,J,61p  (8pl97  )  ,.7] 

fff    1    18    pS 

*(D*0)/L6 

F^'3C<\       ZZ9         >'    QFMTil    7    qKK) 

F-^    1    o2    flF 

FS*H,F 

*L7 
L6:FS'>-oO 

H*17pS  -^ 

FS->rFS.H 

Cf-0 
L9:Cf-C+l 

F*<    I ',(7    5    iPPZCl) 
,.-.       FS-^FS.F 
,57]       -^(C<7)/Lg 
58]    L7:M^M,il2    FS 
.59j       I/-<-D+l 
l60j       J-^1 

62 

63]  .7-C198 

64]  t-LS 

5  5]  N2: Si- 'ERROR;    NORMAL  • 

55]  PP-^EHi?2Z 

57]  -^3 

68]  N3:  Si-' ERROR;    NORM.    W/CC    ' 

70]   'LS 
71]  LU:M 
72]   -»-(iV<26)/Ll 
V 


7  SPViTEST-,N\ALPRA;B\A%C\AA:P',Z\ZZ;ERRZZ;M\ZZC\D;PP;I;E;F;FS;S;J;TT;T',TTC; 
TC : i&i?i?il VE ; ERRTT ; SRRTC ; ERRZZC ;ZC;PCi EHO  ; KK 

fl  THIS  PROGRAM  GENERATES  TABLES  OF  C.D.F.    COMPARISONS  FOR  SPEARMAN'S  R. 

ft  SET  SAMPLE  SIZE  AND  ALPS A   VALUES, 

N-^8 

ALPBAir   0.01    0.02    0.03    0.05    0.08    0.13    0.18 
fl  TEIS  LOOP  INCREMENTS  SAMPLE  SIZE, 

Ll:PPi-pO 

Mi-pO 

KKi-pO 

FS-^pO 

ZZ-«-pO 

ZZC-»-pO 

TT-t-p  0 

TTC-*-p  0 

Wf/V+1 

D*0 

B1-O 

AA-etN*    1+N*2 
a  COMPUTE  CUMULATIVE  DIST  AND  ASSOCIATED  STATS. 

PI :  P-*-PMATSPIn  -  5  :  ;  ] 

PC-<-(PC2:]*0)/PC2:] 
n         THIS  LOOP  CALCULATES  ALPHA   VALUES  AND  APPORXIMATIONS . 
L2:BfB+l 

Ci-ALPHAZBl 

A*-(  +  /lPCiC)) 

RBO-^PtlxAl 

KK->rKK,RHO 

PP*PP,PC2:i5] 
P  COMPUTE  NORMAL  APPROXIMATION . 

Z^RHO^UN-D^O.S) 

ZZi-ZZ,ll-NORMCDF   Z) 
a  COMPUTE   FORMAL   APPROXIMATION  WITH   CONTINUITI   CORRECTION. 

ZC-<-(RBO-AA  }x(  (N-l  >0.5) 

ZZC->-ZZC,.l--'JORMCDF    ZC  ) 
a  COMPUTE   STUDENT  T  APPROXIMATION 

TTi-TT.N  SPAPROXl    REO 
n     COMPUTE   STUDENT  T  APPROXIMATION  WITH   CONTINUITI   CORRECTION. 

TTCi-TTC.N   SPAPROX  RHO 

*(S<7)/i2 
0  COMPUTE  ERROR  DIFFERENCES. 

ERRZZ-PP-ZZ 

ERRZZC*PP-ZZC 

ERRTTi-PP-TT 

ERRTCi-PP-TTC 
p  PRINT  OUT  TABLE  OF   VALUES. 

1*0  . 1 

J-^19  5 

S-'-'TEST  STAT.    VALUE      • 

Mi-M,<  PROBZR   i  fi]:  FOR  SAMPLE  SIZE  EQUAL  TO    ',(2  0  wiV).* 

Tf2 


.881 
89] 


■*L2 
L5:J-<-19u 
L2:M*M.Ln    CU7C  (18pl97  )  .J'.Slp  (8pl97  )  .J'] 

■*{D  =  S)/L^ 
H*   1  18  dS 

F*'i(:uC<T  9.9999  >'  UFMT(.l    1    pKK) 

F*   1    62    pF 

FS-^E.F 

Lo  :FS-»-oO 
3->-l7oS 
FS-^FS,3 

L5:C<-C+1 

F-^'    I'  ,(7  5  9PPZC1  ) 

FS-*-FS.F 

♦(C<7)/L9 
Ll-.M-H. ill    FS 

D-»-D  +  l 

1^1 

N\:S*-^EXACT  C.D.F.  ' 

.7*  1 9  a 

*l3 
N2:S-' ERROR;    NORMAL 

PPfERRZZ 

*L3 
N2:S*-' ERROR;    NORM.    W/CC 

PP-*-ERRZZC 

N^1s->-' ERROR;    T  DIST 
^'°f-ERR'^'" 

Nl-.S*- ERROR;    T  W/CC 

PP-ERRTC 

■*L2 
L^:M 

-UV<10)/L1 

1 


I 

.11. 

,12 

13 

,15 

.16. 

17 

18 

19 

.20. 

.21 

.22, 
23 

.2**. 
25 
26 

,27. 
28 
29 

.30 
1 

!32; 
\^K 

35 

.36. 

,37, 
38 
39 

.'^O 

,'*2, 
>5 


7  WILTEST;NiALPEAxB:A;C'.T;P:NUM:DEN'.Z;ZZiERRZZ;MiZZC;D;PP;NUMC;I;B;F;FSiS 
;J;TTiT;TTC;XC;NUMCiNUMIiDENOMiDENOMC;ERRAVE;ERRTT;ERRTC;ERRZZC;ZC;ERRT 
CZC;KK 

TEIS  PROGRAM  GENERATES  TABLES  OF  C.D.F.    COMPARISONS  FOR  THE  WILCOXON 
SIGNED -RANK  TEST. 


SET  SAMPLE  SIZE  AND  ALPHA   V ALOES. 

01  0.02  0.03  0.05  0.08  0.13  0.18 
TEIS   LOOP  INCREMENTS   SAMPLE  SIZE. 


COMPUTE   CUMULATIVE  DIST  AND  ASSOCIATED   STATS. 

TEIS  LOOP  CALCULATES  ALPEA   VALUES  AND  APPORXIMATIONS . 


N--S 

ALPEAi-   0. 

Ll:PP-^pO 
M-^pO 
KK->-pO 
?S*pO 
Z2-«-p0 
ZZC-^pO 

rr-t-po 

TTC-*-p  0 
N'>-N+l 
D-^0 
Bt-0 

P-^WILP  N 

L2:B*B+1 

CfALPBA LBl 

.?-<-( +/(PiC)) 

T*A-1 

RK*KK,T 

PP-^PP.PiAl 
0  COMPUTE  NORMAL   APPROXIMATION . 

NUM-^T-CiNxiN+D^i-u) 

DEN-{  iNy■^N+l)x(  C2x/V)  +  i;)-i-2'i)*0.5 

Z'>-iWMi-DEN 

ZZ->-ZZ  .NOP.MCDF   Z 
0  COMPUTE  NORMAL   APPROXIMATION  WITE   CONTINUITI   CORRECTION. 

A/raC*(r+0  .  5  )  -  ( (iVx  (W+i )  )  +  u  ) 

ZC-^NUMC^DEN 

ZZC-rZZC  .NORMCDF   ZC 
P  COMPUTE   STUDENT  T  APPROXIMATION 

DENOMi-((  («x  {DEN*2  )  )+  {N-1  ))-(,  (.NUM*2  )+  (W-1  )  )  )*0  .  5 

T^NUMiDENOM 

TT*TT ,(.N-1)    TDISTN  T 
0     COMPUTE  STUDENT  T  APPROXIMATION  WITE   CONTINUITI   CORRECTION . 

NUMT-^{NUM)*0.b 

DENOMC'>rU(Nx(DENi'2))*(N-l))-(.(.(.iNUM)  +  0.5)i'2)*(.N-l)))*0.5 

TC->-NUMT*DENOMC 

TTC*TTC.UN-1)   TDISTN  TO 


113 


*(5<7 )/L2 
p  COMPUTE  ERROR  DIFFERENCES. 

ERRZZ'>rPP-ZZ 

EREZZC-PP-ZZC 

ERRT  T'-i-PP  -  'PT 

ERRTC-i-PP-ITC 

ERRAVE'>-PP-{iZZ+rT)*2)       ^ 

ERRTCZC*PP-l{ZZC+TTC)*2) 
p  PRINT  OUT  TABLE  OF   VALUES. 

1*0.1 

S^'TEST   STAT.    VALUE      ' 

M^M , '  PROS LW  ^   Wit    FOR   SAMPLE  SIZE  EQUAL   TO 

^L3 
L5:Ji-19»- 


,(2    0    v>N),' 


L3:M-<-M.CJ]    CWVC  (18pl97  ) ,  J-,61p  (8pl97  ) ,  J] 
(Z?=8)/r 


>'    QFMTd    7    pKK) 


3    tPPICl ) 
FS 


B*    1    18    qS 

^{D~Q)/Lo 

F*'BG<\       ZZ% 

F-«-    1    62    pF 

FS-R,F 

*L1 
Le:FS*pO 

E-rllpS 

FS^FS.E 

Cf-0 
L9:C*C+1 

F*'     I' .(7 

FS-i-FS.F 

*i.C<7)/L9 
L7 :M*M. :i] 

D-rD  +  1 

1*1 

->-(.Nl  ,N2  .N3  ,N^,N5  .N6  ,N7  ,15  )  IDl 
Nl:S*^EXACT  C.D.F.  ' 

J*198 

N2:S*' ERROR;    NORMAL 

PP*ERRZZ 

■*L3 
N3:S*' ERROR;    NORM.    W/CC 

PP*ERRZZC 

*L3 
Nu:S*' ERROR;    T  VIST 

PP*ERRTT 

*L3 
N5:S*' ERROR;    T  W/CC 

PP*ERRTC 

*L3 
N6:S*' ERROR;    AVE  T/Z 

PP*ERRAVE 

^L3 

N7'.S*' ERROR;    AVE  TC/ZC    ' 
PPi-ERRTCZC 
*L3 

Ln:M 
*(.N<10)/L1 

V 
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